Functional fingerprints of folds: evidence for correlated structure-function evolution

Using structural similarity clustering of protein domains: protein domain universe graph (PDUG), and a hierarchical functional annotation: gene ontology (GO) as two evolutionary lenses, we find that each structural cluster (domain fold) exhibits a distribution of functions that is unique to it. These functional distributions are functional fingerprints that are specific to characteristic structural clusters and vary from cluster to cluster. Furthermore, as structural similarity threshold for domain clustering in the PDUG is relaxed we observe an influx of earlier-diverged domains into clusters. These domains join clusters without destroying the functional fingerprint. These results can be understood in light of a divergent evolution scenario that posits correlated divergence of structural and functional traits in protein domains from one or few progenitors.     

Classification and annotation of the relationship between protein structural change and ligand binding

The causal relationship between protein structural change and ligand binding was classified and annotated for 839 nonredundant pairs of crystal structures in the Protein Data Bank—one with and the other without a bound low-molecular-weight ligand molecule. Protein structural changes were first classified into either domain or local motions depending on the size of the moving protein segments. Whether the protein motion was coupled with ligand binding was then evaluated based on the location of the ligand binding site and by application of the linear response theory of protein structural change. Protein motions coupled with ligand binding were further classified into either closure or opening motions. This classification revealed the following: (i) domain motions coupled with ligand binding are dominated by closure motions, which can be described by the linear response theory; (ii) local motions frequently accompany order-disorder or α-helix-coil conformational transitions; and (iii) transferase activity (Enzyme Commission   number 2) is the predominant function among coupled domain closure motions. This could be explained by the closure motion acting to insulate the reaction site of these enzymes from environmental water.     

miRDB: a microRNA target prediction and functional annotation database with a wiki interface

MicroRNAs (miRNAs) are short noncoding RNAs that are involved in the regulation of thousands of gene targets. Recent studies indicate that miRNAs are likely to be master regulators of many important biological processes. Due to their functional importance, miRNAs are under intense study at present, and many studies have been published in recent years on miRNA functional characterization. The rapid accumulation of miRNA knowledge makes it challenging to properly organize and present miRNA function data. Although several miRNA functional databases have been developed recently, this remains a major bioinformatics challenge to miRNA research community. Here, we describe a new online database system, miRDB, on miRNA target prediction and functional annotation. Flexible web search interface was developed for the retrieval of target prediction results, which were generated with a new bioinformatics algorithm we developed recently. Unlike most other miRNA databases, miRNA functional annotations in miRDB are presented with a primary focus on mature miRNAs, which are the functional carriers of miRNA-mediated gene expression regulation. In addition, a wiki editing interface was established to allow anyone with Internet access to make contributions on miRNA functional annotation. This is a new attempt to develop an interactive community-annotated miRNA functional catalog. All data stored in miRDB are freely accessible at http://mirdb.org.     

Functional information in SWISS-PROT: the basis for large-scale characterisation of protein sequences

With the rapid growth of sequence databases, there is an increasing need for reliable functional characterisation and annotation of newly predicted proteins. To cope with such large data volumes, faster and more effective means of protein sequence characterisation and annotation are required. One promising approach is automatic large-scale functional characterisation and annotation, which is generated with limited human interaction. However, such an approach is heavily dependent on reliable data sources. The SWISS-PROT protein sequence database plays an essential role here owing to its high level of functional information.     

Folding and stability of a primitive protein

We have previously attempted to simulate domain creation in early protein evolution by recombining polypeptide segments from non-homologous proteins, and we have described the structure of one such de novo protein, 1b11, a segment-swapped tetramer with novel architecture. Here, we have analyzed the thermodynamic stability and folding kinetics of the 1b11 tetramer and its monomeric and dimeric intermediates, and of 1b11 mutants with changes at the domain interface. Denatured 1b11 polypeptides fold into transient, folded monomers with marginal stability (DeltaG<1kcalmol(-1)) which convert rapidly ( approximately 6x10(4)M(-1)s(-1)) into dimers (DeltaG=9.8kcal/mol) and then more slowly ( approximately 3M(-1)s(-1)) into tetramers (DeltaG=28kcalmol(-1)). Segment swapping takes place during dimerization, as suggested by mass spectroscopic analysis of covalently linked peptides derived from proteolysis of a disulfide-linked dimer. Our results confirm that segment swapping and associated oligomerization are both powerful ways of stabilizing proteins, and we suggest that this may have been a feature of early protein evolution.     

Protein import into mitochondria: origins and functions today (Review)

