Search

In the protein sequence space, natural proteins form clusters of families which are characterized by their unique native folds whereas the great majority of random polypeptides are neither clustered nor foldable to unique structures. Since a given polypeptide can be either foldable or unfoldable, a kind of "folding transition" is expected at the boundary of a protein family in the sequence space. By Monte Carlo simulations of a statistical mechanical model of protein sequence alignment that coherently incorporates both short-range and long-range interactions as well as variable-length insertions to reproduce the statistics of the multiple sequence alignment of a given protein family, we demonstrate the existence of such transition between natural-like sequences and random sequences in the sequence subspaces for 15 domain families of various folds. The transition was found to be highly cooperative and two-state-like. Furthermore, enforcing or suppressing consensus residues on a few of the well-conserved sites enhanced or diminished, respectively, the natural-like pattern formation over the entire sequence. In most families, the key sites included ligand binding sites. These results suggest some selective pressure on the key residues, such as ligand binding activity, may cooperatively facilitate the emergence of a protein family during evolution. From a more practical aspect, the present results highlight an essential role of long-range effects in precisely defining protein families, which are absent in conventional sequence models., 13 pages, 7 figures, 2 tables (a new subsection added)

We re-examine a site-binding approach independently proposed by Schellman (Schellman, J.A. (1958) Compt. rend. Lab. Carlsberg Ser. Chim. 30, 439-449) and Aune and Tanford (Aune, K.C. and Tanford, D. (1969) Biochemistry, 8, 4586-4590) for explicitly including the denaturant concentration within the protein unfolding equilibrium. We extend and formalize the approach through development of a multi-dimensional analytical model in which the folding reaction coordinate is defined by the number of denaturant molecules bound to sites located on either the initially folded, or unfolded, states of the protein. We use the developed method to re-examine the mechanistic determinants underlying the sigmoidal shape of the unfolding transition. A natural feature of our method is that it presents a landscape picture of the denaturant induced protein unfolding reaction.

The Protein Data Bank Japan (PDBj), a member of the worldwide Protein Data Bank (wwPDB), accepts and processes the deposited data of experimentally determined biological macromolecular structures. In addition to archiving the PDB data in collaboration with the other wwPDB partners, PDBj also provides a wide range of original and unique services and tools, which are continuously improved and updated. Here, we report the new RDB PDBj Mine 2, the WebGL molecular viewer Molmil, the ProMode-Elastic server for normal mode analysis, a virtual reality system for the eF-site protein electrostatic molecular surfaces, the extensions of the Omokage search for molecular shape similarity, and the integration of PDBj and BMRB searches.

A grand canonical Monte Carlo (MC) algorithm is presented for studying the lattice gas model (LGM) of multiple protein sequence alignment, which coherently combines long-range interactions and variable-length insertions. MC simulations are used for both parameter optimization of the model and production runs to explore the sequence subspace around a given protein family. In this Note, I describe the details of the MC algorithm as well as some preliminary results of MC simulations with various temperatures and chemical potentials, and compare them with the mean-field approximation. The existence of a two-state transition in the sequence space is suggested for the SH3 domain family, and inappropriateness of the mean-field approximation for the LGM is demonstrated., Comment: 24 pages, 7 figures; corrected errors

We report on the activities of the 2015 edition of the BioHackathon, an annual event that brings together researchers and developers from around the world to develop tools and technologies that promote the reusability of biological data. We discuss issues surrounding the representation, publication, integration, mining and reuse of biological data and metadata across a wide range of biomedical data types of relevance for the life sciences, including chemistry, genotypes and phenotypes, orthology and phylogeny, proteomics, genomics, glycomics, and metabolomics. We describe our progress to address ongoing challenges to the reusability and reproducibility of research results, and identify outstanding issues that continue to impede the progress of bioinformatics research. We share our perspective on the state of the art, continued challenges, and goals for future research and development for the life sciences Semantic Web.

The multiple sequence alignment (MSA) of a protein family provides a wealth of information in terms of the conservation pattern of amino acid residues not only at each alignment site but also between distant sites. In order to statistically model the MSA incorporating both short-range and long-range correlations as well as insertions, I have derived a lattice gas model of the MSA based on the principle of maximum entropy. The partition function, obtained by the transfer matrix method with a mean-field approximation, accounts for all possible alignments with all possible sequences. The model parameters for short-range and long-range interactions were determined by a self-consistent condition and by a Gaussian approximation, respectively. Using this model with and without long-range interactions, I analyzed the globin and V-set domains by increasing the "temperature" and by "mutating" a site. The correlations between residue conservation and various measures of the system's stability indicate that the long-range interactions make the conservation pattern more specific to the structure, and increasingly stabilize better conserved residues., 21 pages (2-column), 13 figures

