Quantification of the elevated rate of domain rearrangements in metazoa

Most eukaryotic proteins consist of multiple domains created through gene fusions or internal duplications. The most frequent change of a domain architecture (DA) is insertion or deletion of a domain at the N or C terminus. Still, the mechanisms underlying the evolution of multidomain proteins are not very well studied.Here, we have studied the evolution of multidomain architectures (MDA), guided by evolutionary information in the form of a phylogenetic tree. Our results show that Pfam domain families and MDAs have been created with comparable rates (0.1-1 per million years (My)). The major changes in DA evolution have occurred in the process of multicellularization and within the metazoan lineage. In contrast, creation of domains seems to have been frequent already in the early evolution. Furthermore, most of the architectures have been created from older domains or architectures, whereas novel domains are mainly found in single-domain proteins. However, a particular group of exon-bordering domains may have contributed to the rapid evolution of novel multidomain proteins in metazoan organisms. Finally, MDAs have evolved predominantly through insertions of domains, whereas domain deletions are less common.In conclusion, the rate of creation of multidomain proteins has accelerated in the metazoan lineage, which may partly be explained by the frequent insertion of exon-bordering domains into new architectures. However, our results indicate that other factors have contributed as well.     

Intra-chain 3D segment swapping spawns the evolution of new multidomain protein architectures

Multidomain proteins form in evolution through the concatenation of domains, but structural domains may comprise multiple segments of the chain. In this work, we demonstrate that new multidomain architectures can evolve by an apparent three-dimensional swap of segments between structurally similar domains within a single-chain monomer. By a comprehensive structural search of the current Protein Data Bank (PDB), we identified 32 well-defined segment-swapped proteins (SSPs) belonging to 18 structural families. Nearly 13% of all multidomain proteins in the PDB may have a segment-swapped evolutionary precursor as estimated by more permissive searching criteria. The formation of SSPs can be explained by two principal evolutionary mechanisms: (i) domain swapping and fusion (DSF) and (ii) circular permutation (CP). By large-scale comparative analyses using structural alignment and hidden Markov model methods, it was found that the majority of SSPs have evolved via the DSF mechanism, and a much smaller fraction, via CP. Functional analyses further revealed that segment swapping, which results in two linkers connecting the domains, may impart directed flexibility to multidomain proteins and contributes to the development of new functions. Thus, inter-domain segment swapping represents a novel general mechanism by which new protein folds and multidomain architectures arise in evolution, and SSPs have structural and functional properties that make them worth defining as a separate group.     

A simple model for proteins with interacting domains. Applications to scanning calorimetry data

A simple thermodynamic model is formulated for the purpose of interpreting scanning calorimetry data on proteins that have interacting domains. Interactions are quantified by inclusion of an interface free energy, delta GAB, in the thermodynamics of unfolding for multidomain proteins. The assumption is made that delta GAB goes to zero with the unfolding of either domain involved in pairwise interaction, so the interaction term appears to stabilize only the domain with the lower TM. Application of the model to calorimetric data leads to an estimate of -25,000 cal/mol for interactions between the regulatory and catalytic subunits of native aspartate transcarbamoylase and to a value of 0 for delta GAB between the transmembrane and cytoplasmic domains of band 3 of the human erythrocyte membrane. Estimates of changes in delta GAB are also obtained for mutant forms of yeast phosphoglycerate kinase that have been altered in the hinge region between amino-terminal and carboxy-terminal domains. The model is also applied to ligand binding to proteins having domains that communicate through pairwise interaction. It is shown that whenever the delta GAB term is ligand-dependent, then attachment of the ligand to the binding domain will be partially controlled by the other (regulatory) domain. This situation can sometimes be recognized and quantified when calorimetric scans are carried out at varying ligand concentrations. According to the model, the binding of MgATP to the carboxy-terminal domain of phosphoglycerate kinase is strongly stabilized (ca. 20% of the unitary free energy of binding) by participation of the amino-terminal domain, which acts to increase the binding constant 25-fold.(ABSTRACT TRUNCATED AT 250 WORDS)     

