article: Many proteins fold into highly regular and repetitive three dimensional structures. The analysis of structural patterns and repeated elements is fundamental to understand protein function and evolution. We present recent improvements to the CE-Symm tool for systematically detecting and analyzing the internal symmetry and structural repeats in proteins. In addition to the accurate detection of internal symmetry, the tool is now capable of i) reporting the type of symmetry, ii) identifying the smallest repeating unit, iii) describing the arrangement of repeats with transformation operations and symmetry axes, and iv) comparing the similarity of all the internal repeats at the residue level. CE-Symm 2. 0 helps the user investigate proteins with a robust and intuitive sequence-to-structure analysis, with many applications in protein classification, functional annotation and evolutionary studies. We describe the algorithmic extensions of the method and demonstrate its applications to the study of interesting cases of protein evolution.
  Several definitions of internal symmetry and repeats are possible, depending on the biological question of interest. For the purposes of this paper, we define it as the regular arrangement of a common repeating structural unit within a protein chain. Therefore, a repeat is an asymmetric structural motif present multiple times in the same structure. We restrict our consideration of repeats to cases where the orientation between adjacent structural units is regular; that is, where a consistent geometric transformation can be applied to superimpose each repeat onto the next. In other words, CE-Symm focuses on identifying repeats which conserve both the structure and the interface between repeats. Several types of internal symmetry can be derived from this broad definition. The most basic division is between closed symmetry and open symmetry. In proteins with closed symmetry, the repeats are arranged in a point group symmetry. This can be defined mathematically as a set of rotations that superimpose equivalent repeats while keeping at least one point at the center of rotation fixed. In contrast, repeats in proteins with open symmetry are related by transformations with a translational component. Examples of closed and open symmetry can be found in Fig 1a–1d and 1e–1h, respectively. Closed symmetries can be further characterized according to the possible chiral point groups: cyclic (Cn), generated by a single n-fold rotational operator (Fig 1a–1b); dihedral (Dn), which requires an n-fold rotation and n perpendicular 2-fold operators (Fig 1c–1d); and polyhedral point-groups (T, O, and I), which feature non-perpendicular rotation operators. Both cyclic and dihedral internal symmetries are common in proteins, but, although common at the quaternary structure level, polyhedral symmetries have not yet been observed within a single polypeptide chain. Open symmetry can be further subdivided into special cases of helical, translational, and superhelical repeats. Helical symmetry consists of repeats arranged around a screw axis, where each repeat is related to the next by a fixed linear translation combined by a rotation around the central axis (Fig 1e). In cases where the rotation angle is close to an fraction of a turn, we indicate the approximate number of subunits needed per turn (Hn). Proteins with open symmetry that have negligible translation are called rotational repeats (Fig 1f), and those with negligible rotation between repeats are called translational repeats (Fig 1g), both annotated as R. Superhelical symmetry (SH) provides the most general description of repeats with open symmetry, and is reserved for cases which cannot be expressed as a single fixed operator relating each repeat to the next. Instead, the rotation axis between adjacent repeats precesses along a helical path (Fig 1h). Proteins with open symmetry are sometimes referred to as solenoid proteins [28].
  CE-Symm begins with a structural self-alignment (other than the identity alignment) of the input protein structure using the Combinatorial Extension (CE) algorithm [29] (Fig 2b and 2c). Identifying significant self-alignments was the primary focus of the first version of the algorithm [21]. In the self-alignment of structures with closed symmetry the first and last repeats are aligned, forming a circular permutation (CP) of the structure. This is why the structure alignment method used in CE-Symm shares algorithmic primitives with CE-CP [30]. For proteins with open symmetry, the initial self-alignment will always be missing one of the repeats due to the translation component of the symmetry operator. The alignment quality is quantified using TM-Score [31]. Both irregularly arranged repeats and large asymmetric regions in a structure will reduce the score of the self-alignment. In addition, open symmetry will generally have lower scores than closed symmetry, because the terminal repeats are unaligned in the initial self-alignment. The order of symmetry is defined as the number of symmetric units (repeats) in a structure. Extracting the order of symmetry is a key part of symmetry detection, and subsequent steps of the CE-Symm method depend on its correctness. Two methods to automatically determine the order of symmetry in closed structures were described in the previous CE-Symm publication [21]: DeltaPosition and RotationAngle. An error in the distance formula was corrected in CE-Symm 2. 0 (see S1 Text, Supplemental Methods), but the DeltaPosition method still gives better overall performance and is used by default for closed symmetry. A third method for order detection which is able to handle open symmetry has been introduced in the new version of the tool and is named GraphComponent. Conceptually, the self-alignment is treated as a directed graph over the set of aligned residues (Fig 2d). Residues that are aligned in all k repeats will form a path with k nodes. For open symmetry these paths tend to be disjoint, so simply finding the most frequent size of the connected components in the graph can accurately determine the order for open symmetry. For well-aligned cases of closed symmetry, the aligned residues form a cycle of k nodes, so the same method can also work in the general case. Those residues which participate in a path or cycle of the most frequent size form the refined alignment discussed in the following section. For cases of closed symmetry, small alignment errors can lead to a situation where paths of k residues do not form a closed cycle, but rather lead to a different residue at a small offset in the sequence. This can lead to failures of the GraphComponent order detector due to the merging of multiple alignment paths. This case can be handled by the DeltaPosition order detector. Whether a protein is open or closed can be easily determined by looking for a circular permutation in the self-alignment. CE-Symm uses DeltaPosition for closed cases where a permutation is found and GraphComponent for open cases. This can also be overridden by the user when running CE-Symm. The refinement procedure takes as input the self-alignment of the structure and the order of symmetry (k) and returns a multiple alignment of the repeats (Fig 2e). CE-Symm has two implementations of the refinement procedure: GraphComponent and DeltaPosition, which are closely related to their respective order detectors. The GraphComponent refiner combines all connected components of the self-alignment graph with size equal to the order of symmetry. Each connected component contributes one column to the refined alignment of repeats. Care must be taken that the repeat sequences preserve the sequence order of the polypeptide chain; where some pairs of connected components would violate this property, some are discarded in a way that maximizes the total length of the resulting alignment. The DeltaPosition refiner takes the self-alignment graph and modifies it until all remaining nodes are part of k-cycles. The modification heuristic is described in detail in the Supplemental Methods (S1 Text). These cycles each contribute one column of the multiple alignment of the symmetric repeats, as in the GraphComponent refiner. Note that, like the GraphComponent refiner, the multiple alignments obtained at the end of this stage consist of ungapped columns, so all repeats are of the same size. The multiple alignment obtained from the refinement is sometimes far from optimal, and depends very much on the quality of the self-alignment. In addition, the refinement process prioritizes precision over coverage, which means that only the best residue equivalencies will be included, resulting in a shorter multiple alignment. The goal of the optimization is to increase the multiple alignment length while keeping the RMSD low (Fig 2f). Furthermore, the optimization procedure can improve parts of the alignment that were not fully represented in the self-alignment, and thus not captured in the refinement result. The optimization process uses a similar approach to the Combinatorial Extension Monte Carlo (CEMC) multiple structure alignment algorithm [32]. The multiple alignment can be described by a matrix, where the rows represent aligned structures and the columns represent aligned positions (residue equivalencies). Rearranging and modifying the entries of the matrix results in changes of the multiple alignment. There are four possible moves (changes in the multiple alignment): The insertion of gaps allows for partial repeat similarities in the alignment. All moves take into consideration that rows of the alignment occur sequentially in the protein sequence, so unaligned residues between repeats can be considered either at the end of a repeat or the beginning of the following one. In addition, the shrink and insert gap moves have been biased, so that the probability of choosing an alignment column or an equivalent residue, respectively, is proportional to the average inter-residue distance of the given column or the given residue, respectively. A geometric distribution with parameter 0. 5 is chosen to allocate the probability among alignment columns. A schematic representation of the steps and how they affect the multiple alignment is provided in S1 Fig. After each optimization step, an alignment score is calculated. The score function to be optimized has also been smoothed with respect to the original CEMC score to remove discontinuities: S = ∑ i = 0 N ∑ j = 0 L[ C1+ (dijd0) 2−A ]−G (1) N is the number of structures (rows) in the alignment; L is the number of equivalent positions (columns) in the alignment, including gaps; C is the maximum score of an alignment position (by default set to 20); dij is the average distance from aligned residue j in structure i to all its equivalent residues; d0 is the structural similarity function parameter, as defined by the TM-score [31]; A is the distance cutoff penalization, which shifts the function to negative values when the maximum allowed average distance of an aligned position (dc) is higher than dij; and G is a linear gap penalty term. Calculation of A using a distance cutoff parameter dc (by default set to 7Å) is straightforward from the condition that the score S has to be 0 when dij = dc. The shape of the score function for different values of dc is shown in S2 Fig. Moves with positive score changes are always accepted. The acceptance probability of a negative scoring move is proportional to the score difference and decreases proportional to the number of optimization steps as follows: p = [ C - Δ S C m ] (1 - m M) (2) m is the current iteration number, ΔS is the change in alignment score and M is the maximum number of iterations. The maximum number of optimization iterations is proportional to the length of the protein, by default a hundred times the number of residues in the protein. Optimization finishes either because it reaches the maximum number of iterations or in the case that no moves are accepted for a fraction of the total number of iterations (by default M divided by 50). So far, the procedure described can only identify symmetry operations that require a single axis. However, some structures present symmetries represented with more than one axis. This is the case for point groups other than cyclic, like dihedral symmetry, or structures with more than one level of symmetry, what we define as hierarchical symmetries. Multiple CE-Symm iterations are run in a recursive manner, i. e. repeats found in previous rounds are recursively fed into the next run until a non-significant result (no symmetry) is found. The goal is to find all the significant symmetry levels of a structure. At the end of an iteration, repeats are extracted from the internal symmetry result and one of them is chosen as the representative, by default the N-terminal repeat. Results of successive iterations are merged by combining the symmetry axes and multiple alignments, generating a unique result for the query structure. The recursive symmetry detection allows better order determination for difficult cases (e. g. , TIM barrels), because usually fractions of the order of symmetry are initially found (e. g. , 2-fold instead of 8-fold). Continuing the analysis recursively breaks the structure down to the true asymmetric repeating units (e. g. , with three levels of symmetry: 2-fold, 4-fold and finally 8-fold). There are three decision checkpoints in the algorithm flowchart in Fig 2. The first significance criterion for a symmetry result is the self-alignment TM-score. Like in the previous version, the default threshold value is set to 0. 4. The second significance criterion is the order of symmetry. A symmetric structure must have symmetry order greater than one and the refinement of the self-alignment into a multiple repeat alignment has to be successful. The third significance criterion is the average TM-score of the multiple alignment of repeats, defined as the average TM-score of all pairwise repeat alignments. The default threshold value for the average TM-score is set to 0. 36, because a 10% decrease from the original TM-score is allowed after refinement due to the restrictive conditions imposed on it. In addition the number secondary structure elements (SSE) of the final asymmetric repeating unit is considered. If the the number of SSE of each repeat is lower than the threshold, the result will be considered non-significant. For many applications it may be desirable to exclude simple repeat units (e. g. helical bundle proteins), but these are included in CE-Symm analysis by default in order to find the highest possible symmetry in a structure. The recursive symmetry detection identifies a collection of symmetry axes that describe the arrangement of repeats in the query structure. In many cases, several of these axes can be combined to form higher-order symmetries. For example, a two-fold rotation axis can be combined with another orthogonal axis to form dihedral symmetry. Near-identical rotation axes can also be combined to form higher-order rotational symmetry. To determine the point group symmetry, we build on the algorithm described by Levy et al. [33]. The symmetry axes can be found efficiently by first considering only the centroids of each repeat, since they must be in a symmetric configuration if the entire complex is symmetric. To find all possible symmetry axes, the centroids are rotated around axes that go through the centroid of the whole structure using an orientation grid search in quaternion space [34]. For each orientation, the RMSD of the aligned centroids is calculated. If the centroids align within a threshold, then all Cα atoms are superimposed. The symmetry axis is then defined by the rotation matrix of this superposition. If the RMSD is less than a threshold value (i. e. , 5 Å), the symmetry operation is considered valid. Since symmetry operations form a group, only a few are needed to complete the full point group. This procedure allows the combination of axes that have been considered separately by CE-Symm. The point group is included in the final symmetry output and displayed to the user as a polyhedron box around the protein structure. For evaluation purposes we used the manually curated dataset of 1,007 domains selected randomly from the set of SCOP superfamilies, introduced in our previous study [21]. The benchmark is intended to be representative of structural domains in the PDB, and repeats are present in 25. 8% of domains. Domains were curated to require reasonably high coverage by repeats with conserved topology, but allowing for structural divergence and flexibility. A small number of classifications were updated to be more consistent with the new symmetry definitions, especially for the cases of open symmetry. The updated version of the internal symmetry dataset (v2. 0), together with the reasons of the modified annotations, is summarized in S1 Table. The benchmark dataset and results are available in S2 Table or https: //github. com/rcsb/symmetry-benchmark. An important note is that in the evaluation of the previous version the open symmetry cases in the benchmark were part of the asymmetric (negative) set, while they are part of the symmetric (positive) set in this evaluation. RepeatsDB was used as an additional benchmark of positive cases [25]. At time of download, RepeatsDB contained 3,689 manually reviewed entries (accessed October 18,2018). Entries consist of single chains, and may contain multiple structural domains or multiple repeat regions. We selected a total of 3,503 chains with repeats of classes III (solenoid repeats) and IV (closed repeats) that were part of the RepeatsDB-lite benchmarking dataset [24]. Chains which were either annotated in RepeatsDB as having multiple repeat regions or in ECOD [35] as having multiple structural domains were considered multidomain chains in the analysis. The RepeatsDB benchmarking dataset and CE-Symm results are listed in S3 Table.
  In our previous article we compared the performance of CE-Symm against SymD. The performance in symmetry detection has only been affected by the additional order detection and alignment refinement steps. The ROC curves of both versions are very similar, with a slight reduction of false positives in the new one (S3 Fig). At the default TM-score threshold values for result significance, the false positive (FP) rate has decreased from 5. 5% to 2. 5%, while the true positive (TP) rate has been reduced from 81% to 76% on the benchmarking dataset. The bottleneck in symmetry detection continues to be finding a significant self-alignment. The different methods for order detection perform similarly for closed symmetry cases in the benchmark (S4 Table). The simpler GraphComponent method performs worse than the others, but it is the only one that can be used for open symmetries, while the DeltaPosition detector performs better than the RotationAngle method, particularly for difficult cases. On average for symmetric entries in the benchmark where CE-Symm could find symmetry, the optimization step extended the repeat length by 43%, reduced the RMSD by 1. 8% and increased the average TM-Score of the repeat alignment by 19. 6%. Furthermore, using optimization an additional 23 cases (9% of the symmetric structures in the benchmark) were correctly identified as symmetric (209 with optimization, 186 without), which is a 12% improvement in symmetry detection. Because the highest scoring alignment of the simulation trajectory is taken as the result, optimization can only improve the initial alignment. A detailed evaluation of predicted number of repeats by CE-Symm is shown in S3 Fig. Among incorrect predictions, CE-Symm tends to underpredict the number of repeats, typically as a fraction of the true number of repeats (e. g. predicting 4 repeats for a C8 TIM-barrel). Examples with high-order rotational symmetry are also often missed due to the default maximum order of 8 used in the DeltaPosition method. Overall, CE-Symm 2. 0 is able to predict the number of repeats correctly in 89% of cases. High order open symmetries remain challenging due to the prevalence of kinks and structural inhomogeneities in large open structures. The RepeatsDB-lite method was also run on the benchmark for comparison (S5 Fig). RepeatsDB-lite only detects proteins with three or more repeats, so one- and two-repeat cases in the benchmark were binned together for analysis. With this definition, it predicts the correct number of structures for 72% of the benchmark. Since the method is based on a library of known repeat units, it tends to miss repeats in benchmark cases that are not similar to the training dataset. Additionally, the method does not enforce high coverage or symmetric orientations between repeats. This is desirable when identifying candidates that may have repeats, but it leads to a rather high false positive rate (24%) relative to the definitions of repeats used when curating the benchmark. Raw data for the benchmark results is available in S3 Table and at https: //github. com/rcsb/symmetry-benchmark. CE-Symm was run on the RepeatsDB dataset to assess its performance on a larger set of symmetric proteins. It was able to detect repeats in 69% of the dataset of RepeatsDB reviewed entries. Considering only single-domain proteins improves the recall of 77%, indicating that multi-domain proteins are challenging for the method. The database classifies repeat regions according to the Kajava tandem repeat classes [36]. CE-Symm achieves a recall of 64% for single-domain solenoid repeats (class III) and 89% for closed repeats (class IV). Thus, CE-Symm performs best for closed repeats with single-domain input. The new CE-Symm tool is capable of presenting internal symmetry as a multiple structural alignment of repeats, which enables direct association of sequence and structure and can be used for comparative and evolutionary analyses of protein structures. The superposition used for the alignment is constrained by the axes of symmetry found in the structure, so that equivalent residue positions maintain the symmetric orientation. Symmetry-aware alignments are important to study, for example, binding and active sites at the internal symmetry interface, like calcium binding in βγ-crystallins (Fig 3). Structural alignment of the repeats can reveal conserved motifs that have persisted since the duplication event. One example is the βγ-crystallin superfamily, which occurs in a variety of repeat arrangements. Many βγ-crystallins contain a calcium binding site motif [37]. As shown in Fig 3A, the calcium binding motif is structurally conserved after a 2-fold rotation around the symmetry axis, and the residue side-chains preserve their orientation. Furthermore, calcium coordinates residues from both repeats, making the two-fold symmetry an essential feature of the binding site. On the other hand, duplication events allow the appearance of asymmetry by independent sequence and structural divergence of the repeats. An example is the MaoC-like thioesterase/thiol ester dehydrase-isomerase superfamily (SCOP: d. 38. 1. 4). Members of this family fold into a characteristic ‘hot dog fold’ which binds coenzyme A and catalyzes the dehydration of various bound fatty acids. Typically the MaoC-like proteins contain one hot dog domain per chain and assemble into dimers, tetramers, or hexamers [38]. Some members of the family contain a duplication of the hot dog fold [39], accompanied by the loss of the catalytic motif R/N-####-H in one of the domains, in order to accommodate bulkier substrates which would otherwise not fit in a single domain [38]. The structural divergence of the catalytic site in one of the repeats of the double hot dog subunit can be easily observed with CE-Symm (Fig 4). Some proteins contain more than one axis of symmetry. In those cases, the axes of symmetry can be collinear, orthogonal or independent to each other. If the axes are collinear, they can be combined into a single axis with higher symmetry order. If the axes are orthogonal, they can be combined into a point group of higher symmetry order. If the axes are independent to each other, multiple levels of symmetry exist in the structure in a hierarchical organization. This can be an indication of multiple independent duplication events, like in the case of γ-crystallins (Fig 3B), where four repeats are related by two independent 2-fold axes corresponding to two successive duplication events. Additionally, the internal symmetry axes can also combine with the quaternary symmetry axes. Therefore, internal symmetry can increase the order of symmetry of a protein complex. Returning to the previous MaoC-like protein example, the internal two-fold axis of the double hot dog domain in Fig 4B is orthogonal to the three-fold quaternary symmetry axis, combining for an overall dihedral symmetry. This arrangement is structurally similar to the D3 quaternary symmetry of hexameric single hot dog proteins (e. g. 1YLI). Accounting for internal symmetry when comparing the two protein assemblies is therefore important, because proteins can have a similar overall structure despite their different subunit compositions. It would be misleading to say that the structure of the trimeric and hexameric MaoC-like proteins are substantially different. Another well-know example of similar overall arrangement with different subunit composition are DNA clamps, which promote processivity in DNA replication. In archaea and eukaryotes, the clamp is a trimer, while in bacteria it is a dimer [40]. Furthermore, all DNA clamps have further internal symmetry axes leading to an overall D6 symmetry. As a historical note, the homology between bacterial and eurkaryotic DNA clamps was only acknowledged when the structures were solved and the similarity of their complexes was identified [41]. Furthermore, internal symmetry is important in understanding the stoichiometry of protein assemblies. Uneven stoichiometry assemblies are those with an unbalanced number of each entity type in the complex and occur rarely in the biological environment. It was previously reported that up to 40% of all protein assemblies with uneven stoichiometry in the PDB can be explained by the presence of internal symmetry in one or multiple of the subunits in the complex [42]. One such example is the artificial complex of Bowman-Birk inhibitor from snail medic seeds with bovine trypsin, which has an A2B stoichiometry (Fig 5). Although the complex is asymmetric, considering the internal symmetry of the inhibitor shows that the assembly is structurally comparable to an even A2B2 assembly with C2 overall symmetry. This property has also functional consequences, since the binding of two trypsin proteins symmetrically allows the inhibitor to efficiently induce dimerization and block the peptidase activity. Symmetry is characteristic of biological assemblies and can be considered by methods, like EPPIC, in order to predict the biological assembly in the context of crystal latices [43]. Including internal symmetry in these methods could further improve predictions for some known cases like, for example, uneven stoichiometries. The majority of proteins have closed symmetry. In the case of quaternary structures, this is expected since homooligomers with open symmetry are disfavored due to their aggregation potential [44]. However, this is not the case for internal symmetry due to the ability for terminal repeats to diverge and avoid undesirable homotypic interactions. The most general formulation of open repeats in the literature is that of superhelical symmetry, where the repeating unit is simultaneously translated along a helical path (curvature) and rotated around this path (twist) [28]. CE-Symm cannot identify superhelical symmetries, where both curvature and twist are relevant, because of the fundamental limit of the method to find a single symmetry axis (or multiple independent axes). However, we observe that the majority of structures containing tandem repeats that are classified as superhelical in the literature (solenoids) can be approximated with a single axis of symmetry by CE-Symm. They fall in one of the following four conditions: i) the twist is negligible (relative to CE-Symm tolerances); ii) the curvature is negligible; iii) both twist and curvature are negligible; or iv) the twist is much larger than the curvature. In all those cases, CE-Symm can identify the symmetry in the structures and annotate them as helical, translational or open rotational symmetries. For instance, from the 18 solenoid protein representatives from table 1 in Kobe and Kajava [28], in 10 either the twist or the curvature are reported to be small (helical symmetry applies), in 5 both the twist and curvature are annotated as small (translational symmetry applies), and the remaining 3 structure representatives have irregular twist (asymmetric applies). Although many folds are classified as superhelical, only a small number have regular repeats but do not fit into one of the above categories. Therefore, in practice CE-Symm can also be a good method for identifying, classifying and characterizing solenoid and other repeat proteins with open symmetry. We hypothesize that the low prevalence of actual superhelical symmetry in proteins could be a consequence of the benefit in conserving interfaces between adjacent repeats.
  We have extended our internal symmetry analysis tool in order to improve its usability, capabilities and the interpretability of results. In addition to detecting symmetry in protein structures, the tool can identify corresponding residues of the protein from each repeating element and the symmetry operations between them. CE-Symm 2. 0 adds broad capabilities for the detection of all types of internal symmetry, providing information about the type and order of symmetry and the repeat boundaries. The alignments between the repeats are eminently useful in identifying conserved and differential features between repeats, and can be applied to understanding protein function and evolution. The ability to run CE-Symm recursively to detect multiple axes of symmetry allows both higher-order point group symmetries to be identified and non-point group hierarchical symmetries. The simultaneous visualization of this rich information can lead to a better understanding of the structure and provide information about multiple duplication events. One limitation of CE-Symm 2. 0 is that it does not yet integrate quaternary symmetry detection into it’s hierarchy. While it is possible to run the program on biological assemblies, it will have poor performance and may mistakenly fail to detect symmetry. Rather, methods specific to quaternary symmetry detection should be integrated with CE-Symm to provide this feature. CE-Symm has been optimized for finding structures with global symmetry. While it does search for insertions, the length dependence of TM-Score means that structures with large insertions or multi-domain queries may not meet the default score thresholds. For multidomain proteins it may be needed to perform domain decomposition prior to running CE-Symm. Since it is based on CE rigid body alignment, the tool is also unlikely to detect all repeats in structures with conformational changes in some repeats, or with non-sequential rearrangements like circular permutations between repeats. Another limitation is the requirement for a consistent orientation between equivalent repeats. While for some applications preserving a conserved interface between repeats is desirable, there are many cases with large and functionally significant changes in repeat orientation (e. g. in many solenoid proteins). Determining whether the high prevalence of internal symmetry in protein structures is predominantly a consequence of thermodynamic selection or an indication of the history of protein evolution remains an open question. Here, we have presented examples where internal symmetry is a result of evolution and tied to functional consequences, and how our tool can help researchers in the protein evolution, classification and annotation fields. The CE-Symm source code has been integrated into the BioJava library [45] and is freely available on GitHub. In the future, we would like to integrate CE-Symm into leading bioinformatics resources for protein analysis.
