Search

Summary Mating-type switching in yeast occurs through gene conversion between the MAT locus and one of two silent loci ( HML or HMR ) on opposite ends of the chromosome. MAT a cells choose HML as template, whereas MAT α cells use HMR . The recombination enhancer (RE) located on the left arm regulates this process. One long-standing hypothesis is that switching is guided by mating-type-specific and possibly RE-dependent chromosome folding. Here, we use Hi-C, 5C, and live-cell imaging to characterize the conformation of chromosome III in both mating types. We discovered a mating-type-specific conformational difference in the left arm. Deletion of a 1-kb subregion within the RE, which is not necessary during switching, abolished mating-type-dependent chromosome folding. The RE is therefore a composite element with one subregion essential for donor selection during switching and a separate region involved in modulating chromosome conformation. [ABSTRACT FROM AUTHOR]

How to describe the multiple chromosome structures that underlie interactions among genome loci and how to quantify the occurrence of these structures in a cell population remain important challenges to solve, which can be addressed via a proper demultiplexing of chromosome capture conformation related data. Here, we first aim to review two main methodologies that have been proposed to tackle this problem: restrained-based methods, in which the resulting chromosome structures stem from the multiple solutions of a distance satisfaction problem; and thermodynamic-based methods, in which the structures stem from the simulation of polymer models. Next, we propose a novel demultiplexing method based on a matrix decomposition of contact maps. To this end, we extend the notion of topologically associated domains (TADs) by introducing that of statistical interaction domains (SIDs). SIDs can overlap and occur in a cell population at certain frequencies, and we propose a simple method to estimate these frequency values. As an application, we show that SIDs that measure 100 kb to tens of Mb long occur both frequently and specifically in the human genome. [ABSTRACT FROM AUTHOR]

Chromosomes are large polymer molecules composed of nucleotides. In some species, such as humans, this polymer can sum up to meters long and still be properly folded within the nuclear space of few microns in size. The exact mechanisms of how the meters long DNA is folded into the nucleus, as well as how the regulatory machinery can access it, is to a large extend still a mystery. However, and thanks to newly developed molecular, genomic and computational approaches based on the Chromosome Conformation Capture (3C) technology, we are now obtaining insight on how genomes are spatially organized. Here we review a new family of computational approaches that aim at using 3C-based data to obtain spatial restraints for modeling genomes and genomic domains. [ABSTRACT FROM AUTHOR]

RNA is not regarded anymore as a simple transfer molecule between DNA and proteins. Indeed, over the past decades a plethora of new functional roles have been assigned to RNA molecules. Such functions are carried out either by RNA molecules alone or through interactions with DNA, other RNA molecules, or proteins. In all cases, the structure that the RNA molecule adopts will impact its function, as it happens with proteins. Therefore, to fully characterize the function of an RNA molecule, its structure needs to be either determined by experiments or predicted by computation. Unfortunately, our knowledge of the atomic mechanism by which RNA molecules adopt their biological active structures is still limited. Such hurdle is now being addressed by the development of new computational methods for RNA structure prediction, which complement experimental methods such as X-ray crystallography, nuclear magnetic resonance, small-angle X-ray scattering, and cryo-electron microscopy. This software focus is not dedicated to a single computational method but aims at outlining the most adopted methods for computational RNA structure prediction. WIREs Comput Mol Sci 2015, 5:56-61. doi: 10.1002/wcms.1198 For further resources related to this article, please visit the . Conflict of interest: The authors have declared no conflicts of interest for this article. [ABSTRACT FROM AUTHOR]

Mycobacterium tuberculosis, the causative agent of tuberculosis (TB), infects an estimated two billion people worldwide and is the leading cause of mortality due to infectious disease. The development of new anti-TB therapeutics is required, because of the emergence of multi-drug resistance strains as well as co-infection with other pathogens, especially HIV. Recently, the pharmaceutical company GlaxoSmithKline published the results of a high-throughput screen (HTS) of their two million compound library for anti-mycobacterial phenotypes. The screen revealed 776 compounds with significant activity against the M. tuberculosis H37Rv strain, including a subset of 177 prioritized compounds with high potency and low in vitro cytotoxicity. The next major challenge is the identification of the target proteins. Here, we use a computational approach that integrates historical bioassay data, chemical properties and structural comparisons of selected compounds to propose their potential targets in M. tuberculosis. We predicted 139 target - compound links, providing a necessary basis for further studies to characterize the mode of action of these compounds. The results from our analysis, including the predicted structural models, are available to the wider scientific community in the open source mode, to encourage further development of novel TB therapeutics. [ABSTRACT FROM AUTHOR]

