- Open Access
Domain atrophy creates rare cases of functional partial protein domains
Genome Biology volume 16, Article number: 88 (2015)
Protein domains display a range of structural diversity, with numerous additions and deletions of secondary structural elements between related domains. We have observed a small number of cases of surprising large-scale deletions of core elements of structural domains. We propose a new concept called domain atrophy, where protein domains lose a significant number of core structural elements.
Here, we implement a new pipeline to systematically identify new cases of domain atrophy across all known protein sequences. The output of this pipeline was carefully checked by hand, which filtered out partial domain instances that were unlikely to represent true domain atrophy due to misannotations or un-annotated sequence fragments. We identify 75 cases of domain atrophy, of which eight cases are found in a three-dimensional protein structure and 67 cases have been inferred based on mapping to a known homologous structure. Domains with structural variations include ancient folds such as the TIM-barrel and Rossmann folds. Most of these domains are observed to show structural loss that does not affect their functional sites.
Our analysis has significantly increased the known cases of domain atrophy. We discuss specific instances of domain atrophy and see that there has often been a compensatory mechanism that helps to maintain the stability of the partial domain. Our study indicates that although domain atrophy is an extremely rare phenomenon, protein domains under certain circumstances can tolerate extreme mutations giving rise to partial, but functional, domains.
Protein domains are independently evolving structural and function units. These units have been recombined in many ways creating proteins that perform diverse functions. Proteins, in the course of evolution, may acquire (domain gain) or lose (domain loss) one or more entire domains by a variety of means [1,2]. At a smaller scale of modification, protein domains can gain or lose secondary structural elements through insertion or deletion of amino acid residues. Insertions or deletions are known to significantly influence protein folding  and are considered important evolutionary factors which influence protein structure and function . Proteins and their domains are only marginally stable, such that the mutation of a single residue may be sufficient to induce the unfolding of the domain. Single point mutations can lead to a large range of human diseases. However, through evolution a protein’s sequence is slowly changed so that quite divergent sequences can give rise to a similar three-dimensional (3D) structure. This means that proteins are tolerant of small stepwise changes in their sequence.
Analyses of a large number of domain structures have indicated significant variability in domain lengths . These variations are mostly attributed to insertions or deletions in loops, coils or of a few secondary structural elements, leaving the domain core largely unaffected [5-7]. Cases of such domain length variations caused by gain of accessory secondary structural elements or ‘embellishments’ are well documented. Examples of the so-called domain elaborations can be observed in the HUP superfamily , wherein embellishments influence interactions, substrate binding and stability of the core domain.
Similar to domain elaborations, one could envisage events leading to loss or degradation of a domain’s secondary structure elements necessitated by similar evolutionary constraints as above. We propose the term ‘domain atrophy’ for events that lead to significant loss of core secondary structure elements (Figure 1). Such loss of domain-structural elements has been observed in the ‘truncated globins’ , bacterial luciferases  and recently in a domain of unknown function DUF2172 .
It seems possible that the process of domain atrophy occurs by a single large-scale mutation, which raises the question of how can a protein tolerate such a large destabilising event. Such events of loss of >10 amino acid residues are more likely to be very destabilising. Losses at protein termini are more likely to be tolerated, particularly if they are of peripheral secondary structures. For example, partial deletion of the first seven residues at the N-terminus of Ig-domain of the human muscle protein titin only marginally decreased its stability by 2.8 kcal/mol , however mutations within the core are unlikely to be tolerated. Deletions of many amino acids are normally unlikely to be tolerated unless they are found in a loop region. Therefore, natural occurrences of stable protein domains with large-scale deletions of core secondary structures [10,11] are intriguing and suggest that compensatory mechanisms that help stabilise these atrophied domains must exist.
An alternative possibility is that a globular domain that undergoes an atrophy event then becomes disordered and may adopt some new function. We believe that we are unlikely to be able to observe such events as the pattern of conservation is likely to significantly change leading to very rapid sequence divergence. However, the relationship of domain atrophy to protein disorder may be an interesting avenue for future research.
Unlike domain elaboration, which is commonly seen in protein structures, structural data and literature on domain atrophy is very scarce. Using sequence based protein-HMM models of domain families, we have devised an algorithm to identify potential cases of domain atrophy and have identified several new cases in proteins of known structure. Where 3D structures of the atrophied domain were not available, we have used sequence mapping onto reference complete domain structures to infer potential domain atrophy events.
A companion paper by Triant and Pearson  describes several computational artefacts that can produce partial domain instances, which may confound the discovery of cases of domain atrophy. They show that domains can appear short because they have been broken up into pieces by the alignment process, or because they are annotated on incorrectly assembled proteins, or because the Pfam domain models have been built from two or more evolutionarily mobile domains. In this paper, we use several strategies to avoid these biological and computational artefacts, and identify a set of atrophied domains that are complementary to the ‘structural partials’ identified by Triant and Pearson.
Results and discussion
In order to quantify the magnitude of structural loss we have defined a metric called the Atrophy Score (see Methods), which indicates the fraction of the domain that is likely to be missing. It is based on looking at matches of Pfam profile-HMMs to a sequence. For some domains only part of the profile-HMM model may match to the sequence. We use these partial matches to identify where a protein domain is potentially atrophied. The atrophy score is equivalent to the fraction of a Pfam profile-HMM model that is missing from a protein domain. While a positive atrophy score indicates structural loss, zero or negative atrophy scores indicate complete domain boundaries, that is, no atrophy. We did not investigate instances with an atrophy score below 0.15, where less than 15% of the domain is lost. Domains showing this small level of atrophy are very likely to be enriched in false positives where loss of peripheral secondary structures can be tolerated. A feature of the atrophy score is the possibility of scores to reach 1 or higher. This can happen when one domain is nested within another, leading to a negative inter-domain distance. These cases are excluded from our analysis.
A partial match to a Pfam profile-HMM is only weak evidence for domain atrophy. As Triant and Pearson show, even if a match to the whole HMM is not present one can often find that the similarity can be extended to the full domain length with a sequence similarity comparison. Therefore, we look in particular at partial domains that are end bounded by either another domain, or the terminus of the sequence.
We have categorised domain atrophy into five different types that are fully described in the Methods section (Figure 2). The five types are: (1) N-terminal end-bounded atrophy; (2) C-terminal end-bounded atrophy; (3) Upstream domain-bounded atrophy; (4) Downstream domain-bounded atrophy; and (v) Within-domain atrophy.
Domain instances that scored between 0.15 and 1 were manually examined to identify potential false positive cases (see Table 1). Manual inspection helped to discern probable true cases of domain atrophy from possible false positive cases. We were able to confirm eight cases, which have a known structure and a further 67 cases that belonged to a family where a homologous structure was known. In the following section we discuss some examples of domain atrophy that we have discovered. A full list of true and putative domain atrophy cases can be found in Additional file 1.
N-terminal end-bounded atrophy
Bacterial luciferase domain (Pfam: PF00296)
The bacterial luciferase domain of the non-fluorescent flavoprotein (NFP) from Photobacterium phosphoreum LuxF (UniProt: P12745), shows an atrophy score of 0.31, indicating a loss of nearly one-third of the domain’s structure.
The bacterial luciferase domain, homolog of bacterial luciferase subunits , is present mostly among members of gammaproteobacteria. NFP is thought to function as a ‘molecular sponge’ by sequestering myristylated flavine mononucleotide, the side-product of the bioluminescence pathway .
NFP of P. phosphoreum (PDB: 1FVP) resembles a partial TIM-barrel-like fold with one β-strand and three α-helices missing  (Figure 3A). The complete structure of the bacterial luciferase domain from Bacillus cereus (PDB: 2B81) was chosen as the reference domain for comparison. The reference domain has a (β/α)8 TIM-barrel fold (Figure 3B) with a characteristic β-barrel structure consisting of eight alternating β-strands and α-helices.
Although identified as N-terminal end-bounded atrophy, comparison with the reference domain (PDB: 2B81) shows that β1, α1 and β2 of 2B81 structurally aligns with the corresponding structural elements at the N-terminal region of 1FVP. However α2, β3, α3 and α4, of the reference domain 2B81 (residues 61-125, 132-192), have no corresponding matching structural elements on 1FVP. This suggests that P. phosphoreum bacterial luciferase domain is an example of within-domain atrophy rather than N-terminal end-bounded atrophy as first thought. We also identified a similar case of domain atrophy in the orthologous protein of Photobacterium leiognathi (UniProt: P09142, PDB: 1NFP).
A key issue is in determining whether domain atrophy has occurred is whether the ancestral fold of bacterial luciferase domains was indeed a TIM-barrel. Phylogenetic analysis of bacterial luciferase family indicates that the LuxF protein clade is completely enclosed by LuxB (Figure 4) and that the ancestral fold must be the complete TIM-barrel fold observed in the LuxB and LuxA proteins.
Does the structure shed any light on how an atrophied domain can tolerate a large deletion? The TIM-barrel fold has a central hydrophobic β-barrel surrounded by alpha helices, which shield the barrel from the solvent. However, the monomeric structure of the atrophied bacterial luciferase domain shows a barrel with a large hydrophobic cleft apparently exposed to the solvent. The atrophied domain buries this hydrophobic surface by forming a homo-dimeric protein interaction interface (Figure 3C) [16,17]. The hydrophobic core residues thus exposed by the atrophy event are stabilised by interactions at the new dimer interface. We hypothesise that such a compensatory mechanism may be a crucial factor in the selection of the otherwise unstable atrophied domain. Interestingly bacterial luciferase, in its complete form , also functions as a dimer in vivo and forms a homo-dimer with the dimer interface being on the same side of the molecule as observed in the atrophied domain-dimers (Figure 3D).
AMP-binding domain (Pfam: PF00501)
The Adenosine monophosphate (AMP) binding domain of phenylacetate-coenzymeA ligases from Burkholderia cenocepacia (UniProt: B4E7B5) (Figure 5A), exhibits structural loss at the N-terminus with an atrophy score of 0.16 corresponding to a loss of 67 amino acid residues.
Phenylacetate-coenzymeA ligases are adenylate-forming enzymes, which are involved in the metabolism of phenylacetate . The enzyme links the phosphoryl moiety of AMP to the carboxyl group of the substrate thereby activating it before transferring to the acceptor CoA .
The AMP-binding domain is present in proteins such as acyl-CoA synthetases , aryl-CoA ligases , firefly luciferases  and non-ribosomal peptide synthetases . The N-terminal AMP-binding domain of the reference structure, 4-chlorobenzoyl CoA ligase from Alcaligenes sp. (PDB: 3CW9), shows a large, globular, distorted α/β structure comprising three distinct β-sheets with a large cleft containing the binding pocket  (Figure 5B).
Structural superposition of the ‘N-terminal microdomains’ of B. cenocepacia (PDB: 2Y27)  onto the reference structure (PDB: 3CW9) reveals a marked absence of the first sub-domain (residues 1-150) (Figure 5C).
It can be seen from the reference structure that the first sub-domain has minimal interactions with the substrates during adenylation and the subsequent thioester formation processes . The homo-dimeric complex of B. cenocepacia phenylacetate-coenzymeA ligase results in the formation of an intramolecular β-sheet at the dimer interface . Interestingly this interaction, involving the second sub-domains (residues 151-322), forms an arrangement that is structurally similar to the first sub-domain  (Figure 5D). Again a homo-dimeric arrangement is used to ensure that residues exposed by domain atrophy are buried thereby compensating for domain atrophy. The atrophied domain of B. cenocepacia is thought to be functional and is involved in recruiting subsequent enzymes of the phenylacetic acid degradation pathway . In addition to atrophy observed in B. cenocepacia, we find similar domain atrophies in the phenylacetate-coenzymeA ligase orthologue from Bacteroides thetaiotaomicron (UniProt: Q8AAN6, PDB: 3QOV) and the phenazine antibiotic biosynthesis protein from Enterobacter agglomerans (UniProt: Q8GPH0, PDB: 3HGU).
C-terminal end-bounded atrophy
Ral-GTPase-activating protein domain (Pfam: PF02145)
The Ral-GTPase-activating protein (RapGAP) domain of the rat Ral-GTPase-activating subunit α-1 (UniProt: O55007) has a C-terminal atrophy score of 0.46.
The 747 amino acid residue isoform 1 has an atrophied RapGAP domain at the C-terminus with no other detectable domains at the N-terminus. The RapGAP catalytic domain of the human Rap1-GAP protein (PDB: 3BRW) was chosen as reference to infer atrophy. The atrophied rat RapGAP domain and the reference human RapGAP catalytic domain show 30% sequence identity in the aligned region.
GTPase-activating proteins (GAPs) terminate G-protein signalling by inducing hydrolysis of bound GTP to GDP . The human Rap1GAP reference protein comprises an N-terminal dimerisation and a C-terminal catalytic domain and interacts with Rap1B through the catalytic domain (Figure 6). Sequence mapping of the rat RapGAP domain onto the human Rap1-GAP protein (PDB: 3BRW) indicates structural loss in the catalytic domain (residues 301-414). The catalytic domain is an α/β structure with mixed parallel/antiparallel arrangement of β-strands and a conserved C-terminal alpha helix. The catalytic centre comprising Asn290 is close to the nucleotide-binding region and the protein interface [27,28]. The atrophy does not affect the catalytic centre or residues involved in Rap1B interaction suggesting the atrophied domain may be functional.
Rat RapGAP protein isoform 2 (906 amino acids) and isoform 3 (2,035 amino acids), however, show no signs of atrophy, suggesting exon loss mediated by alternative splicing as a probable mechanism for the observed atrophy in isoform 1.
Upstream domain-bounded atrophy
Of the 125 cases of upstream domain-bounded atrophy identified by our pipeline none were confirmed as true atrophy events.
Downstream domain-bounded atrophy
2-hydroxyacid dehydrogenase, NAD binding domain (Pfam: PF02826)
The 2-hydroxyacid dehydrogenase NAD-binding (2-Hacid_dh_C) domain of Staphylococcus aureus PurK (UniProt: A6QFS4) (Figure 7A) shows an atrophy score of 0.23.
2-Hacid_dh_C domains are found in dehydrogenases and oxidoreductases among prokaryotes and eukaryotes. Bacterial PurK and PurE are involved in the two-step conversion of 5-aminoimidazole ribonucleotide to 4-carboxy-5-aminoimidazole ribonucleotide .
The reference 2-Hacid_dh_C domain from Lactobacillus jensenii D-lactate dehydrogenase (PDB: 4PRL) adopts a Rossmann fold with six parallel β-strands bound to NAD  (Figure 7B). The S. aureus (PDB: 3ORQ) 2-Hacid_dh_C domain adopts a partial Rossmann fold and have an additional alpha helix (α2) and β-strand (β3) insertion .
Structural superposition on 4PRL indicates a partial domain at the C-terminus of 3ORQ (Figure 7C). The atrophied domain with four parallel β-strands and four α-helices superimposes well onto the N-terminus of 4PRL, with β-strands β2-β1-β4 aligning on β2-β1-β3 of the reference domain, respectively, but has no structurally equivalent residues for β4-β6 of the reference domain (residues 225-298).
It could also be seen that despite atrophy at the C-terminus, the domain retains the ‘reverse’ Rossmann fold motif (GXXGXG) on the loop, at the N-terminus, connecting β1 with α1. This binding motif and other interaction sites are conserved and are at structurally equivalent locations (3ORQ: 16-21, 38-41) compared to the NAD bound 2-Hacid_dh_C of L. jensenii (4PRL: 153-158, 175-178) and might bind to the adenine moiety, however other residues interacting with the nicotinamide moiety at the C-terminus (4PRL: Asp259 and His295) are absent. His295 located near the nicotinamide ring is thought to be catalytically important [30,31] but is absent in the atrophied domain. The 2-Hacid_dh_C domain of Bacillus anthracis (UniProt: C3PBM5, PDB: 3Q2O) also shows similar atrophy.
RNase_E_G domain (Pfam: PF10150)
The RNase E domain of Pyrococcus furiosus RNA-binding protein AU-1 (UniProt: Q8U4Q7) shows atrophy at the C-terminus with a score of 0.19.
The P. furiosus RNA-binding protein is a large oligomeric complex that binds specifically to AU-rich RNA sequences and may be involved in RNA metabolic processes . The RNA-binding protein AU-1 has a large conserved N-terminal RNase E domain, which shares 25% sequence similarity with E. coli RNase E , and a C-terminal domain of unknown function DUF402. There are no experimental structures available of the P. furiosus RNase E domain and therefore the region of atrophy was inferred through sequence mapping onto the reference structure.
The E. coli RNase E large domain (PDB: 2C0B) is a large multi-domain structure comprising four domains: S1, 5’ sensing region, RNase H and DNase I . Mapping of the P. furiosus RNase E large domain sequences on to the E. coli RNase E reference domain shows the region of atrophy corresponds to the C-terminus of the DNase I domain (residues 339-393) (Figure 8).
The DNase I domain has an α/β arrangement of two α-helices and an antiparallel β-sheet with six antiparallel β-strands in the order β1–2–3–4–6–5. The DNase I domain is the catalytic centre of the complex with active sites Asp303 and Asp346 present on β3 and β4, respectively . Asp303 and Asp346 coordinate a Magnesium ion, which, through nucleophilic attack, cleaves the scissile phosphate on the RNA backbone . Pairwise sequence alignment indicates E. coli and P. furiosus RNase E large domains share 14% sequence identity and atrophy at the C-terminus of the DNase I domain results in the loss of one of the active site residues, Asp346. However, a more accurate multiple sequence alignment of the domain family indicates that both catalytic residues may be present resulting in a potentially functional domain.
The E. coli RNase E catalytic domain is a tetrameric complex formed by a dimer-of-dimers, where the dimer interface of this complex is largely mediated by the interaction of DNase I domains , however, P. furiosus RNA-binding protein AU-1 forms a homo-oligomeric complex that is mainly trimeric . The atrophy of DNase I sub-domain and the presence of a much larger C-terminal domain DUF402, compared to the E. coli RNase E small domain, is thought to cause the difference in the oligomerisation of the complex.
Glycosyl hydrolase family 10 (Pfam: PF00331)
The Glycosyl hydrolase family 10 (Glyco_hydro_10) domain of Cryptococcus albidus Endo-1,4-beta-xylanase (UniProt: P07529) is found with a within-domain atrophy score of 0.24.
Xylanase are found in bacteria, fungi and other microbes [35,36], which degrade hemicellulose by breaking down beta-1,4-xylan into xylose. C. albidus endo-1,4-beta-xylanase is an inducible extracellular enzyme with xylobiose as its natural inducer .
The Glyco_hydro_10 domain of Thermotoga petrophila endo-1,4-beta-xylanase reference protein (PDB: 3NJ3) folds into a TIM-barrel structure . Sequence mapping onto the reference domain indicates a loss of β-strand β4 and two core alpha helices, α3 and α4 (residues 110-176) (Figure 9). Atrophy of the C. albidus Glyco_hydro_10 domain would result in the loss of one of the active sites Glu150, which is present on β4 along with Asn149, which interacts with the ligand xylobiose.
It has been shown that C. albidus endo-1,4-beta-xylanase is inefficient in degrading xylobiose and hydrolyses xylotriose at very slow rates although longer substrates such as xylotetraose are hydrolysed [37,39]. It is surprising that the loss of one of the catalytic residues Glu150 could still make the enzyme functional, albeit with a reduced efficiency. Other open TIM-barrels comprising of only seven of the eight β-strands are observed in the C-terminal domains of yeast nicotinate phosphoribosyltransferase (PDB: 1VLP) and Pyrococcus furiosus Nicotinate-nucleotide pyrophosphorylase (PDB: 2I14). The hydrophobic cores of these partial TIM-barrels are buried by interactions with their N-terminal domains. It would be very interesting to see the structure of C. albidus endo-1,4-beta-xylanase to better understand how the atrophy event is compensated for.
Proteins evolve by a number of different mechanisms that lead to changes in sequence, structure and function. Events such as domain gain or loss have been extensively studied in terms of evolution of prokaryotic and eukaryotic genomes [40,41], other minor events such as gain or loss of secondary structural elements, which cause structural variations in the domain, are less well studied. Domain loss has been observed as a major event during eukaryotic evolution and interestingly hints at an inverse trend with organismal complexity . Unlike domain loss wherein a complete domain is absent, we have investigated events of significant structural loss occurring in a domain, by analysing variation of domain boundaries across 14,831 domain families in Pfam 27.0.
Domain atrophy events were classified into five groups based on the location of the domain within the protein and the region of atrophy in the domain. Domain atrophy events are found to be rare accounting for only a small fraction (0.06%) of the total domain instances analysed. One of the reasons for such rare occurrence is that significant loss of structural elements would usually be detrimental for domain stability. A large number of factors such as hydrophobic  and electrostatic interactions  influence the native fold and stability of the domain, and unlike domain elaborations, loss of core structural elements involved in such interactions would not be favourable.
An intriguing aspect of atrophied domains is the energetics of their folds. Deletions or insertions of individual residues can significantly influence the overall energetics of the fold . Deletion of structural elements, as observed in atrophied domains, would contribute to significant changes in the energetics compared to the native fold and as such warrant major compensatory mechanisms to stabilise fold. It was shown that hydrophobic contacts at the domain-domain interface stabilise the otherwise unstable domains in a multi-domain protein . Formation of multimeric complexes or domain-domain interactions may serve as one of the compensatory mechanisms for stabilising atrophied domains as observed in the cases of P. leiognathi bacterial luciferase domain and the AMP-binding domain of B. cenocepacia.
Atrophied domains should have incomplete/partial folds, but the converse may not be true. For example, the histidine biosynthesis pathway enzymes HisA and HisF consist of two copies of a (β/α)4 half-barrel . It has been proposed that the complete (β/α)8 TIM-barrel may have evolved from duplication and gene fusion of an ancestral (β/α)4 half-barrel .
An important issue is to understand how domain atrophy events arise. Mutational events such as a mutation to a premature stop codon could lead to domains with atrophy at the C-terminus. While the creation of a new downstream start codon could lead to domains with atrophy at the N-terminus. Other mechanisms such as exon loss or deletion could additionally lead to atrophy within the domain. We observe a large number of domain atrophy cases at the termini of proteins: N-termini (n = 40) and C-termini (n = 26), compared to atrophy cases within a protein. We see just four cases of atrophies within a domain, and we only see five cases of atrophy events that are bounded by other domains, and hence also occur within the protein. It is unclear whether this is a reflection of the true proportions or simply due to our methodology. Alternative splicing provides one possible mechanism that enables an organism to explore atrophied domain sequence space while still retaining a functional copy . Analysis of splice variants from the ENCODE dataset indicates a large number of splice sites occur within functional domains , suggesting a plausible mechanism for domain atrophy. It was previously shown that alternative splicing events in human result in fewer partial domains than expected . These partial domains were reported to affect a significant number of functional sites  and may translate into proteins with altered structures and functions [47,48]. Strong selection pressure may also play a role in the evolution of atrophied domains, selecting only atrophied domains with an important structural or functional role for the organism that must be kept.
Unfolded proteins or misfolded proteins are targeted by protein degradation machinery by recognising their solvent exposed hydrophobic residues. Although no clear experimental data are available from known cases of atrophied domains on their possible mechanisms to evade recognition and degradation by the protein degradation machinery, we speculate atrophied domains escape unfolding and degradation by burying their exposed solvent-accessible areas through dimerisation or complex formations.
Disordered or transiently unfolded regions are known to fold cooperatively upon contact with their binding partners  and evade degradation until then. Chaperones, such as BiP, maintain unfolded or misfolded proteins in a folding-competent state . We think atrophied domains with solvent exposed hydrophobic residues may be similarly protected by chaperones until folded stably upon complex formation.
We observe that domain atrophy is a very rare event compared to domain loss or elaborations. Our analyses of a subset of protein sequences, with protein existence value of 1, defined by UniProt, indicates only a minute fraction of domains exhibit atrophy. An independent analysis  of a broader sequence space also corroborates the findings. A large percentage of the protein domains in Pfam, annotated as partial domains, are largely found to be profile-HMM alignment or sequence annotation artefacts . Our work has greatly increased the number of known atrophy cases. As identified by Triant and Pearson it is very difficult to identify true cases of domain atrophy from artefacts. Availability of high quality data with refined protein domain models and sequences, together with experimental structures of such atrophied domains, can further shed light on domain atrophy and its role in evolution of protein domains. Our work presents potential domain atrophy cases for further experimentation that will lead to an improved understanding of how proteins tolerate extreme mutations.
To identify potential cases of domain atrophy we use matches of the UniProt sequence database (release 2012_06) against the profile-HMM models from the Pfam database release 27.0 . This set of matches contains 28,738,352 Pfam-A protein domain instances across 14,831 families in 18,523,877 protein sequences.
To determine which domains may be cases of atrophy we calculate a measure called the Atrophy Score (AS) at both the N- (ASN) and C-terminal (ASC) boundaries of each domain instance. Figure 10 shows a schematic representation of the parameters used in calculating domain atrophy score. The equations used to calculate the atrophy score are shown below:
Where, ASN is the atrophy score at the N-terminus of the domain, DN is the number of unmatched HMM match states at the N-terminus of the domain, dN is the inter-domain distance or domain interval that is the number of amino acid residues between the domain and its adjacent upstream domain or the sequence start site in the case of an N-terminal domain, and L is the HMM model length of the domain family. Similarly ASC, DC and dC correspond to atrophy score, the number of unmatched HMM match states and the inter-domain distance of the C-terminus of the domain, respectively.
Instances of within-domain atrophy can be identified in cases where the profile-HMM matches to a single domain have been split into two profile-HMM matches, with the first corresponding to the N-terminal part of the domain and the second corresponding to the C-terminal part of the domain. We can distinguish these cases from tandem repeats by considering the HMM match states of each domain. The start HMM-match state of the downstream domain is greater than the end HMM-match state of the upstream domain in cases of split domains. The computation of within-domain atrophy score (ASW) is similar to ASN.
Where DW is the number of unmatched HMM match states within the domain, dW is the domain interval within the domain and L is the HMM model length of the domain family.
Alignment co-ordinates of each domain are considered to calculate the inter-domain interval. An intuitive description of the Atrophy Score would be that a score of 0.33 means that one-third of the length of the domain has been lost to domain atrophy.
We classified domain atrophy events into five types, based on the domain location (architecture) in the protein and the region of atrophy in the domain: (1) N-terminal end-bounded atrophy: structural loss at the N-terminal region of the N-terminal domain; (2) C-terminal end-bounded atrophy: structural loss at the C-terminal region of the C-terminal domain; (3) Upstream domain-bounded atrophy: structural loss at the N-terminal region of an inner domain, also including the N-terminal region of the C-terminal domain; (4) Downstream domain-bounded atrophy: structural loss at the C-terminal region of an inner domain, also including the C-terminal region of the N-terminal domain; and (5) Within-domain atrophy: structural loss within the domain.
Initial results from applying the atrophy score to all UniProt proteins showed that there were numerous common failure modes (Figure 11) that would mask our ability to find genuine domain atrophy events. Therefore, we applied a set of filters to reduce the number of false positive matches with high atrophy scores.
Of a total of 23,193,494 sequences in the database, 18,523,877 sequences, which have at least one Pfam-A domain instance, were used in the analysis. Initial filtering was applied to exclude domain models from sequences with protein existence (PE) levels of 2 to 5. These are enriched in gene prediction errors and fragment proteins. This reduced the number of sequence considered from 18,523,877 to 77,305. Proteins with a protein existence level of 1 have clear experimental evidence for the existence of the protein from Edman sequencing, mass spectrometry, X-ray, NMR or other experimental evidence. Although not strictly a measure of protein sequence quality these proteins usually have highly accurate protein sequences. We also removed sequences annotated as fragments in UniProt, which further reduced the set of sequences considered from 77,305 to 75,435.
Adjacent domains that are of the same clan, similar to figure 11B, could lead to ambiguous domain boundary assignments at the interval and hence such cases were filtered out to avoid detecting false atrophy events. The resulting final-set comprising 114,303 domain instances from 75,435 sequences were included in the analysis.
Perl scripts were run to calculate atrophy scores across the set of domains. Domain instances with an atrophy score of 0.15 or more were further investigated.
Domain instances that were obtained after applying the above filters were then selected for manual inspection. Only those domains that had an atrophy score of 0.15 or above were checked manually for identification of false positives. We checked each potential domain atrophy case for evidence of any of the following failure modes (see Figure 11): (1) Gene prediction errors: We checked whether the missing part of a domain was to be found in an adjacent gene or due to an incomplete gene prediction; (2) Nested domains: We checked whether a high atrophy score was due to a domain nesting within another. These were considered as false positives; (3) Multi-domain families: Due to incorrect Pfam domain definitions some Pfam domains actually correspond to multiple structural domains that can be found independently. We checked the structure of each Pfam family to confirm whether this was the cause of a high atrophy score; (4) Small domains: Domains of length less than 30 amino acid residues were not considered since atrophy score of 0.15 and above of small domains correspond to loss of a single secondary structural element or a part thereof, which is not considered true atrophy; (5) Circular permutations: While circular permuted domains are complete domains, the rearrangement of domain HMM start-site and HMM end-site with respect to their domain HMM-model would result in misidentification of such cases as domain atrophies; (6) Short repeats: Domains composed of tandem structural motifs, such as β-propeller, β- or α-helix, are made of short repeating sequence motifs were considered as false positives. Addition or removal of repeats is often tolerated in terms of protein mutation; and (7) Disordered domains: Inferring domain atrophy among intrinsically disordered protein domains is not straightforward mainly owing to their lack of native ordered tertiary structure and such cases were considered false positives.
Apart from the above failure modes domain atrophy cases whose structures were theoretically modelled or had no other reference structures in the family to compare with were also treated as failure modes. Other cases of complete structural domains but identified as domain atrophy were treated as false positives. Table 2 lists the occurrences of different types of failure modes.
From a total of 1,636 domain instances, with atrophy scores between 0.15 and 1, which were manually checked, 1,287 domain instances were classified as failure modes or false positives. The positive predictive value (PPV) of our method is 0.055.
Visualisation of structures
Structures were visualised with Chimera . Where experimental structures of atrophied domains were not available, the shortest full-length domain structure within the domain family was chosen as the reference. The extent of domain atrophy was then inferred by a pairwise sequence alignment guided mapping of unaligned sequence regions onto the reference structure. Instances of putative domain atrophy where no full-length reference structure was available for the family were not considered further.
Evolutionary information was inferred from phylogenetic trees constructed from multiple sequence alignment of domain family seed sequences and homologous sequences from a JackHMMER search . A non-redundant set of sequences of 90% identity or less was aligned with Mafft  with default settings. Alignments were visualised with Belvu  and phylogenetic trees constructed using the neighbour-joining method present in Belvu using default parameters.
Weiner J 3rd, Beaussart F, Bornberg-Bauer E. Domain deletions and substitutions in the modular protein evolution. FEBS J. 2006;273:2037–47.
Buljan M, Frankish A, Bateman A. Quantifying the mechanisms of domain gain in animal proteins. Genome Biol. 2010;11:R74.
Kim R, Guo JT. Systematic analysis of short internal indels and their impact on protein folding. BMC Struct Biol. 2010;10:24.
Hormozdiari F, Salari R, Hsing M, Schonhuth A, Chan SK, Sahinalp SC, et al. The effect of insertions and deletions on wirings in protein-protein interaction networks: a large-scale study. J Comput Biol. 2009;16:159–67.
Sandhya S, Rani SS, Pankaj B, Govind MK, Offmann B, Srinivasan N, et al. Length variations amongst protein domain superfamilies and consequences on structure and function. PLoS One. 2009;4:e4981.
Pascarella S, Argos P. Analysis of insertions/deletions in protein structures. J Mol Biol. 1992;224:461–71.
Taylor MS, Ponting CP, Copley RR. Occurrence and consequences of coding sequence insertions and deletions in Mammalian genomes. Genome Res. 2004;14:555–66.
Dessailly BH, Redfern OC, Cuff AL, Orengo CA. Detailed analysis of function divergence in a large and diverse domain superfamily: toward a refined protocol of function classification. Structure. 2010;18:1522–35.
Nardini M, Pesce A, Milani M, Bolognesi M. Protein fold and structure in the truncated (2/2) globin family. Gene. 2007;398:2–11.
Grishin NV. Fold change in evolution of protein structures. J Struct Biol. 2001;134:167–85.
Das D, Murzin AG, Rawlings ND, Finn RD, Coggill P, Bateman A, et al. Structure and computational analysis of a novel protein with metallopeptidase-like and circularly permuted winged-helix-turn-helix domains reveals a possible role in modified polysaccharide biosynthesis. BMC Bioinformatics. 2014;15:75.
Fowler SB, Best RB, Toca Herrera JL, Rutherford TJ, Steward A, Paci E, et al. Mechanical unfolding of a titin Ig domain: structure of unfolding intermediate revealed by combining AFM, molecular dynamics simulations, NMR and protein engineering. J Mol Biol. 2002;322:841–9.
Triant DA, Pearson WR. Most partial domains in proteins are alignment and annotation artifacts. Genome Biol. 2015.
Moore SA, James MN. Structural refinement of the non-fluorescent flavoprotein from Photobacterium leiognathi at 1.60 A resolution. J Mol Biol. 1995;249:195–214.
Kita A, Kasai S, Miyata M, Miki K. Structure of flavoprotein FP390 from a luminescent bacterium Photobacterium phosphoreum refined at 2.7 A resolution. Acta Crystallogr D Biol Crystallogr. 1996;52:77–86.
Moore SA, James MN, O’Kane DJ, Lee J. Crystal structure of a flavoprotein related to the subunits of bacterial luciferase. EMBO J. 1993;12:1767–74.
Moore SA, James MN. Common structural features of the luxF protein and the subunits of bacterial luciferase: evidence for a (beta alpha)8 fold in luciferase. Protein Sci. 1994;3:1914–26.
Waddle JJ, Johnston TC, Baldwin TO. Polypeptide folding and dimerization in bacterial luciferase occur by a concerted mechanism in vivo. Biochemistry. 1987;26:4917–21.
Martinez-Blanco H, Reglero A, Rodriguez-Aparicio LB, Luengo JM. Purification and biochemical characterization of phenylacetyl-CoA ligase from Pseudomonas putida. A specific enzyme for the catabolism of phenylacetic acid. J Biol Chem. 1990;265:7084–90.
Shah MB, Ingram-Smith C, Cooper LL, Qu J, Meng Y, Smith KS, et al. The 2.1 A crystal structure of an acyl-CoA synthetase from Methanosarcina acetivorans reveals an alternate acyl-binding pocket for small branched acyl substrates. Proteins. 2009;77:685–98.
Gulick AM, Lu X, Dunaway-Mariano D. Crystal structure of 4-chlorobenzoate:CoA ligase/synthetase in the unliganded and aryl substrate-bound states. Biochemistry. 2004;43:8670–9.
Conti E, Franks NP, Brick P. Crystal structure of firefly luciferase throws light on a superfamily of adenylate-forming enzymes. Structure. 1996;4:287–98.
Lee TV, Johnson LJ, Johnson RD, Koulman A, Lane GA, Lott JS, et al. Structure of a eukaryotic nonribosomal peptide synthetase adenylation domain that activates a large hydroxamate amino acid in siderophore biosynthesis. J Biol Chem. 2010;285:2415–27.
Reger AS, Wu R, Dunaway-Mariano D, Gulick AM. Structural characterization of a 140 degrees domain movement in the two-step reaction catalyzed by 4-chlorobenzoate:CoA ligase. Biochemistry. 2008;47:8016–25.
Law A, Boulanger MJ. Defining a structural and kinetic rationale for paralogous copies of phenylacetate-CoA ligases from the cystic fibrosis pathogen Burkholderia cenocepacia J2315. J Biol Chem. 2011;286:15577–85.
Bos JL, Rehmann H, Wittinghofer A. GEFs and GAPs: critical elements in the control of small G proteins. Cell. 2007;129:865–77.
Daumke O, Weyand M, Chakrabarti PP, Vetter IR, Wittinghofer A. The GTPase-activating protein Rap1GAP uses a catalytic asparagine. Nature. 2004;429:197–201.
Scrima A, Thomas C, Deaconescu D, Wittinghofer A. The Rap-RapGAP complex: GTP hydrolysis without catalytic glutamine and arginine residues. EMBO J. 2008;27:1145–53.
Brugarolas P, Duguid EM, Zhang W, Poor CB, He C. Structural and biochemical characterization of N5-carboxyaminoimidazole ribonucleotide synthetase and N5-carboxyaminoimidazole ribonucleotide mutase from Staphylococcus aureus. Acta Crystallogr D Biol Crystallogr. 2011;67:707–15.
Kim S, Gu SA, Kim YH, Kim KJ. Crystal structure and thermodynamic properties of d-lactate dehydrogenase from Lactobacillus jensenii. Int J Biol Macromol. 2014;68:151–7.
Tishkov VI, Matorin AD, Rojkova AM, Fedorchuk VV, Savitsky PA, Dementieva LA, et al. Site-directed mutagenesis of the formate dehydrogenase active centre: role of the His332-Gln313 pair in enzyme catalysis. FEBS Lett. 1996;390:104–8.
Kanai A, Oida H, Matsuura N, Doi H. Expression cloning and characterization of a novel gene that encodes the RNA-binding protein FAU-1 from Pyrococcus furiosus. Biochem J. 2003;372:253–61.
Callaghan AJ, Marcaida MJ, Stead JA, McDowall KJ, Scott WG, Luisi BF. Structure of Escherichia coli RNase E catalytic domain and implications for RNA turnover. Nature. 2005;437:1187–91.
Koslover DJ, Callaghan AJ, Marcaida MJ, Garman EF, Martick M, Scott WG, et al. The crystal structure of the Escherichia coli RNase E apoprotein and a mechanism for RNA degradation. Structure. 2008;16:1238–44.
Beg QK, Kapoor M, Mahajan L, Hoondal GS. Microbial xylanases and their industrial applications: a review. Appl Microbiol Biotechnol. 2001;56:326–38.
Polizeli ML, Rizzatti AC, Monti R, Terenzi HF, Jorge JA, Amorim DS. Xylanases from fungi: properties and industrial applications. Appl Microbiol Biotechnol. 2005;67:577–91.
Biely P, Kratky Z, Vrsanska M, Urmanicova D. Induction and inducers of endo-1,4-beta-xylanase in the yeast Cryptococcus albidus. Eur J Biochem. 1980;108:323–9.
Santos CR, Meza AN, Hoffmam ZB, Silva JC, Alvarez TM, Ruller R, et al. Thermal-induced conformational changes in the product release area drive the enzymatic activity of xylanases 10B: Crystal structure, conformational stability and functional characterization of the xylanase 10B from Thermotoga petrophila RKU-1. Biochem Biophys Res Commun. 2010;403:214–9.
Biely P, Kratky Z, Vrsanska M. Substrate-binding site of endo-1,4-beta-xylanase of the yeast Cryptococcus albidus. Eur J Biochem. 1981;119:559–64.
Zmasek CM, Godzik A. Strong functional patterns in the evolution of eukaryotic genomes revealed by the reconstruction of ancestral protein domain repertoires. Genome Biol. 2011;12:R4.
Nasir A, Kim KM, Caetano-Anolles G. Global patterns of protein domain gain and loss in superkingdoms. PLoS Comput Biol. 2014;10:e1003452.
Pace CN, Fu H, Fryar KL, Landua J, Trevino SR, Shirley BA, et al. Contribution of hydrophobic interactions to protein stability. J Mol Biol. 2011;408:514–28.
Kumar S, Nussinov R. Close-range electrostatic interactions in proteins. Chembiochem. 2002;3:604–17.
Randles LG, Lappalainen I, Fowler SB, Moore B, Hamill SJ, Clarke J. Using model proteins to quantify the effects of pathogenic mutations in Ig-like proteins. J Biol Chem. 2006;281:24216–26.
Bhaskara RM, Srinivasan N. Stability of domain structures in multi-domain proteins. Sci Rep. 2011;1:40.
Lang D, Thoma R, Henn-Sax M, Sterner R, Wilmanns M. Structural evidence for evolution of the beta/alpha barrel scaffold by gene duplication and fusion. Science. 2000;289:1546–50.
Birzele F, Csaba G, Zimmer R. Alternative splicing and protein structure evolution. Nucleic Acids Res. 2008;36:550–8.
Tress ML, Martelli PL, Frankish A, Reeves GA, Wesselink JJ, Yeats C, et al. The implications of alternative splicing in the ENCODE protein complement. Proc Natl Acad Sci U S A. 2007;104:5495–500.
Kriventseva EV, Koch I, Apweiler R, Vingron M, Bork P, Gelfand MS, et al. Increase of functional diversity by alternative splicing. Trends Genet. 2003;19:124–8.
Gruszka DT, Wojdyla JA, Bingham RJ, Turkenburg JP, Manfield IW, Steward A, et al. Staphylococcal biofilm-forming protein has a contiguous rod-like structure. Proc Natl Acad Sci U S A. 2012;109:E1011–8.
Schroder M, Kaufman RJ. The mammalian unfolded protein response. Annu Rev Biochem. 2005;74:739–89.
Finn RD, Bateman A, Clements J, Coggill P, Eberhardt RY, Eddy SR, et al. Pfam: the protein families database. Nucleic Acids Res. 2014;42:D222–30.
Pettersen EF, Goddard TD, Huang CC, Couch GS, Greenblatt DM, Meng EC, et al. UCSF Chimera–a visualization system for exploratory research and analysis. J Comput Chem. 2004;25:1605–12.
Finn RD, Clements J, Eddy SR. HMMER web server: interactive sequence similarity searching. Nucleic Acids Res. 2011;39:W29–37.
Katoh K, Standley DM. MAFFT multiple sequence alignment software version 7: improvements in performance and usability. Mol Biol Evol. 2013;30:772–80.
Sonnhammer EL, Hollich V. Scoredist: a simple and robust protein sequence distance estimator. BMC Bioinformatics. 2005;6:108.
We gratefully acknowledge EMBL for funding this work. We thank Bill Pearson and Deborah Triant for comments on our manuscript.
The authors declare that they have no competing interests.
AB and AP conceived and designed the study. AP wrote the scripts and analysed the data. AP and AB wrote the manuscript. Both authors read and approved the final manuscript.
About this article
Cite this article
Prakash, A., Bateman, A. Domain atrophy creates rare cases of functional partial protein domains. Genome Biol 16, 88 (2015) doi:10.1186/s13059-015-0655-8
- Reference Domain
- Structural Loss
- Domain Loss
- Bacterial Luciferase
- Partial Domain