N-terminal domains of native multidomain proteins have the potential to assist de novo folding of their downstream domains in vivo by acting as solubility enhancers

The fusion of soluble partner to the N terminus of aggregation-prone polypeptide has been popularly used to overcome the formation of inclusion bodies in the E. coli cytosol. The chaperone-like functions of the upstream fusion partner in the artificial multidomain proteins could occur in de novo folding of native multidomain proteins. Here, we show that the N-terminal domains of three E. coli multidomain proteins such as lysyl-tRNA synthetase, threonyl-tRNA synthetase, and aconitase are potent solubility enhancers for various C-terminal heterologous proteins. The results suggest that the N-terminal domains could act as solubility enhancers for the folding of their authentic C-terminal domains in vivo. Tandem repeat of N-terminal domain or insertion of aspartic residues at the C terminus of the N-terminal domain also increased the solubility of fusion proteins, suggesting that the solubilizing ability correlates with the size and charge of N-terminal domains. The solubilizing ability of N-terminal domains would contribute to the autonomous folding of multidomain proteins in vivo, and based on these results, we propose a model of how N-terminal domains solubilize their downstream domains.     

Dynamic protein domains: identification, interdependence, and stability

Existing methods of domain identification in proteins usually provide no information about the degree of domain independence and stability. However, this information is vital for many areas of protein research. The recently developed hierarchical clustering of correlation patterns (HCCP) technique provides machine-based domain identification in a computationally simple and physically consistent way. Here we present the modification of this technique, which not only allows determination of the most plausible number of dynamic domains but also makes it possible to estimate the degree of their independence (the extent of correlated motion) and stability (the range of environmental conditions, where domains remain intact). With this technique we provided domain assignments and calculated intra- and interdomain correlations and interdomain energies for >2500 test proteins. It is shown that mean intradomain correlation of motions can serve as a quantitative criterion of domain independence, and the HCCP stability gap is a measure of their stability. Our data show that the motions of domains with high stability are usually independent. In contrast, the domains with moderate stability usually exhibit a substantial degree of correlated motions. It is shown that in multidomain proteins the domains are most stable if they are of similar size, and this correlates with the observed abundance of such proteins.     

The generation of new protein functions by the combination of domains

During evolution, many new proteins have been formed by the process of gene duplication and combination. The genes involved in this process usually code for whole domains. Small proteins contain one domain; medium and large proteins contain two or more domains. We have compared homologous domains that occur in both one-domain proteins and multidomain proteins. We have determined (1) how the functions of the individual domains in the multidomain proteins combine to produce their overall functions and (2) the extent to which these functions are similar to those in the one-domain homologs. We describe how domain combinations increase the specificity of enzymes; act as links between domains that have functional roles; regulate activity; combine within one chain functions that can act either independently, in concert or in new contexts; and provide the structural framework for the evolution of entirely new functions.     

General trends in the utilization of structural factors contributing to biological complexity

During evolution, proteins containing newly emerged domains and the increasing proportion of multidomain proteins in the full genome-encoded proteome (GEP) have substantially contributed to increasing biological complexity. However, it is not known how these two potential structural factors are preferentially utilized at given physiological states. Here, we classified proteins according to domain number and domain age and explored the general trends across species for the utilization of proteins from GEP to various certain-state proteomes (CSPs, i.e., all the proteins expressed at certain physiological states). We found that multidomain proteins or only older domain-containing proteins are significantly overrepresented in CSPs compared with GEP, which is a trend that is stronger in multicellular organisms than in unicellular organisms. Interestingly, the strengths of overrepresentation decreased during evolution of multicellular eukaryotes. When comparing across CSPs, we found that multidomain proteins are more overrepresented in complex tissues than in simpler ones, whereas no difference among proteins with domains of different ages is evident between complex and simple tissues. Thus, biological complexity under certain conditions is more significantly realized by diverse domain organization than by the emergence of new types of domain. In addition, we found that multidomain or only older domain-containing proteins tend to evolve slowly and generally are under stronger purifying selection, which may partly result from their general overrepresentation trends in CSPs.     