How DNA is organized in three dimensions inside the cell nucleus and how this affects the ways in which cells access, read and interpret genetic information are among the longest standing questions in cell biology. Using newly developed molecular, genomic and computational approaches based on the chromosome conformation capture technology (such as 3C, 4C, 5C and Hi-C), the spatial organization of genomes is being explored at unprecedented resolution. Interpreting the increasingly large chromatin interaction data sets is now posing novel challenges. Here we describe several types of statistical and computational approaches that have recently been developed to analyse chromatin interaction data. [ABSTRACT FROM AUTHOR]

Recent advances in determination of the high-resolution structure of membrane proteins now enable analysis of the main features of amino acids in transmembrane (TM) segments in comparison with amino acids in water-soluble helices. In this work, we conducted a large-scale analysis of the prevalent locations of amino acids by using a data set of 170 structures of integral membrane proteins obtained from the MPtopo database and 930 structures of water-soluble helical proteins obtained from the protein data bank. Large hydrophobic amino acids (Leu, Val, Ile, and Phe) plus Gly were clearly prevalent in TM helices whereas polar amino acids (Glu, Lys, Asp, Arg, and Gln) were less frequent in this type of helix. The distribution of amino acids along TM helices was also examined. As expected, hydrophobic and slightly polar amino acids are commonly found in the hydrophobic core of the membrane whereas aromatic (Trp and Tyr), Pro, and the hydrophilic amino acids (Asn, His, and Gln) occur more frequently in the interface regions. Charged amino acids are also statistically prevalent outside the hydrophobic core of the membrane, and whereas acidic amino acids are frequently found at both cytoplasmic and extra-cytoplasmic interfaces, basic amino acids cluster at the cytoplasmic interface. These results strongly support the experimentally demonstrated biased distribution of positively charged amino acids (that is, the so-called the positive-inside rule) with structural data. [ABSTRACT FROM AUTHOR]