Conrad H. Waddington discovered the phenomenon of genetic assimilation through a series of experiments on fruit flies. In those experiments, artificially exerted environmental stress induced plastic phe-notypic changes in the fruit flies, but after some generations the same phenotypic variant started to appear without the environmental stress. Both the initial state (where the phenotypic changes were environmentally induced and plastic) and the final state (where the phenotypic changes were genetically fixed and constitutive) are experimental facts. However, it remains unclear how the environmentally induced phenotypic change in the first generation becomes genetically fixed in the central process of genetic assimilation itself. We have argued that the key to understanding the mechanism of genetic assimilation lies in epigenetics, and proposed the “cooperative model” in which the evolutionary process depends on both genetic and epigenetic factors. Evolutionary simulations based on the cooperative model reproduced the process of genetic assimilation. Detailed analysis of the trajectories has revealed genetic assimilation is a process in which epigenetically induced phe-notypic changes are incrementally and statistically replaced with multiple minor genetic mutations through natural selection. In this scenario, epigenetic and genetic changes may be considered as mutually independent but equivalent in terms of their effects on phenotypic changes. This finding rejects the common (and confused) hypothesis that epigenetically induced phe-notypic changes depend on genetic mutations. Furthermore, we argue that transgenerational epigenetic inheritance is not required for evolution by genetic assimilation.

Motivation Biological data and knowledge bases increasingly rely on Semantic Web technologies and the use of knowledge graphs for data integration, retrieval and federated queries. In the past years, feature learning methods that are applicable to graph-structured data are becoming available, but have not yet widely been applied and evaluated on structured biological knowledge. Results: We develop a novel method for feature learning on biological knowledge graphs. Our method combines symbolic methods, in particular knowledge representation using symbolic logic and automated reasoning, with neural networks to generate embeddings of nodes that encode for related information within knowledge graphs. Through the use of symbolic logic, these embeddings contain both explicit and implicit information. We apply these embeddings to the prediction of edges in the knowledge graph representing problems of function prediction, finding candidate genes of diseases, protein-protein interactions, or drug target relations, and demonstrate performance that matches and sometimes outperforms traditional approaches based on manually crafted features. Our method can be applied to any biological knowledge graph, and will thereby open up the increasing amount of Semantic Web based knowledge bases in biology to use in machine learning and data analytics. Availability and implementation https://github.com/bio-ontology-research-group/walking-rdf-and-owl Supplementary information Supplementary data are available at Bioinformatics online.

The Protein Data Bank Japan (PDBj), a member of the worldwide Protein Data Bank (wwPDB), accepts and processes the deposited data of experimentally determined biological macromolecular structures. In addition to archiving the PDB data in collaboration with the other wwPDB partners, PDBj also provides a wide range of original and unique services and tools, which are continuously improved and updated. Here, we report the new RDB PDBj Mine 2, the WebGL molecular viewer Molmil, the ProMode‐Elastic server for normal mode analysis, a virtual reality system for the eF‐site protein electrostatic molecular surfaces, the extensions of the Omokage search for molecular shape similarity, and the integration of PDBj and BMRB searches.

We propose the cooperative model of phenotype-driven evolution, in which natural selection operates on a phenotype caused by both genetic and epigenetic factors. The conventional theory of evolutionary synthesis assumes that a phenotypic value (P) is the sum of genotypic value (G) and environmental deviation (E), P=G+E, where E is the fluctuations of the phenotype among individuals in the absence of environmental changes. In contrast, the cooperative model assumes that an evolution is triggered by an environmental change and individuals respond to the change by phenotypic plasticity (epigenetic changes). The phenotypic plasticity, while essentially qualitative, is denoted by a quantitative value F which is modeled as a normal random variable like E, but with a much larger variance. Thus, the fundamental equation of the cooperative model is given as P=G+F where F includes the effect of E. Computer simulations using a genetic algorithm demonstrated that the cooperative model realized much faster evolution than the evolutionary synthesis. This accelerated evolution was found to be due to the cumulative evolution made possible by a ratchet mechanism due to the epigenetic contribution to the phenotypic value. The cooperative model can well account for the phenomenon of genetic assimilation, which, in turn, suggests the mechanism of cumulative selection. The cooperative model may also serve as a theoretical basis to understand various ideas and phenomena of the phenotype-driven evolution such as genetic assimilation, the theory of facilitated phenotypic variation, and epigenetic inheritance over generations.