Accurate prediction of interfacial residues in two‐domain proteins using evolutionary information: Implications for three‐dimensional modeling

With the preponderance of multidomain proteins in eukaryotic genomes, it is essential to recognize the constituent domains and their functions. Often function involves communications across the domain interfaces, and the knowledge of the interacting sites is essential to our understanding of the structure–function relationship. Using evolutionary information extracted from homologous domains in at least two diverse domain architectures (single and multidomain), we predict the interface residues corresponding to domains from the two‐domain proteins. We also use information from the three‐dimensional structures of individual domains of two‐domain proteins to train naïve Bayes classifier model to predict the interfacial residues. Our predictions are highly accurate (～85%) and specific (～95%) to the domain–domain interfaces. This method is specific to multidomain proteins which contain domains in at least more than one protein architectural context. Using predicted residues to constrain domain–domain interaction, rigid‐body docking was able to provide us with accurate full‐length protein structures with correct orientation of domains. We believe that these results can be of considerable interest toward rational protein and interaction design, apart from providing us with valuable information on the nature of interactions. Proteins 2014; 82:1219–1234. © 2013 Wiley Periodicals, Inc.     

SAXSDom: Modeling multidomain protein structures using small-angle X-ray scattering data

Many proteins are composed of several domains that pack together into a complex tertiary structure. Multidomain proteins can be challenging for protein structure modeling, particularly those for which templates can be found for individual domains but not for the entire sequence. In such cases, homology modeling can generate high quality models of the domains but not for the orientations between domains. Small-angle X-ray scattering (SAXS) reports the structural properties of entire proteins and has the potential for guiding homology modeling of multidomain proteins. In this article, we describe a novel multidomain protein assembly modeling method, SAXSDom that integrates experimental knowledge from SAXS with probabilistic Input-Output Hidden Markov model to assemble the structures of individual domains together. Four SAXS-based scoring functions were developed and tested, and the method was evaluated on multidomain proteins from two public datasets. Incorporation of SAXS information improved the accuracy of domain assembly for 40 out of 46 critical assessment of protein structure prediction multidomain protein targets and 45 out of 73 multidomain protein targets from the ab initio domain assembly dataset. The results demonstrate that SAXS data can provide useful information to improve the accuracy of domain-domain assembly. The source code and tool packages are available at https://github.com/jianlin-cheng/SAXSDom .     

MOSC domains: ancient, predicted sulfur-carrier domains, present in diverse metal-sulfur cluster biosynthesis proteins including Molybdenum cofactor sulfurases

Using computational analysis, a novel superfamily of beta-strand-rich domains was identified in the Molybdenum cofactor sulfurase and several other proteins from both prokaryotes and eukaryotes. These MOSC domains contain an absolutely conserved cysteine and occur either as stand-alone forms such as the bacterial YiiM proteins, or fused to other domains such as a NifS-like catalytic domain in Molybdenum cofactor sulfurase. The MOSC domain is predicted to be a sulfur-carrier domain that receives sulfur abstracted by the pyridoxal phosphate-dependent NifS-like enzymes, on its conserved cysteine, and delivers it for the formation of diverse sulfur-metal clusters. The identification of this domain may clarify the mechanism of biogenesis of various metallo-enzymes including Molybdenum cofactor-containing enzymes that are compromised in human type II xanthinuria.     

Complex folding kinetics of a multidomain protein

Spectrin domains are three-helix bundles, commonly found in large tandem arrays. Equilibrium studies have shown that spectrin domains are significantly stabilized by their neighbors. In this work we show that domain:domain interactions can also have profound effects on their kinetic behavior. We have studied the folding of a tandem pair of spectrin domains (R1617) using a combination of single- and double-jump stopped flow experiments (monitoring folding by both circular dichroism and fluorescence). Mutant proteins were also used to investigate the complex folding kinetics. We find that, although the domains fold and unfold individually, there is a single rate-determining step for both folding and unfolding of the protein. This is consistent with the equilibrium observation of cooperative folding of the entire two-domain protein. The results may have important biological implications. Not only will the protein fold more efficiently during cotranslational folding, but the ability of the multidomain protein to withstand thermal unfolding in the cell will be dramatically increased. This study suggests that caution has to be exercised when extrapolating from single domains to larger proteins with a number of independently folding modules arranged in tandem. The multidomain protein spectrin is certainly more than "the sum of its parts".     

