Nucleic Acids Research - current issue

Nucleic Acids Research: VOLUME 37 ISSUE 11 2009

2009-06-21

Nucleic Acids Research

2009-06-21

Subscriptions

2009-06-21

TARGeT: a web-based pipeline for retrieving and characterizing gene and transposable element families from genomic sequences

Han, Y., Burnette, J. M., Wessler, S. R. — 2009-06-21

Gene families compose a large proportion of eukaryotic genomes. The rapidly expanding genomic sequence database provides a good opportunity to study gene family evolution and function. However, most gene family identification programs are restricted to searching protein databases where data are often lagging behind the genomic sequence data. Here, we report a user-friendly web-based pipeline, named TARGeT (Tree Analysis of Related Genes and Transposons), which uses either a DNA or amino acid ‘seed’ query to: (i) automatically identify and retrieve gene family homologs from a genomic database, (ii) characterize gene structure and (iii) perform phylogenetic analysis. Due to its high speed, TARGeT is also able to characterize very large gene families, including transposable elements (TEs). We evaluated TARGeT using well-annotated datasets, including the ascorbate peroxidase gene family of rice, maize and sorghum and several TE families in rice. In all cases, TARGeT rapidly recapitulated the known homologs and predicted new ones. We also demonstrated that TARGeT outperforms similar pipelines and has functionality that is not offered elsewhere.

Text-based over-representation analysis of microarray gene lists with annotation bias

Leong, H. S., Kipling, D. — 2009-06-21

A major challenge in microarray data analysis is the functional interpretation of gene lists. A common approach to address this is over-representation analysis (ORA), which uses the hypergeometric test (or its variants) to evaluate whether a particular functionally defined group of genes is represented more than expected by chance within a gene list. Existing applications of ORA have been largely limited to pre-defined terminologies such as GO and KEGG. We report our explorations of whether ORA can be applied to a wider mining of free-text. We found that a hitherto underappreciated feature of experimentally derived gene lists is that the constituents have substantially more annotation associated with them, as they have been researched upon for a longer period of time. This bias, a result of patterns of research activity within the biomedical community, is a major problem for classical hypergeometric test-based ORA approaches, which cannot account for such bias. We have therefore developed three approaches to overcome this bias, and demonstrate their usability in a wide range of published datasets covering different species. A comparison with existing tools that use GO terms suggests that mining PubMed abstracts can reveal additional biological insight that may not be possible by mining pre-defined ontologies alone.

Sim4cc: a cross-species spliced alignment program

Zhou, L., Pertea, M., Delcher, A. L., Florea, L. — 2009-06-21

Advances in sequencing technologies have accelerated the sequencing of new genomes, far outpacing the generation of gene and protein resources needed to annotate them. Direct comparison and alignment of existing cDNA sequences from a related species is an effective and readily available means to determine genes in the new genomes. Current spliced alignment programs are inadequate for comparing sequences between different species, owing to their low sensitivity and splice junction accuracy. A new spliced alignment tool, sim4cc, overcomes problems in the earlier tools by incorporating three new features: universal spaced seeds, to increase sensitivity and allow comparisons between species at various evolutionary distances, and powerful splice signal models and evolutionarily-aware alignment techniques, to improve the accuracy of gene models. When tested on vertebrate comparisons at diverse evolutionary distances, sim4cc had significantly higher sensitivity compared to existing alignment programs, more than 10% higher than the closest competitor for some comparisons, while being comparable in speed to its predecessor, sim4. Sim4cc can be used in one-to-one or one-to-many comparisons of genomic and cDNA sequences, and can also be effectively incorporated into a high-throughput annotation engine, as demonstrated by the mapping of 64 000 Fagus grandifolia 454 ESTs and unigenes to the poplar genome.

Robust methods for purification of histones from cultured mammalian cells with the preservation of their native modifications

Rodriguez-Collazo, P., Leuba, S. H., Zlatanova, J. — 2009-06-21

Post-translational modifications (PTMs) of histones play a role in modifying chromatin structure for DNA-templated processes in the eukaryotic nucleus, such as transcription, replication, recombination and repair; thus, histone PTMs are considered major players in the epigenetic control of these processes. Linking specific histone PTMs to gene expression is an arduous task requiring large amounts of highly purified and natively modified histones to be analyzed by various techniques. We have developed robust and complementary procedures, which use strong protein denaturing conditions and yield highly purified core and linker histones from unsynchronized proliferating, M-phase arrested and butyrate-treated cells, fully preserving their native PTMs without using enzyme inhibitors. Cell hypotonic swelling and lysis, nuclei isolation/washing and chromatin solubilization under mild conditions are bypassed to avoid compromising the integrity of histone native PTMs. As controls for our procedures, we tested the most widely used conventional methodologies and demonstrated that they indeed lead to drastic histone dephosphorylation. Additionally, we have developed methods for preserving acid-labile histone modifications by performing non-acid extractions to obtain highly purified H3 and H4. Importantly, isolation of histones H3, H4 and H2A/H2B is achieved without the use of HPLC. Functional supercoiling assays reveal that both hyper- and hypo-phosphorylated histones can be efficiently assembled into polynucleosomes. Notably, the preservation of fully phosphorylated mitotic histones and their assembly into polynucleosomes should open new avenues to investigate an important but overlooked question: the impact of mitotic phosphorylation in chromatin structure and function.

Tissue-specific regulatory network extractor (TS-REX): a database and software resource for the tissue and cell type-specific investigation of transcription factor-gene networks

2009-06-21