Abstract: The three-dimensional (3D) architecture of a genome determines the spatial localization of regulatory elements and the genes they regulate. Thus, elucidating the 3D structure of a genome may result in significant insights about how genes are regulated. The current state-of-the art in experimental methods, including light microscopy and cell/molecular biology, are now able to provide detailed information on the position of genes and their interacting partners. However, such methods by themselves are not able to determine the high-resolution 3D structure of genomes or genomic domains. Here we describe a computational module of the Integrative Modeling Platform (IMP, http://www.integrativemodeling.org) that uses chromosome conformation capture data to determine the 3D architecture of genomic domains and entire genomes at unprecedented resolutions. This approach, through the visualization of looping interactions between distal regulatory elements, allows characterizing global chromatin features and their relation to gene expression. We illustrate our work by outlining the determination of the 3D architecture of the α-globin domain in the human genome. [Copyright &y& Elsevier]

The three-dimensional (3D) architecture of a genome is non-random and known to facilitate the spatial colocalization of regulatory elements with the genes they regulate. Determining the 3D structure of a genome may therefore probe an essential step in characterizing how genes are regulated. Currently, there are several experimental and theoretical approaches that aim at determining the 3D structure of genomes and genomic domains; however, approaches integrating experiments and computation to identify the most likely 3D folding of a genome at medium to high resolutions have not been widely explored. Here, we review existing methodologies and propose that the integrative modeling platform (), a computational package developed for structurally characterizing protein assemblies, could be used for integrating diverse experimental data towards the determination of the 3D architecture of genomic domains and entire genomes at unprecedented resolution. Our approach, through the visualization of looping interactions between distal regulatory elements, will allow for the characterization of global chromatin features and their relation to gene expression. We illustrate our work by outlining the recent determination of the 3D architecture of the α-globin domain in the human genome. [ABSTRACT FROM AUTHOR]

We developed a general approach that combines chromosome conformation capture carbon copy (5C) with the Integrated Modeling Platform (IMP) to generate high-resolution three-dimensional models of chromatin at the megabase scale. We applied this approach to the ENm008 domain on human chromosome 16, containing the α-globin locus, which is expressed in K562 cells and silenced in lymphoblastoid cells (GM12878). The models accurately reproduce the known looping interactions between the α-globin genes and their distal regulatory elements. Further, we find using our approach that the domain folds into a single globular conformation in GM12878 cells, whereas two globules are formed in K562 cells. The central cores of these globules are enriched for transcribed genes, whereas nontranscribed chromatin is more peripheral. We propose that globule formation represents a higher-order folding state related to clustering of transcribed genes around shared transcription machineries, as previously observed by microscopy. [ABSTRACT FROM AUTHOR]

Background: In recent years, the number of available RNA structures has rapidly grown reflecting the increased interest on RNA biology. Similarly to the studies carried out two decades ago for proteins, which gave the fundamental grounds for developing comparative protein structure prediction methods, we are now able to quantify the relationship between sequence and structure conservation in RNA. Results: Here we introduce an all-against-all sequence- and three-dimensional (3D) structure-based comparison of a representative set of RNA structures, which have allowed us to quantitatively confirm that: (i) there is a measurable relationship between sequence and structure conservation that weakens for alignments resulting in below 60% sequence identity, (ii) evolution tends to conserve more RNA structure than sequence, and (iii) there is a twilight zone for RNA homology detection. Discussion: The computational analysis here presented quantitatively describes the relationship between sequence and structure for RNA molecules and defines a twilight zone region for detecting RNA homology. Our work could represent the theoretical basis and limitations for future developments in comparative RNA 3D structure prediction. [ABSTRACT FROM AUTHOR]

Background: Conventional patent-based drug development incentives work badly for the developing world, where commercial markets are usually small to non-existent. For this reason, the past decade has seen extensive experimentation with alternative R&D institutions ranging from private–public partnerships to development prizes. Despite extensive discussion, however, one of the most promising avenues—open source drug discovery—has remained elusive. We argue that the stumbling block has been the absence of a critical mass of preexisting work that volunteers can improve through a series of granular contributions. Historically, open source software collaborations have almost never succeeded without such "kernels". Methodology/Principal Findings: Here, we use a computational pipeline for: (i) comparative structure modeling of target proteins, (ii) predicting the localization of ligand binding sites on their surfaces, and (iii) assessing the similarity of the predicted ligands to known drugs. Our kernel currently contains 143 and 297 protein targets from ten pathogen genomes that are predicted to bind a known drug or a molecule similar to a known drug, respectively. The kernel provides a source of potential drug targets and drug candidates around which an online open source community can nucleate. Using NMR spectroscopy, we have experimentally tested our predictions for two of these targets, confirming one and invalidating the other. Conclusions/Significance: The TDI kernel, which is being offered under the Creative Commons attribution share-alike license for free and unrestricted use, can be accessed on the World Wide Web at http://www.tropicaldisease.org. We hope that the kernel will facilitate collaborative efforts towards the discovery of new drugs against parasites that cause tropical diseases. Author Summary: Open source drug discovery, a promising alternative avenue to conventional patent-based drug development, has so far remained elusive with few exceptions. A major stumbling block has been the absence of a critical mass of preexisting work that volunteers can improve through a series of granular contributions. This paper introduces the results from a newly assembled computational pipeline for identifying protein targets for drug discovery in ten organisms that cause tropical diseases. We have also experimentally tested two promising targets for their binding to commercially available drugs, validating one and invalidating the other. The resulting kernel provides a base of drug targets and lead candidates around which an open source community can nucleate. We invite readers to donate their judgment and in silico and in vitro experiments to develop these targets to the point where drug optimization can begin. [ABSTRACT FROM AUTHOR]

The accuracy of an alignment between two protein sequences can be improved by including other detectably related sequences in the comparison. We optimize and benchmark such an approach that relies on aligning two multiple sequence alignments, each one including one of the two protein sequences. Thirteen different protocols for creating and comparing profiles corresponding to the multiple sequence alignments are implemented in the SALIGN command of MODELLER. A test set of 200 pairwise, structure-based alignments with sequence identities below 40% is used to benchmark the 13 protocols as well as a number of previously described sequence alignment methods, including heuristic pairwise sequence alignment by BLAST, pairwise sequence alignment by global dynamic programming with an affine gap penalty function by the ALIGN command of MODELLER, sequence-profile alignment by PSI-BLAST, Hidden Markov Model methods implemented in SAM and LOBSTER, pairwise sequence alignment relying on predicted local structure by SEA, and multiple sequence alignment by CLUSTALW and COMPASS. The alignment accuracies of the best new protocols were significantly better than those of the other tested methods. For example, the fraction of the correctly aligned residues relative to the structure-based alignment by the best protocol is 56%, which can be compared with the accuracies of 26%, 42%, 43%, 48%, 50%, 49%, 43%, and 43% for the other methods, respectively. The new method is currently applied to large-scale comparative protein structure modeling of all known sequences. [ABSTRACT FROM AUTHOR]