Mitochondria are organelles derived from α-proteobacteria over the course of one to two billion years. Mitochondria from the major eukaryotic lineages display some variation in functions and coding capacity but sequence analysis demonstrates them to be derived from a single common ancestral endosymbiont. The loss of assorted functions, the transfer of genes to the nucleus, and the acquisition of various ‘eukaryotic’ proteins have resulted in an organelle that contains approximately 1000 different proteins, with most of these proteins imported into the organelle across one or two membranes. A single translocase in the outer membrane and two translocases in the inner membrane mediate protein import. Comparative sequence analysis and functional complementation experiments suggest some components of the import pathways to be directly derived from the eubacterial endosymbiont's own proteins, and some to have arisen ‘de novo’ at the earliest stages of ‘mitochondrification’ of the endosymbiont. A third class of components appears lineage-specific, suggesting they were incorporated into the process of protein import long after mitochondria was established as an organelle and after the divergence of the various eukaryotic lineages. Protein sorting pathways inherited from the endosymbiont have been co-opted and play roles in intraorganelle protein sorting after import. The import apparatus of animals and fungi show significant similarity to one another, but vary considerably to the plant apparatus. Increasing complexity in the eukaryotic lineage, i.e., from single celled to multi-cellular life forms, has been accompanied by an expansion in genes encoding each component, resulting in small gene families encoding many components. The functional differences in these gene families remain to be elucidated, but point to a mosaic import apparatus that can be regulated by a variety of signals.     

Modeling the evolution of protein domain architectures using maximum parsimony

Domains are basic evolutionary units of proteins and most proteins have more than one domain. Advances in domain modeling and collection are making it possible to annotate a large fraction of known protein sequences by a linear ordering of their domains, yielding their architecture. Protein domain architectures link evolutionarily related proteins and underscore their shared functions. Here, we attempt to better understand this association by identifying the evolutionary pathways by which extant architectures may have evolved. We propose a model of evolution in which architectures arise through rearrangements of inferred precursor architectures and acquisition of new domains. These pathways are ranked using a parsimony principle, whereby scenarios requiring the fewest number of independent recombination events, namely fission and fusion operations, are assumed to be more likely. Using a data set of domain architectures present in 159 proteomes that represent all three major branches of the tree of life allows us to estimate the history of over 85% of all architectures in the sequence database. We find that the distribution of rearrangement classes is robust with respect to alternative parsimony rules for inferring the presence of precursor architectures in ancestral species. Analyzing the most parsimonious pathways, we find 87% of architectures to gain complexity over time through simple changes, among which fusion events account for 5.6 times as many architectures as fission. Our results may be used to compute domain architecture similarities, for example, based on the number of historical recombination events separating them. Domain architecture "neighbors" identified in this way may lead to new insights about the evolution of protein function.     

On adaptation: A reduction of the Kauffman-Levin model to a problem in graph theory and its consequences

It is shown that complex adaptations are best modelled as discrete processes represented on directed weighted graphs. Such a representation captures the idea that problems of adaptation in evolutionary biology are problems in a discrete space, something that the conventional representations using continuous adaptive landscapes does not. Further, this representation allows the utilization of well-known algorithms for the computation of several biologically interesting results such as the accessibility of one allele from another by a specified number of point mutations, the accessibility of alleles at a local maximum of fitness, the accessibility of the allele with the globally maximum fitness, etc. A reduction of a model due to Kauffman and Levin to such a representation is explicitly carried out and it is shown how this reduction clarifies the biological questions that are of interest.Thanks are due to William Wimsatt, James F. Crow, and the referees for Biology and Philosophy for comments on an earlier version of this paper. Remarks by members of the audience, especially Abner Shimony, of a seminar at Boston University, February 19, 1988, were also very helpful. The diagrams were prepared with the assistance of Tracy Lubas.     

Structure and functional annotation of hypothetical proteins having putative Rubisco activase function from Vitis vinifera

Rubisco is a very large, complex and one of the most abundant proteins in the world and comprises up to 50% of all soluble protein in plants. The activity of Rubisco, the enzyme that catalyzes CO₂ assimilation in photosynthesis, is regulated by Rubisco activase (Rca). In the present study, we searched for hypothetical protein of Vitis vinifera which has putative Rubisco activase function. The Arabidopsis and tobacco Rubisco activase protein sequences were used as seed sequences to search against Vitis vinifera in UniprotKB database. The selected hypothetical proteins of Vitis vinifera were subjected to sequence, structural and functional annotation. Subcellular localization predictions suggested it to be cytoplasmic protein. Homology modelling was used to define the three-dimensional (3D) structure of selected hypothetical proteins of Vitis vinifera. Template search revealed that all the hypothetical proteins share more than 80% sequence identity with structure of green-type Rubisco activase from tobacco, indicating proteins are evolutionary conserved. The homology modelling was generated using SWISS-MODEL. Several quality assessment and validation parameters computed indicated that homology models are reliable. Further, functional annotation through PFAM, CATH, SUPERFAMILY, CDART suggested that selected hypothetical proteins of Vitis vinifera contain ATPase family associated with various cellular activities (AAA) and belong to the AAA+ super family of ring-shaped P-loop containing nucleoside triphosphate hydrolases. This study will lead to research in the optimization of the functionality of Rubisco which has large implication in the improvement of plant productivity and resource use efficiency.     