The prediction of transcription factor binding sites in genomic sequences is in principle very useful to identify upstream regulatory factors. However, when applying this concept to genomes of multicellular organisms such as mammals, one has to deal with a large number of false positive predictions since many transcription factor genes are only expressed in specific tissues or cell types. We developed TS-REX, a database/software system that supports the analysis of tissue and cell type-specific transcription factor-gene networks based on expressed sequence tag abundance of transcription factor-encoding genes in UniGene EST libraries. The use of expression levels of transcription factor-encoding genes according to hierarchical anatomical classifications covering different tissues and cell types makes it possible to filter out irrelevant binding site predictions and to identify candidates of potential functional importance for further experimental testing. TS-REX covers ESTs from H. sapiens and M. musculus, and allows the characterization of both presence and specificity of transcription factors in user-specified tissues or cell types. The software allows users to interactively visualize transcription factor-gene networks, as well as to export data for further processing. TS-REX was applied to predict regulators of Polycomb group genes in six human tumor tissues and in human embryonic stem cells.

MM-align: a quick algorithm for aligning multiple-chain protein complex structures using iterative dynamic programming

Mukherjee, S., Zhang, Y. — 2009-06-21

Structural comparison of multiple-chain protein complexes is essential in many studies of protein–protein interactions. We develop a new algorithm, MM-align, for sequence-independent alignment of protein complex structures. The algorithm is built on a heuristic iteration of a modified Needleman–Wunsch dynamic programming (DP) algorithm, with the alignment score specified by the inter-complex residue distances. The multiple chains in each complex are first joined, in every possible order, and then simultaneously aligned with cross-chain alignments prevented. The alignments of interface residues are enhanced by an interface-specific weighting factor. MM-align is tested on a large-scale benchmark set of 205 x 3897 non-homologous multiple-chain complex pairs. Compared with a naïve extension of the monomer alignment program of TM-align, the alignment accuracy of MM-align is significantly higher as judged by the average TM-score of the physically-aligned residues. MM-align is about two times faster than TM-align because of omitting the cross-alignment zone of the DP matrix. It also shows that the enhanced alignment of the interfaces helps in identifying biologically relevant protein complex pairs.

Replication fork reversal and the maintenance of genome stability

Atkinson, J., McGlynn, P. — 2009-06-21

The progress of replication forks is often threatened in vivo, both by DNA damage and by proteins bound to the template. Blocked forks must somehow be restarted, and the original blockage cleared, in order to complete genome duplication, implying that blocked fork processing may be critical for genome stability. One possible pathway that might allow processing and restart of blocked forks, replication fork reversal, involves the unwinding of blocked forks to form four-stranded structures resembling Holliday junctions. This concept has gained increasing popularity recently based on the ability of such processing to explain many genetic observations, the detection of unwound fork structures in vivo and the identification of enzymes that have the capacity to catalyse fork regression in vitro. Here, we discuss the contexts in which fork regression might occur, the factors that may promote such a reaction and the possible roles of replication fork unwinding in normal DNA metabolism.

Nontarget DNA binding shapes the dynamic landscape for enzymatic recognition of DNA damage

Friedman, J. I., Majumdar, A., Stivers, J. T. — 2009-06-21

The DNA repair enzyme human uracil DNA glycosylase (UNG) scans short stretches of genomic DNA and captures rare uracil bases as they transiently emerge from the DNA duplex via spontaneous base pair breathing motions. The process of DNA scanning requires that the enzyme transiently loosen its grip on DNA to allow stochastic movement along the DNA contour, while engaging extrahelical bases requires motions on a more rapid timescale. Here, we use NMR dynamic measurements to show that free UNG has no intrinsic dynamic properties in the millisecond to microsecond and subnanosecond time regimes, and that the act of binding to nontarget DNA reshapes the dynamic landscape to allow productive millisecond motions for scanning and damage recognition. These results suggest that DNA structure and the spontaneous dynamics of base pairs may drive the evolution of a protein sequence that is tuned to respond to this dynamic regime.

Poly(ADP-ribosyl)ation of heterogeneous nuclear ribonucleoproteins modulates splicing

Ji, Y., Tulin, A. V. — 2009-06-21

The biological functions of poly(ADP-ribosyl)ation of heterogeneous nuclear ribonucleoproteins (hnRNPs) are not well understood. However, it is known that hnRNPs are involved in the regulation of alternative splicing for many genes, including the Ddc gene in Drosophila. Therefore, we first confirmed that poly(ADP-ribose) (pADPr) interacts with two Drosophila hnRNPs, Squid/hrp40 and Hrb98DE/hrp38, and that this function is regulated by Poly(ADP-ribose) Polymerase 1 (PARP1) and Poly(ADP-ribose) Glycohydrolase (PARG) in vivo. These findings then provided a basis for analyzing the role of pADPr binding to these two hnRNPs in terms of alternative splicing regulation. Our results showed that Parg null mutation does cause poly(ADP-ribosyl)ation of Squid and hrp38 protein, as well as their dissociation from active chromatin. Our data also indicated that pADPr binding to hnRNPs inhibits the RNA-binding ability of hnRNPs. Following that, we demonstrated that poly(ADP-ribosyl)ation of Squid and hrp38 proteins inhibits splicing of the intron in the Hsr-RC transcript, but enhances splicing of the intron in the Ddc pre-mRNA. Taken together, these findings suggest that poly(ADP-ribosyl)ation regulates the interaction between hnRNPs and RNA and thus modulates the splicing pathways.

The amino terminal domain from Mrt4 protein can functionally replace the RNA binding domain of the ribosomal P0 protein

2009-06-21

In Saccharomyces cerevisiae, the Mrt4 protein is a component of the ribosome assembly machinery that shares notable sequence homology to the P0 ribosomal stalk protein. Here, we show that these proteins can not bind simultaneously to ribosomes and moreover, a chimera containing the first 137 amino acids of Mrt4 and the last 190 amino acids from P0 can partially complement the absence of the ribosomal protein in a conditional P0 null mutant. This chimera is associated with ribosomes isolated from this strain when grown under restrictive conditions, although its binding is weaker than that of P0. These ribosomes contain less P1 and P2 proteins, the other ribosomal stalk components. Similarly, the interaction of the L12 protein, a stalk base component, is affected by the presence of the chimera. These results indicate that Mrt4 and P0 bind to the same site in the 25S rRNA. Indeed, molecular dynamics simulations using modelled Mrt4 and P0 complexes provide further evidence that both proteins bind similarly to rRNA, although their interaction with L12 displays notable differences. Together, these data support the participation of the Mrt4 protein in the assembly of the P0 protein into the ribosome and probably, that also of the L12 protein.