A New Protein Superfamily: TPPP-Like Proteins

The introduction of the term ‘Tubulin Polymerization Promoting Protein (TPPP)-like proteins’ is suggested. They constitute a eukaryotic protein superfamily, characterized by the presence of the p25alpha domain (Pfam05517, IPR008907), and named after the first identified member, TPPP/p25, exhibiting microtubule stabilizing function. TPPP-like proteins can be grouped on the basis of two characteristics: the length of their p25alpha domain, which can be long, short, truncated or partial, and the presence or absence of additional domain(s). TPPPs, in the strict sense, contain no other domains but one long or short p25alpha one (long- and short-type TPPPs, respectively). Proteins possessing truncated p25alpha domain are first described in this paper. They evolved from the long-type TPPPs and can be considered as arthropod-specific paralogs of long-type TPPPs. Phylogenetic analysis shows that the two groups (long-type and truncated TPPPs) split in the common ancestor of arthropods. Incomplete p25alpha domains can be found in multidomain TPPP-like proteins as well. The various subfamilies occur with a characteristic phyletic distribution: e. g., animal genomes/proteomes contain almost without exception long-type TPPPs; the multidomain apicortins occur almost exclusively in apicomplexan parasites. There are no data about the physiological function of these proteins except two human long-type TPPP paralogs which are involved in developmental processes of the brain and the musculoskeletal system, respectively. I predict that the superfamily members containing long or partial p25alpha domain are often intrinsically disordered proteins, while those with short or truncated domain(s) are structurally ordered. Interestingly, members of this superfamily connected or maybe connected to diseases are intrinsically disordered proteins.     

β-Strand-mediated interactions of protein domains

Protein domains exist by themselves or in combination with other domains to form complex multidomain proteins. Defining domain boundaries in proteins is essential for understanding their evolution and function but is not trivial. More specifically, partitioning domains that interact by forming a single β-sheet is known to be particularly troublesome for automatic structure-based domain decomposition pipelines. Here, we study edge-to-edge β-strand interactions between domains in a protein chain, to help define the boundaries for some more difficult cases where a single β-sheet spanning over two domains gives an appearance of one. We give a number of examples where β-strands belonging to a single β-sheet do not belong to a single domain and highlight the difficulties of automatic domain parsers on these examples. This work can be used as a baseline for defining domain boundaries in homologous proteins or proteins with similar domain interactions in the future.     

Folding of multidomain proteins: Biophysical consequences of tethering even in apparently independent folding

Most eukaryotic and a substantial fraction of prokaryotic proteins are composed of more than one domain. The tethering of these evolutionary, structural, and functional units raises, among others, questions regarding the folding process of conjugated domains. Studying the folding of multidomain proteins in silico enables one to identify and isolate the tethering‐induced biophysical determinants that govern crosstalks generated between neighboring domains. For this purpose, we carried out coarse‐grained and atomistic molecular dynamics simulations of two two‐domain constructs from the immunoglobulin‐like β‐sandwich fold. Each of these was experimentally shown to behave as the “sum of its parts,” that is, the thermodynamic and kinetic folding behavior of the constituent domains of these constructs seems to occur independently, with the folding of each domain uncoupled from the folding of its partner in the two‐domain construct. We show that the properties of the individual domains can be significantly affected by conjugation to another domain. The tethering may be accompanied by stabilizing as well as destabilizing factors whose magnitude depends on the size of the interface, the length, and the flexibility of the linker, and the relative stability of the domains. Accordingly, the folding of a multidomain protein should not be viewed as the sum of the folding patterns of each of its parts, but rather, it involves abrogating several effects that lead to this outcome. An imbalance between these effects may result in either stabilization or destabilization owing to the tethering. Proteins 2012; © 2012 Wiley Periodicals, Inc.     