summary: Many protein structures show a great deal of regularity. Even within single polypeptide chains, about 25% of proteins contain self-similar repeating structures, which can be organized in ring-like symmetric arrangements or linear open repeats. The repeats are often related, and thus comparing the sequence and structure of repeats can give an idea as to the early evolutionary history of a protein family. Additionally, the conservation and divergence of repeats can lead to insights about the function of the proteins. This work describes CE-Symm 2. 0, a tool for the analysis of protein symmetry. The method automatically detects internal symmetry in protein structures and produces a multiple alignment of structural repeats. The algorithm is able to detect the geometric relationships between the repeats, including cyclic, dihedral, and polyhedral symmetries, translational repeats, and cases where multiple symmetry operators are applicable in a hierarchical manner. These complex relationships can then be visualized in a graphical interface as a complete structure, as a superposition of repeats, or as a multiple alignment of the protein sequence. CE-Symm 2. 0 can be systematically used for the automatic detection of internal symmetry in protein structures, or as an interactive tool for the analysis of structural repeats.
keywords: split-decomposition method
  geometry
  multiple alignment calculation
  mathematics
  protein structure
  sequence motif analysis
  research and analysis methods
  sequence analysis
  sequence alignment
  bioinformatics
  proteins
  molecular evolution
  molecular biology
  protein structure comparison
  biochemistry
  computational techniques
  database and informatics methods
  protein domains
  biology and life sciences
  physical sciences
  evolutionary biology
  macromolecular structure analysis
  symmetry
title: Analyzing the symmetrical arrangement of structural repeats in proteins with CE-Symm