Colicin import into E. coli cells: A model system for insights into the import mechanisms of bacteriocins

Bacteriocins are a diverse group of ribosomally synthesized protein antibiotics produced by most bacteria. They range from small lanthipeptides produced by lactic acid bacteria to much larger multi domain proteins of Gram negative bacteria such as the colicins from Escherichia coli. For activity bacteriocins must be released from the producing cell and then bind to the surface of a sensitive cell to instigate the import process leading to cell death. For over 50 years, colicins have provided a working platform for elucidating the structure/function studies of bacteriocin import and modes of action. An understanding of the processes that contribute to the delivery of a colicin molecule across two lipid membranes of the cell envelope has advanced our knowledge of protein–protein interactions (PPI), protein–lipid interactions and the role of order–disorder transitions of protein domains pertinent to protein transport. In this review, we provide an overview of the arrangement of genes that controls the synthesis and release of the mature protein. We examine the uptake processes of colicins from initial binding and sequestration of binding partners to crossing of the outer membrane, and then discuss the translocation of colicins through the cell periplasm and across the inner membrane to their cytotoxic site of action. This article is part of a Special Issue entitled: Protein trafficking and secretion in bacteria. Guest Editors: Anastassios Economou and Ross Dalbey.     

Quantifying the impact of protein tertiary structure on molecular evolution

To investigate the evolutionary impact of protein structure, the experimentally determined tertiary structure and the protein-coding DNA sequence were collected for each of 1,195 genes. These genes were studied via a model of sequence change that explicitly incorporates effects on evolutionary rates due to protein tertiary structure. In the model, these effects act via the solvent accessibility environments and pairwise amino acid interactions that are induced by tertiary structure. To compare the hypotheses that structure does and does not have a strong influence on evolution, Bayes factors were estimated for each of the 1,195 sequences. Most of the Bayes factors strongly support the hypothesis that protein structure affects protein evolution. Furthermore, both solvent accessibility and pairwise interactions among amino acids are inferred to have important roles in protein evolution. Our results also indicate that the strength of the relationship between tertiary structure and evolution has a weak but real correlation to the annotation information in the Gene Ontology database. Although their influences on rates of evolution vary among protein families, we find that the mean impacts of solvent accessibility and pairwise interactions are about the same.     

Unravelling the binding affinity between model transport protein and a prospective tuberculosis therapeutic agent: a spectroscopic and theoretical simulation exploration

Abstract

Haloxyfop was reported to exhibit inhibition effect targeting Mycobacterium tuberculosis and pathogenic parasites. To pave its way for drug development, more research is required to determine the affinities interacting with biological receptors in vivo. In this work, the interactions of Haloxyfop with two model transport proteins were investigated by spectroscopic techniques and theoretical simulation. The interaction mechanism, thermodynamic properties and the impact of Haloxyfop-induced conformational change in serum albumins were revealed by series of fluorescence, UV-Vis absorption and circular dichroic spectroscopy. The specific binding sites were determined by site-competitive replacement experiment. Molecular docking and dynamic simulation provided a visual screening in the microscopic binding mode. The structure of Haloxyfop was roughly divided into three parts that exhibit different covalent interaction affinities. The two isomers of Haloxyfop showed a certain degree of affinity difference. Hydrophobic, polar interaction and π-effect were analyzed in detail, and the surface electrostatic potential energy maps were simulated to provide references. The free energy, calculated by the molecular mechanics-generalized born surface area (MM-GBSA) and molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) methods, was decomposed to per-residues, which intuitively revealed relevant contributions in binding process. The role of water existence was explored through molecular dynamic refinement, and the frontier molecular orbital analysis explored the ionic interaction mechanism in electronic level. In general, multiple chemistry method was adopted to fully unravel the properties of Haloxyfop-binding for the sake of rationalizing the applicability as a therapeutic agent.

Communicated by Ramaswamy H. Sarma     

Molecular basis of the Tfs1/Ira2 interaction: a combined protein engineering and molecular modelling study

