Correlation between nucleotide composition and folding energy of coding sequences with special attention to wobble bases

Background: The secondary structure and complexity of mRNA influences its accessibility to regulatory molecules (proteins, micro-RNAs), its stability and its level of expression. The mobile elements of the RNA sequence, the wobble bases, are expected to regulate the formation of structures encompassing coding sequences. Results: The sequence/folding energy (FE) relationship was studied by statistical, bioinformatic methods in 90 CDS containing 26,370 codons. I found that the FE (dG) associated with coding sequences is significant and negative (407 kcal/1000 bases, mean +/- S.E.M.) indicating that these sequences are able to form structures. However, the FE has only a small free component, less than 10% of the total. The contribution of the 1st and 3rd codon bases to the FE is larger than the contribution of the 2nd (central) bases. It is possible to achieve a ~ 4-fold change in FE by altering the wobble bases in synonymous codons. The sequence/FE relationship can be described with a simple algorithm, and the total FE can be predicted solely from the sequence composition of the nucleic acid. The contributions of different synonymous codons to the FE are additive and one codon cannot replace another. The accumulated contributions of synonymous codons of an amino acid to the total folding energy of an mRNA is strongly correlated to the relative amount of that amino acid in the translated protein. Conclusion: Synonymous codons are not interchangable with regard to their role in determining the mRNA FE and the relative amounts of amino acids in the translated protein, even if they are indistinguishable in respect of amino acid coding.Comment: 14 pages including 6 figures and 1 tabl

The Set2/Rpd3S Pathway Suppresses Cryptic Transcription without Regard to Gene Length or Transcription Frequency

In cells lacking the histone methyltransferase Set2, initiation of RNA polymerase II transcription occurs inappropriately within the protein-coding regions of genes, rather than being restricted to the proximal promoter. It was previously reported that this “cryptic” transcription occurs preferentially in long genes, and in genes that are infrequently transcribed. Here, we mapped the transcripts produced in an S. cerevisiae strain lacking Set2, and applied rigorous statistical methods to identify sites of cryptic transcription at high resolution. We find that suppression of cryptic transcription occurs independent of gene length or transcriptional frequency. Our conclusions differ with those reported previously because we obtained a higher-resolution dataset, we accounted for the fact that gene length and transcriptional frequency are not independent variables, and we accounted for several ascertainment biases that make cryptic transcription easier to detect in long, infrequently transcribed genes. These new results and conclusions have implications for many commonly used genomic analysis approaches, and for the evolution of high-fidelity RNA polymerase II transcriptional initiation in eukaryotes

Intragenic DNA methylation: implications of this epigenetic mechanism for cancer research

Epigenetics is the study of all mechanisms that regulate gene transcription and genome stability that are maintained throughout the cell division, but do not include the DNA sequence itself. The best-studied epigenetic mechanism to date is DNA methylation, where methyl groups are added to the cytosine base within cytosine–guanine dinucleotides (CpG sites). CpGs are frequently clustered in high density (CpG islands (CGIs)) at the promoter of over half of all genes. Current knowledge of transcriptional regulation by DNA methylation centres on its role at the promoter where unmethylated CGIs are present at most actively transcribed genes, whereas hypermethylation of the promoter results in gene repression. Over the last 5 years, research has gradually incorporated a broader understanding that methylation patterns across the gene (so-called intragenic or gene body methylation) may have a role in transcriptional regulation and efficiency. Numerous genome-wide DNA methylation profiling studies now support this notion, although whether DNA methylation patterns are a cause or consequence of other regulatory mechanisms is not yet clear. This review will examine the evidence for the function of intragenic methylation in gene transcription, and discuss the significance of this in carcinogenesis and for the future use of therapies targeted against DNA methylation

Topology by Design in Magnetic nano-Materials: Artificial Spin Ice