PROCAIN: protein profile comparison with assisting information

Wang, Y., Sadreyev, R. I., Grishin, N. V. — 2009-06-21

Detection of remote sequence homology is essential for the accurate inference of protein structure, function and evolution. The most sensitive detection methods involve the comparison of evolutionary patterns reflected in multiple sequence alignments (MSAs) of protein families. We present PROCAIN, a new method for MSA comparison based on the combination of ‘vertical’ MSA context (substitution constraints at individual sequence positions) and ‘horizontal’ context (patterns of residue content at multiple positions). Based on a simple and tractable profile methodology and primitive measures for the similarity of horizontal MSA patterns, the method achieves the quality of homology detection comparable to a more complex advanced method employing hidden Markov models (HMMs) and secondary structure (SS) prediction. Adding SS information further improves PROCAIN performance beyond the capabilities of current state-of-the-art tools. The potential value of the method for structure/function predictions is illustrated by the detection of subtle homology between evolutionary distant yet structurally similar protein domains. ProCAIn, relevant databases and tools can be downloaded from: http://prodata.swmed.edu/procain/download. The web server can be accessed at http://prodata.swmed.edu/procain/procain.php.

Cleavage of a model DNA replication fork by a Type I restriction endonuclease

Ishikawa, K., Handa, N., Kobayashi, I. — 2009-06-21

Cleavage of a DNA replication fork leads to fork restoration by recombination repair. In prokaryote cells carrying restriction–modification systems, fork passage reduces genome methylation by the modification enzyme and exposes the chromosome to attack by the restriction enzyme. Various observations have suggested a relationship between the fork and Type I restriction enzymes, which cleave DNA at a distance from a recognition sequence. Here, we demonstrate that a Type I restriction enzyme preparation cleaves a model replication fork at its branch. The enzyme probably tracks along the DNA from an unmethylated recognition site on the daughter DNA and cuts the fork upon encountering the branch point. Our finding suggests that these restriction–modification systems contribute to genome maintenance through cell death and indicates that DNA replication fork cleavage represents a critical point in genome maintenance to choose between the restoration pathway and the destruction pathway.

A genotype-to-phenotype map of in vitro selected RNA-cleaving DNAzymes: implications for accessing the target phenotype

Schlosser, K., Lam, J. C. F., Li, Y. — 2009-06-21

Herein, we describe a case study into the population dynamics of in vitro selection, using RNA-cleaving DNAzymes as a model system. We sought to understand how the composition of the population can change over time in response to different levels of selection pressure, and how well these changes are correlated with selection of the target phenotype. The model population is composed of 857 DNAzyme clones representing 215 discrete sequence classes, which had previously been identified from two parallel selection experiments, conducted under an increasingly stringent, or permissive and constant selection time pressure. In this report, we determined the principal phenotypic properties (i.e. k_obs, maximum cleavage yield and PCR efficiency) from a sample of 58 clones representing 46 different major and minor sequence classes from various rounds of each selection experiment. Interestingly, a positive correlation between the catalytic rate constant and the corresponding frequency and temporal position of a given DNAzyme was not consistently observed; however, the strength of the correlation was qualitatively higher under conditions of more stringent selection time pressure. These results suggest that the selective sampling paradigm on which in vitro selection is based, may underestimate the true functional capacity of any given random-sequence library.

Non-specific interactions are sufficient to explain the position of heterochromatic chromocenters and nucleoli in interphase nuclei

de Nooijer, S., Wellink, J., Mulder, B., Bisseling, T. — 2009-06-21

The organization of the eukaryote nucleus into functional compartments arises by self-organization both through specific protein–protein and protein–DNA interactions and non-specific interactions that lead to entropic effects, such as e.g. depletion attraction. While many specific interactions have so far been demonstrated, the contributions of non-specific interactions are still unclear. We used coarse-grained molecular dynamics simulations of previously published models for Arabidopsis thaliana chromatin organization to show that non-specific interactions can explain the in vivo localization of nucleoli and chromocenters. Also, we quantitatively demonstrate that chromatin looping contributes to the formation of chromosome territories. Our results are consistent with the previously published Rosette model for Arabidopsis chromatin organization and suggest that chromocenter-associated loops play a role in suppressing chromocenter clustering.

Accurate prediction of NAGNAG alternative splicing

2009-06-21

Alternative splicing (AS) involving NAGNAG tandem acceptors is an evolutionarily widespread class of AS. Recent predictions of alternative acceptor usage reported better results for acceptors separated by larger distances, than for NAGNAGs. To improve the latter, we aimed at the use of Bayesian networks (BN), and extensive experimental validation of the predictions. Using carefully constructed training and test datasets, a balanced sensitivity and specificity of ≥92% was achieved. A BN trained on the combined dataset was then used to make predictions, and 81% (38/47) of the experimentally tested predictions were verified. Using a BN learned on human data on six other genomes, we show that while the performance for the vertebrate genomes matches that achieved on human data, there is a slight drop for Drosophila and worm. Lastly, using the prediction accuracy according to experimental validation, we estimate the number of yet undiscovered alternative NAGNAGs. State of the art classifiers can produce highly accurate prediction of AS at NAGNAGs, indicating that we have identified the major features of the ‘NAGNAG-splicing code’ within the splice site and its immediate neighborhood. Our results suggest that the mechanism behind NAGNAG AS is simple, stochastic, and conserved among vertebrates and beyond.

mwr Xer site-specific recombination is hypersensitive to DNA supercoiling

2009-06-21