This Déjà Vu Feeling—Analysis of Multidomain Protein Evolution in Eukaryotic Genomes

Evolutionary innovation in eukaryotes and especially animals is at least partially driven by genome rearrangements and the resulting emergence of proteins with new domain combinations, and thus potentially novel functionality. Given the random nature of such rearrangements, one could expect that proteins with particularly useful multidomain combinations may have been rediscovered multiple times by parallel evolution. However, existing reports suggest a minimal role of this phenomenon in the overall evolution of eukaryotic proteomes. We assembled a collection of 172 complete eukaryotic genomes that is not only the largest, but also the most phylogenetically complete set of genomes analyzed so far. By employing a maximum parsimony approach to compare repertoires of Pfam domains and their combinations, we show that independent evolution of domain combinations is significantly more prevalent than previously thought. Our results indicate that about 25% of all currently observed domain combinations have evolved multiple times. Interestingly, this percentage is even higher for sets of domain combinations in individual species, with, for instance, 70% of the domain combinations found in the human genome having evolved independently at least once in other species. We also show that previous, much lower estimates of this rate are most likely due to the small number and biased phylogenetic distribution of the genomes analyzed. The process of independent emergence of identical domain combination is widespread, not limited to domains with specific functional categories. Besides data from large-scale analyses, we also present individual examples of independent domain combination evolution. The surprisingly large contribution of parallel evolution to the development of the domain combination repertoire in extant genomes has profound consequences for our understanding of the evolution of pathways and cellular processes in eukaryotes and for comparative functional genomics.     

Engineered multidomain cysteine protease inhibitors yield resistance against western flower thrips (Frankliniella occidentalis) in greenhouse trials

Western flower thrips, Frankliniella occidentalis (Pergande) (Thysanoptera: Thripidae), cause very large economic damage on a variety of field and greenhouse crops. In this study, plant resistance against thrips was introduced into transgenic potato plants through the expression of novel, custom-made, multidomain protease inhibitors. Representative classes of inhibitors of cysteine and aspartic proteases [kininogen domain 3 (K), stefin A (A), cystatin C (C), potato cystatin (P) and equistatin (EIM)] were fused into reading frames consisting of four (K-A-C-P) to five (EIM-K-A-C-P) proteins, and were shown to fold into functional inhibitors in the yeast Pichia pastoris. The multidomain proteins were expressed in potato and found to be more resistant to degradation by plant proteases than the individual domains. In a time span of 14-16 days, transgenic potato plants expressing EIMKACP and KACP at a similar concentration reduced the number of larvae and adults to less than 20% of the control. Leaf damage on protected plants was minimal. Engineered multidomain cysteine protease inhibitors thus provide a novel way of controlling western flower thrips in greenhouse and field crops, and open up possibilities for novel insect resistance applications in transgenic crops.     

Single and multiple CH (calponin homology) domain containing multidomain proteins in Arabidopsis and Saccharomyces: an inventory

Genes for individual domains such as CH, lim, ankyrin, PH and RhoGAP, IQ motif, Ig_FLMN, spectrin, and EF hand probably existed in early evolution before there were plants, fungi or animals so that when we examine multidomain proteins in Arabidopsis, Saccharomyces, Dictyostelium or Homo Sapiens we encounter various combinations of such domains. While all of these four species express Fimbrin and EB1, the lists of CH containing multidomain proteins, however, differ in number and in type for each of them. There was no further great increase in the number of new single domain proteins. Still many new multidomain genes evolved—but far more so in metazoans—than in plants or fungi. In both plants and fungi only singlet CH domains but no doublets (other than those forming the Fimbrin quadruplet) were incorporated. That is in these two branches one finds no alpha actinin, dystrophin or filamin even though the individual building blocks (i.e. domains such as spectrin or IG-FLMN) were available in Arabidopsis. Possibly transposons create new chimeric multidomain genes by mixing and matching genes or gene fragments.