The epigenome, i.e. the whole of chromatin modifications, is transferred from mother to daughter cells during cell differentiation. When de novo chromatin modifications (establishment or erasure of, respectively, new or pre-existing DNA methylations and/or histone modifications) are made in a daughter cell, however, it has a different epigenome than its mother cell. Although de novo chromatin modifications are an important event that comprises elementary processes of cell differentiation, its molecular mechanism remains poorly understood. We argue in this Letter that a key to solving this problem lies in understanding the role of long non-coding RNAs (lncRNAs)- a type of RNA that is becoming increasingly prominent in epigenetic studies. Many studies show that lncRNAs form ribonucleo-protein complexes in the nucleus and are involved in chromatin modifications. However, chromatin-modifying enzymes lack the information about genomic positions on which they act. It is known, on the other hand, that a single-stranded RNA in general can bind to a double-stranded DNA to form a triple helix. If each lncRNA forms a ribonucleo-protein complex with chromatin-modifying enzymes on one hand and, at the same time, a triple helix with a genomic region based on its specific nucleotide sequence on the other hand, it can induce de novo chromatin modifications at specific sites. Thus, the great variety of lncRNAs can be explained by the requirement for the diversity of "genomic address codes" specific to their cognate genomic regions where de novo chromatin modifications take place., Comment: 12 pages, 1 figure

Life science research now heavily relies on all sorts of databases for genome sequences, transcription, protein three-dimensional (3D) structures, protein-protein interactions, phenotypes and so forth. The knowledge accumulated by all the omics research is so vast that a computer-aided search of data is now a prerequisite for starting a new study. In addition, a combinatory search throughout these databases has a chance to extract new ideas and new hypotheses that can be examined by wet-lab experiments. By virtually integrating the related databases on the Internet, we have built a new web application that facilitates life science researchers for retrieving experts' knowledge stored in the databases and for building a new hypothesis of the research target. This web application, named VaProS, puts stress on the interconnection between the functional information of genome sequences and protein 3D structures, such as structural effect of the gene mutation. In this manuscript, we present the notion of VaProS, the databases and tools that can be accessed without any knowledge of database locations and data formats, and the power of search exemplified in quest of the molecular mechanisms of lysosomal storage disease. VaProS can be freely accessed at http://p4d-info.nig.ac.jp/vapros/ .

Date hub proteins are a type of proteins that show multispecificity in a time-dependent manner. To understand dynamic aspects of such multispecificity we studied Ubiquitin as a typical example of a date hub protein. Here we analyzed 9 biologically relevant Ubiquitin-protein (ligand) heterodimer structures by using normal mode analysis based on an elastic network model. Our result showed that the self-coupled motion of Ubiquitin in the complex, rather than its ligand-coupled motion, is similar to the motion of Ubiquitin in the unbound condition. The ligand-coupled motions are correlated to the conformational change between the unbound and bound conditions of Ubiquitin. Moreover, ligand-coupled motions favor the formation of the bound states, due to its in-phase movements of the contacting atoms at the interface. The self-coupled motions at the interface indicated loss of conformational entropy due to binding. Therefore, such motions disfavor the formation of the bound state. We observed that the ligand-coupled motions are embedded in the motions of unbound Ubiquitin. In conclusion, multispecificity of Ubiquitin can be characterized by an intricate balance of the ligand- and self-coupled motions, both of which are embedded in the motions of the unbound form.