The multiresistance plasmid pJHCMW1, first identified in a Klebsiella pneumoniae strain isolated from a neonate with meningitis, includes a Xer recombination site, mwr, with unique characteristics. Efficiency of resolution of mwr-containing plasmid dimers is strongly dependent on the osmotic pressure of the growth medium. An increase in supercoiling density of plasmid DNA was observed as the osmotic pressure of the growth culture decreased. Reporter plasmids containing directly repeated mwr, or the related cer sites were used to test if DNA topological changes were correlated with significant changes in efficiency of Xer recombination. Quantification of Holliday junctions showed that while recombination at cer was efficient at all levels of negative supercoiling, recombination at mwr became markedly less efficient as the level of supercoiling was reduced. These results support a model in which modifications at the level of supercoiling density caused by changes in the osmotic pressure of the culture medium affects resolution of mwr-containing plasmid dimers, a property that separates mwr from other Xer recombination target sites.

Protein-coding gene promoters in Methanocaldococcus (Methanococcus) jannaschii

Zhang, J., Li, E., Olsen, G. J. — 2009-06-21

Although Methanocaldococcus (Methanococcus) jannaschii was the first archaeon to have its genome sequenced, little is known about the promoters of its protein-coding genes. To expand our knowledge, we have experimentally identified 131 promoters for 107 protein-coding genes in this genome by mapping their transcription start sites. Compared to previously identified promoters, more than half of which are from genes for stable RNAs, the protein-coding gene promoters are qualitatively similar in overall sequence pattern, but statistically different at several positions due to greater variation among their sequences. Relative binding affinity for general transcription factors was measured for 12 of these promoters by competition electrophoretic mobility shift assays. These promoters bind the factors less tightly than do most tRNA gene promoters. When a position weight matrix (PWM) was constructed from the protein gene promoters, factor binding affinities correlated with corresponding promoter PWM scores. We show that the PWM based on our data more accurately predicts promoters in the genome and transcription start sites than could be done with the previously available data. We also introduce a PWM logo, which visually displays the implications of observing a given base at a position in a sequence.

Rapid Cdc13 turnover and telomere length homeostasis are controlled by Cdk1-mediated phosphorylation of Cdc13

Tseng, S.-F., Shen, Z.-J., Tsai, H.-J., Lin, Y.-H., Teng, S.-C. — 2009-06-21

Budding yeast telomerase is mainly activated by Tel1/Mec1 (yeast ATM/ATR) on Cdc13 from late S to G2 phase of the cell cycle. Here, we demonstrated that the telomerase-recruitment domain of Cdc13 is also phosphorylated by Cdk1 at the same cell cycle stage as the Tel1/Mec1-dependent regulation. Phosphor-specific gel analysis demonstrated that Cdk1 phosphorylates residues 308 and 336 of Cdc13. The residue T308 of Cdc13 is critical for efficient Mec1-mediated S306 phosphorylation in vitro. Phenotypic analysis in vivo revealed that the mutations in the Cdc13 S/TP motifs phosphorylated by Cdk1 caused cell cycle delay and telomere shortening and these phenotypes could be partially restored by the replacement with a negative charge residue. In the absence of Ku or Tel1, Cdk1-mediated phosphorylation of Cdc13 showed no effect on telomere length maintenance. Moreover, this Cdk1-mediated phosphorylation was required to promote the regular turnover of Cdc13. Together these results demonstrate that Cdk1 phosphorylates the telomerase recruitment domain of Cdc13, thereby preserves optimal function and expression level of Cdc13 for precise telomere replication and cell cycle progression.

Alternative polyadenylation variants of the RNA binding protein, HuR: abundance, role of AU-rich elements and auto-Regulation

Al-Ahmadi, W., Al-Ghamdi, M., Al-Haj, L., Al-Saif, M., Khabar, K. S. A. — 2009-06-21

The RNA-binding protein, HuR, is involved in the stabilization of AU-rich element-containing mRNAs with products that are involved in cell-cycle progression, cell differentiation and inflammation. We show that there are multiple polyadenylation variants of HuR mRNA that differ in their abundance, using both bioinformatics and experimental approaches. A polyadenylation variant with distal poly(A) signal is a rare transcript that harbors functional AU-rich elements (ARE) in the 3'UTR. A minimal 60-nt region, but not a mutant form, fused to reporter-3'UTR constructs was able to downregulate the reporter activity. The most predominant and alternatively polyadenylated mature transcript does not contain the ARE. HuR itself binds HuR mRNA, and upregulated the activity of reporter from constructs fused with ARE-isoform and the HuR ARE. Wild-type tristetraprolin (TTP), but not the zinc finger mutant TTP, competes for HuR binding and upregulation of HuR mRNA. The study shows that the HuR gene codes for several polyadenylation variants differentially regulated by AU-rich elements, and demonstrates an auto-regulatory role of HuR.

Dramatic effect of single-base mutation on the conformational dynamics of human telomeric G-quadruplex

Lee, J. Y., Kim, D. S. — 2009-06-21

Guanine-rich DNA sequences can form G-quadruplexes. These four-stranded structures are known to form in several genomic regions and to influence certain biological activities. Sometimes, the instability of G-quadruplexes causes the abnormal biological processes. Mutation is a culprit for the destabilization of G-quadruplexes, but the details of mutated G-quadruplexes are poorly understood. In this article, we investigated the conformational dynamics of single-base mutated human telomeric G-quadruplexes in the presence of K⁺ with single-molecule FRET spectroscopy. We observed that the replacement of single guanine by thymine in a G-track induces various folded structures, i.e. structural polymorphism. Moreover, direct observation of their dynamics revealed that a single-base mutation causes fast unfolding of folded states under physiological conditions. Furthermore, we found that the degree of destabilization varies according to mutation positions. When the central guanine of a G-track is replaced, the G-quadruplexes unfold quickly at any K⁺ concentrations and temperature. Meanwhile, outer-quartet mutated G-quadruplexes have heterogeneous dynamics at intermediate K⁺ concentrations and longstanding folded states at high K⁺ concentrations. Several factors such as base-stacking interaction and K⁺ coordination are responsible for the different dynamics according to the mutation position.