Tfs1p and Ylr179cp are yeast proteins belonging to the PEBP family. Tfs1p, but not Ylr179cp, has been shown to interact with and inhibit Ira2p, a GTPase-activating protein of Ras. Tfs1p has been shown to be a specific inhibitor of the CPY protease and the 3D structure of the complex has been resolved. To shed light on the molecular determinants of Tfs1p involved in the Tfs1/Ira2 interaction, the 3D structure of Ylr179cp has been modelled and compared to that of Tfs1p. Tfs1p point mutants and Tfs1 hybrid proteins combining regions of Tfs1p and Ylr179cp were also designed and their function was tested. Results, interpreted from a structural point of view, show that the accessibility of the surface pocket of Tfs1p, its N-terminal region and the specific electrostatic properties of a large surface region containing these two elements, play a crucial role in this interaction.     

The prediction of protein subcellular localization from sequence: a shortcut to functional genome annotation.

Automated sequence annotation is a major goal of post-genomic era with hundreds of genomes in the databases, from both prokaryotes and eukaryotes. While the number of fully sequenced chromosomes from microbial organisms exponentially increased in the last decade above 600, presently we know the whole DNA content of only 25 eukaryotic organisms, including Homo sapiens. However, the process of genome annotation is far from being completed. This is particularly relevant in eukaryotes, whose cells contain several subcellular compartments, or organelles, enclosed by membranes, where different relevant functions are performed. Translocation across the membrane into the organelles is a highly regulated and complex cellular process. Indeed different proteins and/or protein isoforms, originated from genes by alternative splicing, may be conveyed to different cell compartments, depending on their specific role in the cell. During recent years the prediction of subcellular localization (SL) by computational means has been an active research area. Several methods are presently available based on different notions and addressing different aspects of SL. This review provides a short overview of the most well performing methods described in the literature, highlighting their predictive capabilities and different applications.     

Sequencing,annotation and comparative analysis of nine BACs of the giant panda(Ailuropoda melanoleuca)

A 10-fold BAC library for the giant panda was constructed and nine BACs were selected to generate finish sequences.These BACs could be used as a validation resource for the de novo assembly accuracy of the whole genome shotgun sequencing reads of the giant panda newly generated by Illumina GA sequencing technology.Complete Sanger sequencing,assembly,annotation and comparative analysis were carried out on the selected BACs of a joint length 878 kb.Homologue search and de novo prediction methods were used to ...     

Evolutionary graph theory: breaking the symmetry between interaction and replacement

We study evolutionary dynamics in a population whose structure is given by two graphs: the interaction graph determines who plays with whom in an evolutionary game; the replacement graph specifies the geometry of evolutionary competition and updating. First, we calculate the fixation probabilities of frequency dependent selection between two strategies or phenotypes. We consider three different update mechanisms: birth-death, death-birth and imitation. Then, as a particular example, we explore the evolution of cooperation. Suppose the interaction graph is a regular graph of degree h, the replacement graph is a regular graph of degree g and the overlap between the two graphs is a regular graph of degree l. We show that cooperation is favored by natural selection if b/c>hg/l. Here, b and c denote the benefit and cost of the altruistic act. This result holds for death-birth updating, weak-selection and large population size. Note that the optimum population structure for cooperators is given by maximum overlap between the interaction and the replacement graph (g=h=l), which means that the two graphs are identical. We also prove that a modified replicator equation can describe how the expected values of the frequencies of an arbitrary number of strategies change on replacement and interaction graphs: the two graphs induce a transformation of the payoff matrix.     

Regulators of G protein signaling: a bestiary of modular protein binding domains

Members of the newly discovered regulator of G protein signaling (RGS) families of proteins have a common RGS domain. This RGS domain is necessary for conferring upon RGS proteins the capacity to regulate negatively a variety of Galpha protein subunits. However, RGS proteins are more than simply negative regulators of signaling. RGS proteins can function as effector antagonists, and recent evidence suggests that RGS proteins can have positive effects on signaling as well. Many RGS proteins possess additional C- and N-terminal modular protein-binding domains and motifs. The presence of these additional modules within the RGS proteins provides for multiple novel regulatory interactions performed by these molecules. These regions are involved in conferring regulatory selectivity to specific Galpha-coupled signaling pathways, enhancing the efficacy of the RGS domain, and the translocation or targeting of RGS proteins to intracellular membranes. In other instances, these domains are involved in cross-talk between different Galpha-coupled signaling pathways and, in some cases, likely serve to integrate small GTPases with these G protein signaling pathways. This review discusses these C- and N-terminal domains and their roles in the biology of the brain-enriched RGS proteins. Methods that can be used to investigate the function of these domains are also discussed.