1
|
Blacklock KM, Yang L, Mulligan VK, Khare SD. A computational method for the design of nested proteins by loop-directed domain insertion. Proteins 2018; 86:354-369. [PMID: 29250820 DOI: 10.1002/prot.25445] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2017] [Revised: 12/04/2017] [Accepted: 12/15/2017] [Indexed: 12/23/2022]
Abstract
The computational design of novel nested proteins-in which the primary structure of one protein domain (insert) is flanked by the primary structure segments of another (parent)-would enable the generation of multifunctional proteins. Here we present a new algorithm, called Loop-Directed Domain Insertion (LooDo), implemented within the Rosetta software suite, for the purpose of designing nested protein domain combinations connected by flexible linker regions. Conformational space for the insert domain is sampled using large libraries of linker fragments for linker-to-parent domain superimposition followed by insert-to-linker superimposition. The relative positioning of the two domains (treated as rigid bodies) is sampled efficiently by a grid-based, mutual placement compatibility search. The conformations of the loop residues, and the identities of loop as well as interface residues, are simultaneously optimized using a generalized kinematic loop closure algorithm and Rosetta EnzymeDesign, respectively, to minimize interface energy. The algorithm was found to consistently sample near-native conformations and interface sequences for a benchmark set of structurally similar but functionally divergent domain-inserted enzymes from the α/β hydrolase superfamily, and discriminates well between native and nonnative conformations and sequences, although loop conformations tended to deviate from the native conformations. Furthermore, in cross-domain placement tests, native insert-parent domain combinations were ranked as the best-scoring structures compared to nonnative domain combinations. This algorithm should be broadly applicable to the design of multi-domain protein complexes with any combination of inserted or tandem domain connections.
Collapse
Affiliation(s)
- Kristin M Blacklock
- Institute for Quantitative Biomedicine, Rutgers The State University of New Jersey, Piscataway, New Jersey.,Department of Chemistry and Chemical Biology, Rutgers The State University of New Jersey, Piscataway, New Jersey.,Center for Integrative Proteomics Research, Rutgers The State University of New Jersey, Piscataway, New Jersey
| | - Lu Yang
- Department of Chemistry and Chemical Biology, Rutgers The State University of New Jersey, Piscataway, New Jersey.,Center for Integrative Proteomics Research, Rutgers The State University of New Jersey, Piscataway, New Jersey
| | - Vikram K Mulligan
- Institute for Protein Design and Department of Biochemistry, University of Washington, Seattle, Washington
| | - Sagar D Khare
- Institute for Quantitative Biomedicine, Rutgers The State University of New Jersey, Piscataway, New Jersey.,Department of Chemistry and Chemical Biology, Rutgers The State University of New Jersey, Piscataway, New Jersey.,Center for Integrative Proteomics Research, Rutgers The State University of New Jersey, Piscataway, New Jersey
| |
Collapse
|
2
|
Abstract
To explore protein space from a global perspective, we consider 9,710 SCOP (Structural Classification of Proteins) domains with up to 70% sequence identity and present all similarities among them as networks: In the "domain network," nodes represent domains, and edges connect domains that share "motifs," i.e., significantly sized segments of similar sequence and structure. We explore the dependence of the network on the thresholds that define the evolutionary relatedness of the domains. At excessively strict thresholds the network falls apart completely; for very lax thresholds, there are network paths between virtually all domains. Interestingly, at intermediate thresholds the network constitutes two regions that can be described as "continuous" versus "discrete." The continuous region comprises a large connected component, dominated by domains with alternating alpha and beta elements, and the discrete region includes the rest of the domains in isolated islands, each generally corresponding to a fold. We also construct the "motif network," in which nodes represent recurring motifs, and edges connect motifs that appear in the same domain. This network also features a large and highly connected component of motifs that originate from domains with alternating alpha/beta elements (and some all-alpha domains), and smaller isolated islands. Indeed, the motif network suggests that nature reuses such motifs extensively. The networks suggest evolutionary paths between domains and give hints about protein evolution and the underlying biophysics. They provide natural means of organizing protein space, and could be useful for the development of strategies for protein search and design.
Collapse
|
3
|
Dimitriadis D, Koumandou VL, Trimpalis P, Kossida S. Protein functional links in Trypanosoma brucei, identified by gene fusion analysis. BMC Evol Biol 2011; 11:193. [PMID: 21729286 PMCID: PMC3155505 DOI: 10.1186/1471-2148-11-193] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2011] [Accepted: 07/05/2011] [Indexed: 01/03/2023] Open
Abstract
BACKGROUND Domain or gene fusion analysis is a bioinformatics method for detecting gene fusions in one organism by comparing its genome to that of other organisms. The occurrence of gene fusions suggests that the two original genes that participated in the fusion are functionally linked, i.e. their gene products interact either as part of a multi-subunit protein complex, or in a metabolic pathway. Gene fusion analysis has been used to identify protein functional links in prokaryotes as well as in eukaryotic model organisms, such as yeast and Drosophila. RESULTS In this study we have extended this approach to include a number of recently sequenced protists, four of which are pathogenic, to identify fusion linked proteins in Trypanosoma brucei, the causative agent of African sleeping sickness. We have also examined the evolution of the gene fusion events identified, to determine whether they can be attributed to fusion or fission, by looking at the conservation of the fused genes and of the individual component genes across the major eukaryotic and prokaryotic lineages. We find relatively limited occurrence of gene fusions/fissions within the protist lineages examined. Our results point to two trypanosome-specific gene fissions, which have recently been experimentally confirmed, one fusion involving proteins involved in the same metabolic pathway, as well as two novel putative functional links between fusion-linked protein pairs. CONCLUSIONS This is the first study of protein functional links in T. brucei identified by gene fusion analysis. We have used strict thresholds and only discuss results which are highly likely to be genuine and which either have already been or can be experimentally verified. We discuss the possible impact of the identification of these novel putative protein-protein interactions, to the development of new trypanosome therapeutic drugs.
Collapse
|
4
|
Durrens P, Nikolski M, Sherman D. Fusion and fission of genes define a metric between fungal genomes. PLoS Comput Biol 2008; 4:e1000200. [PMID: 18949021 PMCID: PMC2557144 DOI: 10.1371/journal.pcbi.1000200] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2008] [Accepted: 09/05/2008] [Indexed: 12/19/2022] Open
Abstract
Gene fusion and fission events are key mechanisms in the evolution of gene architecture, whose effects are visible in protein architecture when they occur in coding sequences. Until now, the detection of fusion and fission events has been performed at the level of protein sequences with a post facto removal of supernumerary links due to paralogy, and often did not include looking for events defined only in single genomes. We propose a method for the detection of these events, defined on groups of paralogs to compensate for the gene redundancy of eukaryotic genomes, and apply it to the proteomes of 12 fungal species. We collected an inventory of 1,680 elementary fusion and fission events. In half the cases, both composite and element genes are found in the same species. Per-species counts of events correlate with the species genome size, suggesting a random mechanism of occurrence. Some biological functions of the genes involved in fusion and fission events are slightly over- or under-represented. As already noted in previous studies, the genes involved in an event tend to belong to the same functional category. We inferred the position of each event in the evolution tree of the 12 fungal species. The event localization counts for all the segments of the tree provide a metric that depicts the “recombinational” phylogeny among fungi. A possible interpretation of this metric as distance in adaptation space is proposed. One consequence of genome remodelling in evolution is the modification of genes, either by fusion with other genes, or by fission into several parts. By tracking the mathematical relations between groups of similar genes, rather than between individual genes, we can paint a global picture of remodelling across many species simultaneously. The strengths of our method are that it allows us to include highly redundant eukaryote genomes, and that it avoids alignment artifacts by representing each group of similar genes by a mathematical model. Applying our method to a set of fungal genomes, we confirmed first that the number of fusion/fission events is correlated with genome size, second that the fusion to fission ratio favors fusions, third that the set of events is not saturated, and fourth that while genes assembled in a fusion tend to have the same biochemical function, there appears to be little bias for the functions that are involved. Indeed, fusion and fission events are landmarks of random remodelling, independent of mutation rate: they define a metric of “recombination distance.” This distance lets us build a genome evolution history of species and may well be a better measure than mutation distance of the process of adaptation.
Collapse
Affiliation(s)
- Pascal Durrens
- MAGNOME Team, INRIA Centre de Recherche Bordeaux- Sud-Ouest, Laboratoire Bordelais de Recherche en Informatique, UMR 5800 CNRS, Domaine Universitaire, Talence Cedex, France.
| | | | | |
Collapse
|
5
|
Kamburov A, Goldovsky L, Freilich S, Kapazoglou A, Kunin V, Enright AJ, Tsaftaris A, Ouzounis CA. Denoising inferred functional association networks obtained by gene fusion analysis. BMC Genomics 2007; 8:460. [PMID: 18081932 PMCID: PMC2248599 DOI: 10.1186/1471-2164-8-460] [Citation(s) in RCA: 16] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2006] [Accepted: 12/14/2007] [Indexed: 12/04/2022] Open
Abstract
Background Gene fusion detection – also known as the 'Rosetta Stone' method – involves the identification of fused composite genes in a set of reference genomes, which indicates potential interactions between its un-fused counterpart genes in query genomes. The precision of this method typically improves with an ever-increasing number of reference genomes. Results In order to explore the usefulness and scope of this approach for protein interaction prediction and generate a high-quality, non-redundant set of interacting pairs of proteins across a wide taxonomic range, we have exhaustively performed gene fusion analysis for 184 genomes using an efficient variant of a previously developed protocol. By analyzing interaction graphs and applying a threshold that limits the maximum number of possible interactions within the largest graph components, we show that we can reduce the number of implausible interactions due to the detection of promiscuous domains. With this generally applicable approach, we generate a robust set of over 2 million distinct and testable interactions encompassing 696,894 proteins in 184 species or strains, most of which have never been the subject of high-throughput experimental proteomics. We investigate the cumulative effect of increasing numbers of genomes on the fidelity and quantity of predictions, and show that, for large numbers of genomes, predictions do not become saturated but continue to grow linearly, for the majority of the species. We also examine the percentage of component (and composite) proteins with relation to the number of genes and further validate the functional categories that are highly represented in this robust set of detected genome-wide interactions. Conclusion We illustrate the phylogenetic and functional diversity of gene fusion events across genomes, and their usefulness for accurate prediction of protein interaction and function.
Collapse
Affiliation(s)
- Atanas Kamburov
- Computational Genomics Group, The European Bioinformatics Institute, EMBL Cambridge Outstation, Cambridge CB10 1SD, UK.
| | | | | | | | | | | | | | | |
Collapse
|
6
|
Abstract
Many essential cellular processes such as signal transduction, transport, cellular motion and most regulatory mechanisms are mediated by protein-protein interactions. In recent years, new experimental techniques have been developed to discover the protein-protein interaction networks of several organisms. However, the accuracy and coverage of these techniques have proven to be limited, and computational approaches remain essential both to assist in the design and validation of experimental studies and for the prediction of interaction partners and detailed structures of protein complexes. Here, we provide a critical overview of existing structure-independent and structure-based computational methods. Although these techniques have significantly advanced in the past few years, we find that most of them are still in their infancy. We also provide an overview of experimental techniques for the detection of protein-protein interactions. Although the developments are promising, false positive and false negative results are common, and reliable detection is possible only by taking a consensus of different experimental approaches. The shortcomings of experimental techniques affect both the further development and the fair evaluation of computational prediction methods. For an adequate comparative evaluation of prediction and high-throughput experimental methods, an appropriately large benchmark set of biophysically characterized protein complexes would be needed, but is sorely lacking.
Collapse
Affiliation(s)
- András Szilágyi
- Center of Excellence in Bioinformatics, University at Buffalo, State University of New York, 901 Washington St, Buffalo, NY 14203, USA
| | | | | | | |
Collapse
|
7
|
Lindner HA, Nadeau G, Matte A, Michel G, Ménard R, Cygler M. Site-directed Mutagenesis of the Active Site Region in the Quinate/Shikimate 5-Dehydrogenase YdiB of Escherichia coli. J Biol Chem 2005; 280:7162-9. [PMID: 15596430 DOI: 10.1074/jbc.m412028200] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022] Open
Abstract
YdiB and its paralog AroE are members of the quinate/shikimate 5-dehdrogenase family. Enzymes from this family function in the shikimate pathway that is essential for survival of microorganisms and plants and represent potential drug targets. Recent YdiB and AroE crystal structures revealed the presence of a NAD(P)-binding and a catalytic domain. We carried out site-directed mutagenesis of 8 putative active site residues in YdiB from Escherichia coli and analyzed structural and kinetic properties of the mutant enzymes. Our data indicate critical roles for an invariant lysine and aspartate residue in substrate binding and allowed us to differentiate between two previously proposed models for the binding of the substrate in the active site. Comparison of several YdiB and AroE structures led us to conclude that, upon cofactor binding and domain closure, the 2 identified binding residues are repositioned to bind to the substrate. Although the lysine residue contributes to some extent to the stabilization of the transition state, we did not identify any residue as catalytically essential. This indicates that catalysis does not operate through a general acid-base mechanism, as thought originally. Our improved understanding of the medically and agriculturally important quinate/shikimate 5-dehydrogenase family at the molecular level may prove useful in the development of novel herbicides and antimicrobial agents.
Collapse
Affiliation(s)
- Holger A Lindner
- Biotechnology Research Institute, National Research Council of Canada (NRCC), Montréal, Québec H4P 2R2, Canada.
| | | | | | | | | | | |
Collapse
|
8
|
Aroul-Selvam R, Hubbard T, Sasidharan R. Domain insertions in protein structures. J Mol Biol 2004; 338:633-41. [PMID: 15099733 PMCID: PMC2665287 DOI: 10.1016/j.jmb.2004.03.039] [Citation(s) in RCA: 53] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2003] [Revised: 03/07/2004] [Accepted: 03/10/2004] [Indexed: 10/26/2022]
Abstract
Domains are the structural, functional or evolutionary units of proteins. Proteins can comprise a single domain or a combination of domains. In multi-domain proteins, the domains almost always occur end-to-end, i.e., one domain follows the C-terminal end of another domain. However, there are exceptions to this common pattern, where multi-domain proteins are formed by insertion of one domain (insert) into another domain (parent). Here, we provide a quantitative description of known insertions in the Protein Data Bank (PDB). We found that 9% of domain combinations observed in non-redundant PDB are insertions. Although 90% of all insertions involve only one insert, proteins can clearly have multiple (nested, two-domain and three-domain) inserts. We also observed correlations between the structure and function of a domain and its tendency to be found as a parent or an insert. There is a bias in insert position towards the C terminus of parents. We observed that the atomic distance between the N and C terminus of an insert is significantly smaller when compared to the N-to-C distance in a parent context or a single domain context. Insertions are found always to occur in loop regions of parent domains. Our observations regarding the relationship between domain insertions and the structure, function and evolution of proteins have implications for protein engineering.
Collapse
Affiliation(s)
- R. Aroul-Selvam
- The Wellcome Trust Sanger Institute, Genome Campus Hinxton, Cambridge CB10 1SA UK
| | - Tim Hubbard
- The Wellcome Trust Sanger Institute, Genome Campus Hinxton, Cambridge CB10 1SA UK
| | - Rajkumar Sasidharan
- MRC Laboratory of Molecular Biology, Hills Road, Cambridge CB2 2QH, UK
- Corresponding author E-mail address of the corresponding author:
| |
Collapse
|