Targeted correction of a thalassemia-associated {beta}-globin mutation induced by pseudo-complementary peptide nucleic acids

2009-06-21

β-Thalassemia is a genetic disorder caused by mutations in the β-globin gene. Triplex-forming oligonucleotides and triplex-forming peptide nucleic acids (PNAs) have been shown to stimulate recombination in mammalian cells via site-specific binding and creation of altered helical structures that provoke DNA repair. However, the use of these molecules for gene targeting requires homopurine tracts to facilitate triple helix formation. Alternatively, to achieve binding to mixed-sequence target sites for the induced gene correction, we have used pseudo-complementary PNAs (pcPNAs). Due to steric hindrance, pcPNAs are unable to form pcPNA–pcPNA duplexes but can bind to complementary DNA sequences via double duplex-invasion complexes. We demonstrate here that pcPNAs, when co-transfected with donor DNA fragments, can promote single base pair modification at the start of the second intron of the beta-globin gene. This was detected by the restoration of proper splicing of transcripts produced from a green fluorescent protein-beta globin fusion gene. We also demonstrate that pcPNAs are effective in stimulating recombination in human fibroblast cells in a manner dependent on the nucleotide excision repair factor, XPA. These results suggest that pcPNAs can be effective tools to induce heritable, site-specific modification of disease-related genes in human cells without purine sequence restriction.

Coupled termination/reinitiation for translation of the downstream open reading frame B of the prototypic hypovirus CHV1-EP713

Guo, L.-h., Sun, L., Chiba, S., Araki, H., Suzuki, N. — 2009-06-21

Cryphonectria hypovirus 1 (CHV1), associated with the picorna-like superfamily, infects the chestnut blight fungus and attenuates the virulence of the host fungus. The genomic RNA of the virus has two continuous open reading frames, A and B, separated by the pentanucleotide UAAUG. We present here evidence suggesting that ORF B is translated from genome-sized virus mRNA by a coupled termination/reinitiation mechanism mediated by the pentamer. In the coupled translation, the overlapping UAA and AUG triplets serve as the stop codon of ORF A and the initiator of ORF B, respectively. This was established by the use of a luciferase assay with a basic construct containing the ORF A sequence and the firefly luciferase gene while retaining the pentamer between the two coding sequences. The proportion of ribosomes reinitiating translation after terminating was determined to be 2.5–4.4% by three independent assay systems in fungal and insect cells. Use of a series of mutant constructs identified two sequence elements, the pentamer and the p40 sequence, that affect the efficiency of coupled translation and virus replication. Together, these results provide the first example of coupled translation facilitated by the pentanucleotide UAAUG in the kingdom Fungi. The mechanism by which the preceding p40-coding sequence promotes reinitiation is discussed.

Mutational analysis of the HIV-1 auxiliary protein Vif identifies independent domains important for the physical and functional interaction with HIV-1 reverse transcriptase

2009-06-21

The HIV-1 accessory protein Vif plays a dual role: it counteracts the natural restriction factors APOBEC3G and 3F and ensures efficient retrotranscription of the HIV-1 RNA genome. We have previously shown that Vif can act as an auxiliary factor for HIV-1 reverse transcriptase (RT), increasing its rate of association to RNA or DNA templates. Here, by using seven different Vif mutants, we provide in vitro evidences that Vif stimulates HIV-1 RT through direct protein–protein interaction, which is mediated by its C-terminal domain. Physical interaction appears to require the proline-rich region comprised between amino acid (aa) 161 and 164 of Vif, whereas the RT stimulatory activity requires, in addition, the extreme C-terminal region (aa 169–192) of the Vif protein. Neither the RNA interaction domain, nor the Zn⁺⁺-binding domain of Vif are required for its interaction with the viral RT. Pseudotyped HIV-1 lentiviral vectors bearing Vif mutants deleted in the RNA- or RT-binding domains show defects in retrotranscription/integration processes in both permissive and nonpermissive cells. Our results broaden our knowledge on how three important functions of Vif (RNA binding, RT binding and stimulation and Zn⁺⁺ binding), are coordinated by different domains.

Constructing RNA dynamical ensembles by combining MD and motionally decoupled NMR RDCs: new insights into RNA dynamics and adaptive ligand recognition

Frank, A. T., Stelzer, A. C., Al-Hashimi, H. M., Andricioaei, I. — 2009-06-21

We describe a strategy for constructing atomic resolution dynamical ensembles of RNA molecules, spanning up to millisecond timescales, that combines molecular dynamics (MD) simulations with NMR residual dipolar couplings (RDC) measured in elongated RNA. The ensembles are generated via a Monte Carlo procedure by selecting snap-shot from an MD trajectory that reproduce experimentally measured RDCs. Using this approach, we construct ensembles for two variants of the transactivation response element (TAR) containing three (HIV-1) and two (HIV-2) nucleotide bulges. The HIV-1 TAR ensemble reveals significant mobility in bulge residues C24 and U25 and to a lesser extent U23 and neighboring helical residue A22 that give rise to large amplitude spatially correlated twisting and bending helical motions. Omission of bulge residue C24 in HIV-2 TAR leads to a significant reduction in both the local mobility in and around the bulge and amplitude of inter-helical bending motions. In contrast, twisting motions of the helices remain comparable in amplitude to HIV-1 TAR and spatial correlations between them increase significantly. Comparison of the HIV-1 TAR dynamical ensemble and ligand bound TAR conformations reveals that several features of the binding pocket and global conformation are dynamically preformed, providing support for adaptive recognition via a ‘conformational selection’ type mechanism.

Transcriptional regulation shapes the organization of genes on bacterial chromosomes

Janga, S. C., Salgado, H., Martinez-Antonio, A. — 2009-06-21