Artificial Spin Ices are two dimensional arrays of magnetic, interacting nano-structures whose geometry can be chosen at will, and whose elementary degrees of freedom can be characterized directly. They were introduced at first to study frustration in a controllable setting, to mimic the behavior of spin ice rare earth pyrochlores, but at more useful temperature and field ranges and with direct characterization, and to provide practical implementation to celebrated, exactly solvable models of statistical mechanics previously devised to gain an understanding of degenerate ensembles with residual entropy. With the evolution of nano--fabrication and of experimental protocols it is now possible to characterize the material in real-time, real-space, and to realize virtually any geometry, for direct control over the collective dynamics. This has recently opened a path toward the deliberate design of novel, exotic states, not found in natural materials, and often characterized by topological properties. Without any pretense of exhaustiveness, we will provide an introduction to the material, the early works, and then, by reporting on more recent results, we will proceed to describe the new direction, which includes the design of desired topological states and their implications to kinetics.Comment: 29 pages, 13 figures, 116 references, Book Chapte

A mathematical and computational review of Hartree-Fock SCF methods in Quantum Chemistry

We present here a review of the fundamental topics of Hartree-Fock theory in Quantum Chemistry. From the molecular Hamiltonian, using and discussing the Born-Oppenheimer approximation, we arrive to the Hartree and Hartree-Fock equations for the electronic problem. Special emphasis is placed in the most relevant mathematical aspects of the theoretical derivation of the final equations, as well as in the results regarding the existence and uniqueness of their solutions. All Hartree-Fock versions with different spin restrictions are systematically extracted from the general case, thus providing a unifying framework. Then, the discretization of the one-electron orbitals space is reviewed and the Roothaan-Hall formalism introduced. This leads to a exposition of the basic underlying concepts related to the construction and selection of Gaussian basis sets, focusing in algorithmic efficiency issues. Finally, we close the review with a section in which the most relevant modern developments (specially those related to the design of linear-scaling methods) are commented and linked to the issues discussed. The whole work is intentionally introductory and rather self-contained, so that it may be useful for non experts that aim to use quantum chemical methods in interdisciplinary applications. Moreover, much material that is found scattered in the literature has been put together here to facilitate comprehension and to serve as a handy reference.Comment: 64 pages, 3 figures, tMPH2e.cls style file, doublesp, mathbbol and subeqn package

The role of multiple marks in epigenetic silencing and the emergence of a stable bivalent chromatin state

We introduce and analyze a minimal model of epigenetic silencing in budding yeast, built upon known biomolecular interactions in the system. Doing so, we identify the epigenetic marks essential for the bistability of epigenetic states. The model explicitly incorporates two key chromatin marks, namely H4K16 acetylation and H3K79 methylation, and explores whether the presence of multiple marks lead to a qualitatively different systems behavior. We find that having both modifications is important for the robustness of epigenetic silencing. Besides the silenced and transcriptionally active fate of chromatin, our model leads to a novel state with bivalent (i.e., both active and silencing) marks under certain perturbations (knock-out mutations, inhibition or enhancement of enzymatic activity). The bivalent state appears under several perturbations and is shown to result in patchy silencing. We also show that the titration effect, owing to a limited supply of silencing proteins, can result in counter-intuitive responses. The design principles of the silencing system is systematically investigated and disparate experimental observations are assessed within a single theoretical framework. Specifically, we discuss the behavior of Sir protein recruitment, spreading and stability of silenced regions in commonly-studied mutants (e.g., sas2, dot1) illuminating the controversial role of Dot1 in the systems biology of yeast silencing.Comment: Supplementary Material, 14 page

From DNA sequence to application: possibilities and complications