The Protein Data Bank Japan (PDBj, http://pdbj.org) is a member of the worldwide Protein Data Bank (wwPDB) and accepts and processes the deposited data of experimentally determined macromolecular structures. While maintaining the archive in collaboration with other wwPDB partners, PDBj also provides a wide range of services and tools for analyzing structures and functions of proteins, which are summarized in this article. To enhance the interoperability of the PDB data, we have recently developed PDB/RDF, PDB data in the Resource Description Framework (RDF) format, along with its ontology in the Web Ontology Language (OWL) based on the PDB mmCIF Exchange Dictionary. Being in the standard format for the Semantic Web, the PDB/RDF data provide a means to integrate the PDB with other biological information resources.

Cellular functions of an organism are maintained by protein–protein interactions. Those proteins that bind multiple partners asynchronously (date hub proteins) are important to make the interaction network coordinated. It is known that many date hub proteins bind different partners at overlapping (OV) interfaces. To understand how OV interfaces of date hub proteins can recognize multiple partners, we analyzed the difference between OV and non-overlapping (Non-OV) regions of interfaces involved in the binding of different partners. By using the structures of 16 date hub proteins with various interaction partners (ranging from 5 to 33), we compared buried surface area, compositions of amino acid residues and secondary structures, and side-chain orientations. It was found that buried interface residues are important for recognizing multiple partners, while exposed interface residues are important for determining specificity to a particular ligand. In addition, our analyses reveal that residue compositions in OV and Non-OV regions are different and that residues in OV region show diverse side-chain torsion angles to accommodate binding to multiple targets.

To elucidate the structural basis of the diversity and universality in protein-protein interactions, an exhaustive all-against-all structural comparison of all known protein interfaces in the Protein Data Bank was performed at atomic resolution. After similar interfaces were clustered, approximately 20,000 structural motifs with at least two members were identified, out of which 3678 motifs consisted of at least 10 interfaces. Except for some trivial interfaces involving single alpha helices, almost all motifs were found to be confined within single protein families. Furthermore, the interaction partners of each motif were found to be very limited, and, accordingly, the interaction networks of the motifs tend to be small and are much more restricted than the binding sites for small ligand molecules. These findings suggest that, at the level of atomic structures, protein-protein interactions are precisely designed; hence, protein interfaces with multiple interacting partners should involve incompletely overlapping multiple interfaces and/or accommodate structural changes upon binding to their targets.

A thermionic RF gun usually suffers from a rapid decrease of electron beam energy caused by a rapid increase of beam current due to the back-bombardment effect. The authors propose a new method named ldquocavity detuningrdquo to keep the electron beam energy constant. This method detunes the resonant frequency of the gun cavity to a frequency that is lower (by several hundred kilohertz) than the driving frequency. A proof-of-principle experiment was carried out using a 4.5 cell S-band RF gun. The beam energy was successfully kept constant during a macro-pulse duration of 7.5 mus using a cavity detuning of -590 kHz, even with a significant current increase from 240 to 660 mA. A numerical simulation was also conducted to evaluate the cathode temperature, the current density on the cathode surface, and the phase stability of the output beam. All results demonstrated that cavity detuning provides an easy and simple way to compensate for the energy decrease in thermionic RF guns.

Comprehensive knowledge of protein-ligand interactions should provide a useful basis for annotating protein functions, studying protein evolution, engineering enzymatic activity, and designing drugs. To investigate the diversity and universality of ligand binding sites in protein structures, we conducted the all-against-all atomic-level structural comparison of over 180,000 ligand binding sites found in all the known structures in the Protein Data Bank by using a recently developed database search and alignment algorithm. By applying a hybrid top-down-bottom-up clustering analysis to the comparison results, we determined approximately 3000 well-defined structural motifs of ligand binding sites. Apart from a handful of exceptions, most structural motifs were found to be confined within single families or superfamilies, and to be associated with particular ligands. Furthermore, we analyzed the components of the similarity network and enumerated more than 4000 pairs of ligand binding sites that were shared across different protein folds., Comment: 13 pages, 8 figures

A method to search for local structural similarities in proteins at atomic resolution is presented. It is demonstrated that a huge amount of structural data can be handled within a reasonable CPU time by using a conventional relational database management system with appropriate indexing of geometric data. This method, which we call geometric indexing, can enumerate ligand binding sites that are structurally similar to sub-structures of a query protein among more than 160,000 possible candidates within a few hours of CPU time on an ordinary desktop computer. After detecting a set of high scoring ligand binding sites by the geometric indexing search, structural alignments at atomic resolution are constructed by iteratively applying the Hungarian algorithm, and the statistical significance of the final score is estimated from an empirical model based on a gamma distribution. Applications of this method to several protein structures clearly shows that significant similarities can be detected between local structures of non-homologous as well as homologous proteins., Comment: 29 pages, 8 figures, 3 tables

Background The nuclear magnetic resonance (NMR) spectroscopic data for biological macromolecules archived at the BioMagResBank (BMRB) provide a rich resource of biophysical information at atomic resolution. The NMR data archived in NMR-STAR ASCII format have been implemented in a relational database. However, it is still fairly difficult for users to retrieve data from the NMR-STAR files or the relational database in association with data from other biological databases. Findings To enhance the interoperability of the BMRB database, we present a full conversion of BMRB entries to two standard structured data formats, XML and RDF, as common open representations of the NMR-STAR data. Moreover, a SPARQL endpoint has been deployed. The described case study demonstrates that a simple query of the SPARQL endpoints of the BMRB, UniProt, and Online Mendelian Inheritance in Man (OMIM), can be used in NMR and structure-based analysis of proteins combined with information of single nucleotide polymorphisms (SNPs) and their phenotypes. Conclusions We have developed BMRB/XML and BMRB/RDF and demonstrate their use in performing a federated SPARQL query linking the BMRB to other databases through standard semantic web technologies. This will facilitate data exchange across diverse information resources. Electronic supplementary material The online version of this article (doi:10.1186/s13326-016-0057-1) contains supplementary material, which is available to authorized users.

Human transcriptional regulation factors, such as activators, repressors, and enhancer-binding factors are quite different from their prokaryotic counterparts in two respects: the average sequence in human is more than twice as long as that in prokaryotes, while the fraction of sequence aligned to domains of known structure is 31% in human transcription factors (TFs), less than half of that in bacterial TFs (72%). Intrinsically disordered (ID) regions were identified by a disorder-prediction program, and were found to be in good agreement with available experimental data. Analysis of 401 human TFs with experimental evidence from the Swiss-Prot database showed that as high as 49% of the entire sequence of human TFs is occupied by ID regions. More than half of the human TFs consist of a small DNA binding domain (DBD) and long ID regions frequently sandwiching unassigned regions. The remaining TFs have structural domains in addition to DBDs and ID regions. Experimental studies, particularly those with NMR, revealed that the transactivation domains in unbound TFs are usually unstructured, but become structured upon binding to their partners. The sequences of human and mouse TF orthologues are 90.5% identical despite a high incidence of ID regions, probably reflecting important functional roles played by ID regions. In general ID regions occupy a high fraction in TFs of eukaryotes, but not in prokaryotes. Implications of this dichotomy are discussed in connection with their functional roles in transcriptional regulation and evolution.

The living cell is inherently crowded with proteins and macromolecules. To avoid aggregation of denatured proteins in the living cell, molecular chaperones play important roles. Here we introduce a simple model to describe crowded protein solutions with chaperone-like species based on a dynamic density functional theory. As predicted by others, our simulations show that macromolecular crowding enhances the association of proteins and chaperones. However, when the intrinsic folding rate of the protein is slow, it is possible that crowding also enhances aggregation of proteins. The results of simulation suggest that, when the concentration of the crowding agent is as high as that in the cell, the association of the protein and unbound chaperone becomes correlated with the aggregation process, and that the protein-bound chaperones efficiently destroy the potential nuclei of aggregates and thus prevent the aggregation.

The direct-coupling analysis is a powerful method for protein contact prediction, and enables us to extract "direct" correlations between distant sites that are latent in "indirect" correlations observed in a protein multiple-sequence alignment. I show that the direct correlation can be obtained by using a formulation analogous to the Ornstein-Zernike integral equation in liquid theory. This formulation intuitively illustrates how the indirect or apparent correlation arises from an infinite series of direct correlations, and provides interesting insights into protein structure prediction., Comment: 3 pages, 1 figure

Background Linked Data has gained some attention recently in the life sciences as an effective way to provide and share data. As a part of the Semantic Web, data are linked so that a person or machine can explore the web of data. Resource Description Framework (RDF) is the standard means of implementing Linked Data. In the process of generating RDF data, not only are data simply linked to one another, the links themselves are characterized by ontologies, thereby allowing the types of links to be distinguished. Although there is a high labor cost to define an ontology for data providers, the merit lies in the higher level of interoperability with data analysis and visualization software. This increase in interoperability facilitates the multi-faceted retrieval of data, and the appropriate data can be quickly extracted and visualized. Such retrieval is usually performed using the SPARQL (SPARQL Protocol and RDF Query Language) query language, which is used to query RDF data stores. For the database provider, such interoperability will surely lead to an increase in the number of users. Results This manuscript describes the experiences and discussions shared among participants of the week-long BioHackathon 2011 who went through the development of RDF representations of their own data and developed specific RDF and SPARQL use cases. Advice regarding considerations to take when developing RDF representations of their data are provided for bioinformaticians considering making data available and interoperable. Conclusions Participants of the BioHackathon 2011 were able to produce RDF representations of their data and gain a better understanding of the requirements for producing such data in a period of just five days. We summarize the work accomplished with the hope that it will be useful for researchers involved in developing laboratory databases or data analysis, and those who are considering such technologies as RDF and Linked Data.