Transcription factors (TFs) are the key elements responsible for controlling the expression of genes in bacterial genomes and when visualized on a genomic scale form a dense network of transcriptional interactions among themselves and with other protein coding genes. Although the structure of transcriptional regulatory networks (TRNs) is well understood, it is not clear what constrains govern them. Here, we explore this question using the TRNs of model prokaryotes and provide a link between the transcriptional hierarchy of regulons and their genome organization. We show that, to drive the kinetics and concentration gradients, TFs belonging to big and small regulons, depending on the number of genes they regulate, organize themselves differently on the genome with respect to their targets. We then propose a conceptual model that can explain how the hierarchical structure of TRNs might be ultimately governed by the dynamic biophysical requirements for targeting DNA-binding sites by TFs. Our results suggest that the main parameters defining the position of a TF in the network hierarchy are the number and chromosomal distances of the genes they regulate and their protein concentration gradients. These observations give insights into how the hierarchical structure of transcriptional networks can be encoded on the chromosome to drive the kinetics and concentration gradients of TFs depending on the number of genes they regulate and could be a common theme valid for other prokaryotes, proposing the role of transcriptional regulation in shaping the organization of genes on a chromosome.

The C. elegans Snail homolog CES-1 can activate gene expression in vivo and share targets with bHLH transcription factors

2009-06-21

Snail-type transcription factors (TFs) are found in numerous metazoan organisms and function in a plethora of cellular and developmental processes including mesoderm and neuronal development, apoptosis and cancer. So far, Snail-type TFs are exclusively known as transcriptional repressors. They repress gene expression by recruiting transcriptional co-repressors and/or by preventing DNA binding of activators from the basic helix-loop-helix (bHLH) family of TFs to CAGGTG E-box sequences. Here we report that the Caenorhabditis elegans Snail-type TF CES-1 can activate transcription in vivo. Moreover, we provide results that suggest that CES-1 can share its binding site with bHLH TFs, in different tissues, rather than only occluding bHLH DNA binding. Together, our data indicate that there are at least two types of CES-1 target genes and, therefore, that the molecular function of Snail-type TFs is more plastic than previously appreciated.

Histone deacetylase Rpd3 antagonizes Sir2-dependent silent chromatin propagation

Zhou, J., Zhou, B. O., Lenzmeier, B. A., Zhou, J.-Q. — 2009-06-21

In the eukaryotic genome, transcriptionally silent chromatin tends to propagate along a chromosome and encroach upon adjacent active chromatin. The silencing machinery can be stopped by chromatin boundary elements. We performed a screen in Saccharomyces cerevisiae for proteins that may contribute to the establishment of a chromatin boundary. We found that disruption of histone deacetylase Rpd3p results in defective boundary activity, leading to a Sir-dependent local propagation of transcriptional repression. In rpd3 cells, the amount of Sir2p that was normally found in the nucleolus decreased and the amount of Sir2p found at telomeres and at HM and its adjacent loci increased, leading to an extension of silent chromatin in those areas. In addition, Rpd3p interacted directly with chromatin at boundary regions to deacetylate histone H4 at lysine 5 and at lysine 12. Either the mutation of histone H4 at lysine 5 or a decrease in the histone acetyltransferase (HAT) activity of Esa1p abrogated the silencing phenotype associated with rpd3 mutation, suggesting a novel role for the H4 amino terminus in Rpd3p-mediated heterochromatin boundary regulation. Together, these data provide insight into the molecular mechanisms for the anti-silencing functions of Rpd3p during the formation of heterochromatin boundaries.

Nucleotide analogs and molecular modeling studies reveal key interactions involved in substrate recognition by the yeast RNA triphosphatase

Issur, M., Despins, S., Bougie, I., Bisaillon, M. — 2009-06-21

RNA triphosphatases (RTPases) are involved in the addition of the distinctive cap structure found at the 5' ends of eukaryotic mRNAs. Fungi, protozoa and some DNA viruses possess an RTPase that belongs to the triphosphate tunnel metalloenzyme family of enzymes that can also hydrolyze nucleoside triphosphates. Previous crystallization studies revealed that the phosphohydrolase catalytic core is located in a hydrophilic tunnel composed of antiparallel β-strands. However, all past efforts to obtain structural information on the interaction between RTPases and their substrates were unsuccessful. In the present study, we used computational molecular docking to model the binding of a nucleotide substrate into the yeast RTPase active site. In order to confirm the docking model and to gain additional insights into the molecular determinants involved in substrate recognition, we also evaluated both the phosphohydrolysis and the inhibitory potential of an important number of nucleotide analogs. Our study highlights the importance of specific amino acids for the binding of the sugar, base and triphosphate moieties of the nucleotide substrate, and reveals both the structural flexibility and complexity of the active site. These data illustrate the functional features required for the interaction of an RTPase with a ligand and pave the way to the use of nucleotide analogs as potential inhibitors of RTPases of pathogenic importance.

Molecular mechanism of poly(ADP-ribosyl)ation by PARP1 and identification of lysine residues as ADP-ribose acceptor sites

Altmeyer, M., Messner, S., Hassa, P. O., Fey, M., Hottiger, M. O. — 2009-06-21

Poly(ADP-ribose) polymerase 1 (PARP1) synthesizes poly(ADP-ribose) (PAR) using nicotinamide adenine dinucleotide (NAD) as a substrate. Despite intensive research on the cellular functions of PARP1, the molecular mechanism of PAR formation has not been comprehensively understood. In this study, we elucidate the molecular mechanisms of poly(ADP-ribosyl)ation and identify PAR acceptor sites. Generation of different chimera proteins revealed that the amino-terminal domains of PARP1, PARP2 and PARP3 cooperate tightly with their corresponding catalytic domains. The DNA-dependent interaction between the amino-terminal DNA-binding domain and the catalytic domain of PARP1 increased V_max and decreased the K_m for NAD. Furthermore, we show that glutamic acid residues in the auto-modification domain of PARP1 are not required for PAR formation. Instead, we identify individual lysine residues as acceptor sites for ADP-ribosylation. Together, our findings provide novel mechanistic insights into PAR synthesis with significant relevance for the different biological functions of PARP family members.