The development of sophisticated genetic tools during the past 15 years have facilitated a tremendous increase of fundamental and application-oriented knowledge of lactic acid bacteria (LAB) and their bacteriophages. This knowledge relates both to the assignments of open reading frames (ORF’s) and the function of non-coding DNA sequences. Comparison of the complete nucleotide sequences of several LAB bacteriophages has revealed that their chromosomes have a fixed, modular structure, each module having a set of genes involved in a specific phase of the bacteriophage life cycle. LAB bacteriophage genes and DNA sequences have been used for the construction of temperature-inducible gene expression systems, gene-integration systems, and bacteriophage defence systems. The function of several LAB open reading frames and transcriptional units have been identified and characterized in detail. Many of these could find practical applications, such as induced lysis of LAB to enhance cheese ripening and re-routing of carbon fluxes for the production of a specific amino acid enantiomer. More knowledge has also become available concerning the function and structure of non-coding DNA positioned at or in the vicinity of promoters. In several cases the mRNA produced from this DNA contains a transcriptional terminator-antiterminator pair, in which the antiterminator can be stabilized either by uncharged tRNA or by interaction with a regulatory protein, thus preventing formation of the terminator so that mRNA elongation can proceed. Evidence has accumulated showing that also in LAB carbon catabolite repression in LAB is mediated by specific DNA elements in the vicinity of promoters governing the transcription of catabolic operons. Although some biological barriers have yet to be solved, the vast body of scientific information presently available allows the construction of tailor-made genetically modified LAB. Today, it appears that societal constraints rather than biological hurdles impede the use of genetically modified LAB.

Silent but Not Static: Accelerated Base-Pair Substitution in Silenced Chromatin of Budding Yeasts

Subtelomeric DNA in budding yeasts, like metazoan heterochromatin, is gene poor, repetitive, transiently silenced, and highly dynamic. The rapid evolution of subtelomeric regions is commonly thought to arise from transposon activity and increased recombination between repetitive elements. However, we found evidence of an additional factor in this diversification. We observed a surprising level of nucleotide divergence in transcriptionally silenced regions in inter-species comparisons of Saccharomyces yeasts. Likewise, intra-species analysis of polymorphisms also revealed increased SNP frequencies in both intergenic and synonymous coding positions of silenced DNA. This analysis suggested that silenced DNA in Saccharomyces cerevisiae and closely related species had increased single base-pair substitution that was likely due to the effects of the silencing machinery on DNA replication or repair

Two-Particle-Self-Consistent Approach for the Hubbard Model

Even at weak to intermediate coupling, the Hubbard model poses a formidable challenge. In two dimensions in particular, standard methods such as the Random Phase Approximation are no longer valid since they predict a finite temperature antiferromagnetic phase transition prohibited by the Mermin-Wagner theorem. The Two-Particle-Self-Consistent (TPSC) approach satisfies that theorem as well as particle conservation, the Pauli principle, the local moment and local charge sum rules. The self-energy formula does not assume a Migdal theorem. There is consistency between one- and two-particle quantities. Internal accuracy checks allow one to test the limits of validity of TPSC. Here I present a pedagogical review of TPSC along with a short summary of existing results and two case studies: a) the opening of a pseudogap in two dimensions when the correlation length is larger than the thermal de Broglie wavelength, and b) the conditions for the appearance of d-wave superconductivity in the two-dimensional Hubbard model.Comment: Chapter in "Theoretical methods for Strongly Correlated Systems", Edited by A. Avella and F. Mancini, Springer Verlag, (2011) 55 pages. Misprint in Eq.(23) corrected (thanks D. Bergeron

Improved ChIP-chip analysis by a mixture model approach

<p>Abstract</p> <p>Background</p> <p>Microarray analysis of immunoprecipitated chromatin (ChIP-chip) has evolved from a novel technique to a standard approach for the systematic study of protein-DNA interactions. In ChIP-chip, sites of protein-DNA interactions are identified by signals from the hybridization of selected DNA to tiled oligomers and are graphically represented as peaks. Most existing methods were designed for the identification of relatively sparse peaks, in the presence of replicates.</p> <p>Results</p> <p>We propose a data normalization method and a statistical method for peak identification from ChIP-chip data based on a mixture model approach. In contrast to many existing methods, including methods that also employ mixture model approaches, our method is more flexible by imposing less restrictive assumptions and allowing a relatively large proportion of peak regions. In addition, our method does not require experimental replicates and is computationally efficient. We compared the performance of our method with several representative existing methods on three datasets, including a spike-in dataset. These comparisons demonstrate that our approach is more robust and has comparable or higher power than the other methods, especially in the context of abundant peak regions.</p> <p>Conclusion</p> <p>Our data normalization and peak detection methods have improved performance to detect peak regions in ChIP-chip data.</p