Silencing of a plant gene by transcriptional interference

Hedtke, B., Grimm, B. — 2009-06-21

Integration of foreign DNA into eukaryotic genomes results frequently in a total or partial loss of gene function, caused by the interruption of indispensable structures of the gene itself. Using T-DNA insertions in Arabidopsis we screened for mutants with deregulated chlorophyll precursor accumulation in etiolated seedlings. A mutant designated rfd1 (red fluorescent in darkness) with increased protochlorophyllide accumulation showed a fluorescent phenotype that was associated with a lack of transcript initiation from the AtRibA1 promoter situated downstream of the integrated T-DNA. Complementation experiments confirmed rfd1 to be a knockout phenotype. Comparison with two SALK insertion lines bearing T-DNA integrations in the 5'UTR of AtRibA1 demonstrated that the insertion event in rfd1 itself does not explain the complete lack of transcript initiation. A 35S tetrameric enhancer sequence present on the rfd1 T-DNA causes the overaccumulation of a large polycistronic transcript originating inside the T-DNA. This 5.5-kb RNA runs over the downstream situated AtRibA1 promoter, which was shown by 5'RACE analyses to be consequently silenced. Hence, a transcription process that starts upstream and overlaps AtRibA1 blocks the initiation at the AtRibA1 promoter in rfd1. This regulatory mechanism has recently been introduced in yeast as transcriptional interference and is described here for the first time in a plant system.

Low modularity of aminoacyl-tRNA substrates in polymerization by the ribosome

Forster, A. C. — 2009-06-21

Aminoacyl-transfer RNAs contain four standardized units: amino acids, an invariant 3'-terminal CCA, trinucleotide anticodons and tRNA bodies. The degree of interchangeability of the three variable modules is poorly understood, despite its role in evolution and the engineering of translation to incorporate unnatural amino acids. Here, a purified translation system is used to investigate effects of various module swaps on the efficiency of multiple ribosomal incorporations of unnatural aminoacyl-tRNA substrates per peptide product. The yields of products containing three to five adjacent l-amino acids with unnatural side chains are low and cannot be improved by optimization or explained simply by any single factor tested. Though combinations of modules that allow quantitative single unnatural incorporations are found readily, finding combinations that enable efficient synthesis of products containing multiple unnatural amino acids is challenging. This implies that assaying multiple, as opposed to single, incorporations per product is a more stringent assay of substrate activity. The unpredictability of most results illustrates the multifactorial nature of substrate recognition and the value of synthetic biology for testing our understanding of translation. Data indicate that the degree of interchangeability of the modules of aminoacyl-tRNAs is low.

Kinetics and thermodynamics of DNA hybridization on gold nanoparticles

Chen, C., Wang, W., Ge, J., Zhao, X. S. — 2009-06-21

Hybridization of single-stranded DNA immobilized on the surface of gold nanoparticles (GNPs) into double stranded DNA and its subsequent dissociation into ssDNA were investigated. Melting curves and rates of dissociation and hybridization were measured using fluorescence detection based on hybridization-induced fluorescence change. Two distribution functions, namely the state distribution and the rate distribution, were proposed in order to take interfacial heterogeneity into account and to quantitatively analyze the data. Reaction and activation enthalpies and entropies of DNA hybridization and dissociation on GNPs were derived and compared with the same quantities in solution. Our results show that the interaction between GNPs and DNA reduces the energetic barrier and accelerates the dissociation of adhered DNA. At low surface densities of ssDNA adhered to GNP surface, the primary reaction pathway is that ssDNA in solution first adsorbs onto the GNP, and then diffuses along the surface until hybridizing with an immobilized DNA. We also found that the secondary structure of a DNA hairpin inhibits the interaction between GNPs and DNA and enhances the stability of the DNA hairpin adhered to GNPs.

Local and global effects of strong DNA bending induced during molecular dynamics simulations

Curuksu, J., Zacharias, M., Lavery, R., Zakrzewska, K. — 2009-06-21

DNA bending plays an important role in many biological processes, but its molecular and energetic details as a function of base sequence remain to be fully understood. Using a recently developed restraint, we have studied the controlled bending of four different B-DNA oligomers using molecular dynamics simulations. Umbrella sampling with the AMBER program and the recent parmbsc0 force field yield free energy curves for bending. Bending 15-base pair oligomers by 90° requires roughly 5 kcal mol^–1, while reaching 150° requires of the order of 12 kcal mol^–1. Moderate bending occurs mainly through coupled base pair step rolls. Strong bending generally leads to local kinks. The kinks we observe all involve two consecutive base pair steps, with disruption of the central base pair (termed Type II kinks in earlier work). A detailed analysis of each oligomer shows that the free energy of bending only varies quadratically with the bending angle for moderate bending. Beyond this point, in agreement with recent experiments, the variation becomes linear. An harmonic analysis of each base step yields force constants that not only vary with sequence, but also with the degree of bending. Both these observations suggest that DNA is mechanically more complex than simple elastic rod models would imply.

Low-fidelity DNA synthesis by the L979F mutator derivative of Saccharomyces cerevisiae DNA polymerase {zeta}

2009-06-21

To probe Pol functions in vivo via its error signature, here we report the properties of Saccharomyces cerevisiae Pol in which phenyalanine was substituted for the conserved Leu-979 in the catalytic (Rev3) subunit. We show that purified L979F Pol is 30% as active as wild-type Pol when replicating undamaged DNA. L979F Pol shares with wild-type Pol the ability to perform moderately processive DNA synthesis. When copying undamaged DNA, L979F Pol is error-prone compared to wild-type Pol , providing a biochemical rationale for the observed mutator phenotype of rev3-L979F yeast strains. Errors generated by L979F Pol in vitro include single-base insertions, deletions and substitutions, with the highest error rates involving stable misincorporation of dAMP and dGMP. L979F Pol also generates multiple errors in close proximity to each other. The frequency of these events far exceeds that expected for independent single changes, indicating that the first error increases the probability of additional errors within 10 nucleotides. Thus L979F Pol , and perhaps wild-type Pol , which also generates clustered mutations at a lower but significant rate, performs short patches of processive, error-prone DNA synthesis. This may explain the origin of some multiple clustered mutations observed in vivo.

Molecular genetics and comparative genomics reveal RNAi is not functional in malaria parasites

2009-06-21

Techniques for targeted genetic disruption in Plasmodium, the causative agent of malaria, are currently intractable for those genes that are essential for blood stage development. The ability to use RNA interference (RNAi) to silence gene expression would provide a powerful means to gain valuable insight into the pathogenic blood stages but its functionality in Plasmodium remains controversial. Here we have used various RNA-based gene silencing approaches to test the utility of RNAi in malaria parasites and have undertaken an extensive comparative genomics search using profile hidden Markov models to clarify whether RNAi machinery exists in malaria. These investigative approaches revealed that Plasmodium lacks the enzymology required for RNAi-based ablation of gene expression and indeed no experimental evidence for RNAi was observed. In its absence, the most likely explanations for previously reported RNAi-mediated knockdown are either the general toxicity of introduced RNA (with global down-regulation of gene expression) or a specific antisense effect mechanistically distinct from RNAi, which will need systematic analysis if it is to be of use as a molecular genetic tool for malaria parasites.

Crystal structure of the {beta}{beta}{alpha}-Me type II restriction endonuclease Hpy99I with target DNA

Sokolowska, M., Czapinska, H., Bochtler, M. — 2009-06-21

The ββ-Me restriction endonuclease (REase) Hpy99I recognizes the CGWCG target sequence and cleaves it with unusual stagger (five nucleotide 5'-recessed ends). Here we present the crystal structure of the specific complex of the dimeric enzyme with DNA. The Hpy99I protomer consists of an antiparallel β-barrel and two β42 repeats. Each repeat coordinates a structural zinc ion with four cysteine thiolates in two CXXC motifs. The ββ-Me region of the second β42 repeat holds the catalytic metal ion (or its sodium surrogate) via Asp148 and Asn165 and activates a water molecule with the general base His149. In the specific complex, Hpy99I forms a ring-like structure around the DNA that contacts DNA bases on the major and minor groove sides via the first and second β42 repeats, respectively. Hpy99I interacts with the central base pair of the recognition sequence only on the minor groove side, where A:T resembles T:A and G:C is similar to C:G. The Hpy99I–DNA co-crystal structure provides the first detailed illustration of the ββ-Me site in REases and complements structural information on the use of this active site motif in other groups of endonucleases such as homing endonucleases (e.g. I-PpoI) and Holliday junction resolvases (e.g. T4 endonuclease VII).

Transcriptionally active TFIIH of the early-diverged eukaryote Trypanosoma brucei harbors two novel core subunits but not a cyclin-activating kinase complex

Lee, J. H., Jung, H. S., Gunzl, A. — 2009-06-21

Trypanosoma brucei is a member of the early-diverged, protistan family Trypanosomatidae and a lethal parasite causing African Sleeping Sickness in humans. Recent studies revealed that T. brucei harbors extremely divergent orthologues of the general transcription factors TBP, TFIIA, TFIIB and TFIIH and showed that these factors are essential for initiating RNA polymerase II-mediated synthesis of spliced leader (SL) RNA, a trans splicing substrate and key molecule in trypanosome mRNA maturation. In yeast and metazoans, TFIIH is composed of a core of seven conserved subunits and the ternary cyclin-activating kinase (CAK) complex. Conversely, only four TFIIH subunits have been identified in T. brucei. Here, we characterize the first protistan TFIIH which was purified in its transcriptionally active form from T. brucei extracts. The complex consisted of all seven core subunits but lacked the CAK sub-complex; instead it contained two trypanosomatid-specific subunits, which were indispensable for parasite viability and SL RNA gene transcription. These findings were corroborated by comparing the molecular structures of trypanosome and human TFIIH. While the ring-shaped core domain was surprisingly congruent between the two structures, trypanosome TFIIH lacked the knob-like CAK moiety and exhibited extra densities on either side of the ring, presumably due to the specific subunits.

MicroRNA-338-3p and microRNA-451 contribute to the formation of basolateral polarity in epithelial cells

2009-06-21

MicroRNAs are small noncoding RNA species, some of which are playing important roles in cell differentiation. However, the level of participations of microRNAs in epithelial cell differentiation is largely unknown. Here, utilizing an epithelial differentiation model with T84 cells, we demonstrate that miR-338-3p and miR-451 contribute to the formation of epithelial basolateral polarity by facilitating translocalization of β1 integrin to the basolateral membrane. Among 250 microRNAs screened in this study, the expression levels of four microRNAs (miR-33a, 210, 338-3p and 451) were significantly elevated in the differentiated stage of T84 cells, when epithelial cell polarity was established. To investigate the involvement of these microRNAs in terms of epithelial cell polarity, we executed loss-of- and gain-of-function analyses of these microRNAs. The blockade of endogenous miR-338-3p or miR-451 via each microRNA-specific antisense oligonucleotides inhibited the translocalization of β1 integrin to the basolateral membrane, whereas inhibition of miR-210 or miR-33a had no effect on it. On the other hand, simultaneous transfection of synthetic miR-338-3p and miR-451 accelerated the translocalization of β1 integrin to the basolateral membrane, although the introduction of individual synthetic microRNAs exhibited no effect. Therefore, we concluded that both miR-338-3p and miR-451 are necessary for the development of epithelial cell polarity.