Rapid evolution in conformational space: a study of loop regions in a ubiquitous GTP binding domain

Generation of a Soluble African Horse Sickness Virus VP7 Protein Capable of Forming Core-like Particles

A unique characteristic of the African horse sickness virus (AHSV) major core protein VP7 is that it is highly insoluble, and spontaneously forms crystalline particles in AHSV-infected cells and when expressed in vitro. The aggregation of AHSV VP7 into these crystals presents many problems in AHSV vaccine development, and it is unclear whether VP7 aggregation affects AHSV assembly or contributes to AHSV pathogenesis. Here, we set out to abolish VP7 self-assembly by targeting candidate amino acid regions on the surface of the VP7 trimer via site-directed mutagenesis. It was found that the substitution of seven amino acids resulted in the complete disruption of AHSV VP7 self-assembly, which abolished the formation of VP7 crystalline particles and converted VP7 to a fully soluble protein still capable of interacting with VP3 to form core-like particles. This work provides further insight into the formation of AHSV VP7 crystalline particles and the successful development of AHSV vaccines. It also paves the way for future research by drawing comparisons with similar viral phenomena observed in human virology.

Multidimensional Phylogenetic Metrics Identify Class I Aminoacyl-tRNA Synthetase Evolutionary Mosaicity and Inter-Modular Coupling

The role of aminoacyl-tRNA synthetases (aaRS) in the emergence and evolution of genetic coding poses challenging questions concerning their provenance. We seek evidence about their ancestry from curated structure-based multiple sequence alignments of a structurally invariant “scaffold” shared by all 10 canonical Class I aaRS. Three uncorrelated phylogenetic metrics—mutation frequency, its uniformity, and row-by-row cladistic congruence—imply that the Class I scaffold is a mosaic assembled from successive genetic sources. Metrics for different modules vary in accordance with their presumed functionality. Sequences derived from the ATP– and amino acid– binding sites exhibit specific two-way coupling to those derived from Connecting Peptide 1, a third module whose metrics suggest later acquisition. The data help validate: (i) experimental fragmentations of the canonical Class I structure into three partitions that retain catalytic activities in proportion to their length; and (ii) evidence that the ancestral Class I aaRS gene also encoded a Class II ancestor in frame on the opposite strand. A 46-residue Class I “protozyme” roots the Class I tree prior to the adaptive radiation of the Rossmann dinucleotide binding fold that refined substrate discrimination. Such rooting implies near simultaneous emergence of genetic coding and the origin of the proteome, resolving a conundrum posed by previous inferences that Class I aaRS evolved after the genetic code had been implemented in an RNA world. Further, pinpointing discontinuous enhancements of aaRS fidelity establishes a timeline for the growth of coding from a binary amino acid alphabet.

Intergenic ORFs as elementary structural modules of de novo gene birth and protein evolution

The noncoding genome plays an important role in de novo gene birth and in the emergence of genetic novelty. Nevertheless, how noncoding sequences' properties could promote the birth of novel genes and shape the evolution and the structural diversity of proteins remains unclear. Therefore, by combining different bioinformatic approaches, we characterized the fold potential diversity of the amino acid sequences encoded by all intergenic open reading frames (ORFs) of S. cerevisiae with the aim of (1) exploring whether the structural states' diversity of proteomes is already present in noncoding sequences, and (2) estimating the potential of the noncoding genome to produce novel protein bricks that could either give rise to novel genes or be integrated into pre-existing proteins, thus participating in protein structure diversity and evolution. We showed that amino acid sequences encoded by most yeast intergenic ORFs contain the elementary building blocks of protein structures. Moreover, they encompass the large structural state diversity of canonical proteins, with the majority predicted as foldable. Then, we investigated the early stages of de novo gene birth by reconstructing the ancestral sequences of 70 yeast de novo genes and characterized the sequence and structural properties of intergenic ORFs with a strong translation signal. This enabled us to highlight sequence and structural factors determining de novo gene emergence. Finally, we showed a strong correlation between the fold potential of de novo proteins and one of their ancestral amino acid sequences, reflecting the relationship between the noncoding genome and the protein structure universe.

Diversity of sialidases found in the human body - A review

Sialidases are enzymes essential for numerous organisms including humans. Hydrolytic sialidases (EC 3.2.1.18), trans-sialidases and anhydrosialidases (intramolecular trans-sialidases, EC 4.2.2.15) are glycoside hydrolase enzymes that cleave the glycosidic linkage and release sialic acid residues from sialyl substrates. The paper summarizes diverse sialidases present in the human body and their potential impact on development of antiviral compounds - inhibitors of viral neuraminidases. It includes a brief overview of catalytic mechanisms of action of sialidases and describes the origin of sialidases in the human body. This is followed by description of the structure and function of sialidase families with a special focus on the GH33 and GH34 families. Various effects of sialidases on human body are also briefly described. Modulation of sialidase activity may be considered a useful tool for effective treatment of various diseases. In some cases, it is desired to completely suppress the activity of sialidases by suitable inhibitors. Specific sialidase inhibitors are useful for the treatment of influenza, epilepsy, Alzheimer's disease, diabetes, different types of cancer, or heart defects. Challenges and future directions are shortly depicted in the final part of the paper.

The interdigitating loop of the enolase superfamily as a specificity binding determinant or 'flying buttress'

BACKGROUND

Enzymes of the enolase superfamily (ENS) are mechanistically diverse, yet share a common partial reaction (abstraction of the α-proton from a carboxylate substrate). While the catalytic machinery responsible for the deprotonation reaction has been conserved, divergent evolution has led to numerous ENS members that catalyze different overall reactions. This rich functional diversity has made the ENS an excellent model system for developing the approaches necessary to validate enzyme function. However, enzymes of the ENS also share a common bidomain structure ((β/α)₇β-barrel domain and α+β capping domain) which makes validation of function from structural information challenging.

SCOPE OF THE REVIEW

This review presents a comparative survey of the structural data obtained over the past decade for enzymes from all seven subgroups that comprise the ENS.

MAJOR CONCLUSIONS

Of the seven ENS subgroups (enolase, mandelate racemase (MR), muconate lactonizing enzyme, β-methylaspartate ammonia lyase, d-glucarate dehydratase, d-mannonate dehydratase (ManD), and galactarate dehydratase 2), only enzymes of the MR and ManD subgroups exhibit an additional feature of structural complexity-an interdigitating loop. This loop emanates from one protomer of a homodimeric pair and penetrates into the adjacent, symmetry-related protomer to either contribute a binding determinant to the active site of the adjacent protomer, or act as a 'flying buttress' to support residues of the active site.

GENERAL SIGNIFICANCE

The analysis presented in this review suggests that the interdigitating loop is the only gross structural element that permits functional distinction between ENS subgroups at the tertiary level of protein structure.

Three-dimensional (3D) structure prediction of the American and African oil-palms β-ketoacyl-[ACP] synthase-II protein by comparative modelling

BACKGROUND

The fatty-acid profile of the vegetable oils determines its properties and nutritional value. Palm-oil obtained from the African oil-palm [Elaeis guineensis Jacq. (Tenera)] contains 44% palmitic acid (C16:0), but, palm-oil obtained from the American oilpalm [Elaeis oleifera] contains only 25% C16:0. In part, the b-ketoacyl-[ACP] synthase II (KASII) [EC: 2.3.1.179] protein is responsible for the high level of C16:0 in palm-oil derived from the African oil-palm. To understand more about E. guineensis KASII (EgKASII) and E. oleifera KASII (EoKASII) proteins, it is essential to know its structures. Hence, this study was undertaken.

OBJECTIVE

The objective of this study was to predict three-dimensional (3D) structure of EgKASII and EoKASII proteins using molecular modelling tools.

MATERIALS AND METHODS

The amino-acid sequences for KASII proteins were retrieved from the protein database of National Center for Biotechnology Information (NCBI), USA. The 3D structures were predicted for both proteins using homology modelling and ab-initio technique approach of protein structure prediction. The molecular dynamics (MD) simulation was performed to refine the predicted structures. The predicted structure models were evaluated and root mean square deviation (RMSD) and root mean square fluctuation (RMSF) values were calculated.

RESULTS

The homology modelling showed that EgKASII and EoKASII proteins are 78% and 74% similar with Streptococcus pneumonia KASII and Brucella melitensis KASII, respectively. The EgKASII and EoKASII structures predicted by using ab-initio technique approach shows 6% and 9% deviation to its structures predicted by homology modelling, respectively. The structure refinement and validation confirmed that the predicted structures are accurate.

CONCLUSION

The 3D structures for EgKASII and EoKASII proteins were predicted. However, further research is essential to understand the interaction of EgKASII and EoKASII proteins with its substrates.

A proline-rich loop mediates specific functions of human sialidase NEU4 in SK-N-BE neuronal differentiation

Human sialidase NEU4 long (N4L) is a membrane-associated enzyme that has been shown to be localized in the outer mitochondrial membrane. A role in different cellular processes has been suggested for this enzyme, such as apoptosis, neuronal differentiation and tumorigenesis. However, the molecular bases for these roles, not found in any of the other highly similar human sialidases, are not understood. We have found that a proline-rich sequence of 81 amino acids, unique to NEU4 sequence, contains potential Akt and Erk1 kinase motifs. Molecular modeling, based on the experimentally determined three-dimensional structure of cytosolic human NEU2, showed that the proline-rich sequence is accommodated in a loop, thus preserving the typical beta-barrel structure of sialidases. In order to investigate the role of this loop in neuronal differentiation, we obtained SK-N-BE neuroblastoma cells stably overexpressing either human wild-type N4L or a deletion mutant lacking the proline-rich loop. Our results demonstrate that the proline-rich region can also enhance cell proliferation and retinoic acid (RA)-induced neuronal differentiation and it is also involved in NEU4 interaction with Akt, as well as in substrate recognition, modifying directly or through the interaction with other protein(s) the enzyme specificity toward sialylated glycoprotein(s). On the whole, our results suggest that N4L could be a downstream component of the PI3K/Akt signaling pathway required for RA-induced differentiation of neuroblastoma SK-N-BE cells.

A role for indels in the evolution of Cro protein folds

Insertions and deletions in protein sequences, or indels, can disrupt structure and may result in changes in protein folds during evolution or in association with alternative splicing. Pfl 6 and Xfaso 1 are two proteins in the Cro family that share a common ancestor but have different folds. Sequence alignments of the two proteins show two gaps, one at the N terminus, where the sequence of Xfaso 1 is two residues shorter, and one near the center of the sequence, where the sequence of Pfl 6 is five residues shorter. To test the potential importance of indels in Cro protein evolution, we generated hybrid variants of Pfl 6 and Xfaso 1 with indels in one or both regions, chosen according to several plausible sequence alignments. All but one deletion variant completely unfolded both proteins, showing that a longer N-terminal sequence was critical for Pfl 6 folding and a longer central region sequence was critical for Xfaso 1 folding. By contrast, Xfaso 1 tolerated a longer N-terminal sequence with little destabilization, and Pfl 6 tolerated central region insertions, albeit with substantial effects on thermal stability and some perturbation of the surrounding structure. None of the mutations appeared to convert one stable fold into the other. On the basis of this two-protein comparison, short insertion and deletion mutations probably played a role in evolutionary fold change in the Cro family, but were also not the only factors.

A modular perspective of protein structures: application to fragment based loop modeling

Proteins can be decomposed into supersecondary structure modules. We used a generic definition of supersecondary structure elements, so-called Smotifs, which are composed of two flanking regular secondary structures connected by a loop, to explore the evolution and current variety of structure building blocks. Here, we discuss recent observations about the saturation of Smotif geometries in protein structures and how it opens new avenues in protein structure modeling and design. As a first application of these observations we describe our loop conformation modeling algorithm, ArchPred that takes advantage of Smotifs classification. In this application, instead of focusing on specific loop properties the method narrows down possible template conformations in other, often not homologous structures, by identifying the most likely supersecondary structure environment that cradles the loop. Beyond identifying the correct starting supersecondary structure geometry, it takes into account information of fit of anchor residues, sterical clashes, match of predicted and observed dihedral angle preferences, and local sequence signal.

New insights on the sialidase protein family revealed by a phylogenetic analysis in metazoa

Sialidases are glycohydrolytic enzymes present from virus to mammals that remove sialic acid from oligosaccharide chains. Four different sialidase forms are known in vertebrates: the lysosomal NEU1, the cytosolic NEU2 and the membrane-associated NEU3 and NEU4. These enzymes modulate the cell sialic acid content and are involved in several cellular processes and pathological conditions. Molecular defects in NEU1 are responsible for sialidosis, an inherited disease characterized by lysosomal storage disorder and neurodegeneration. The studies on the biology of sialic acids and sialyltransferases, the anabolic counterparts of sialidases, have revealed a complex picture with more than 50 sialic acid variants selectively present in the different branches of the tree of life. The gain/loss of specific sialoconjugates have been proposed as key events in the evolution of deuterostomes and Homo sapiens, as well as in the host-pathogen interactions. To date, less attention has been paid to the evolution of sialidases. Thus we have conducted a survey on the state of the sialidase family in metazoan. Using an in silico approach, we identified and characterized sialidase orthologs from 21 different organisms distributed among the evolutionary tree: Metazoa relative (Monosiga brevicollis), early Deuterostomia, precursor of Chordata and Vertebrata (teleost fishes, amphibians, reptiles, avians and early and recent mammals). We were able to reconstruct the evolution of the sialidase protein family from the ancestral sialidase NEU1 and identify a new form of the enzyme, NEU5, representing an intermediate step in the evolution leading to the modern NEU3, NEU4 and NEU2. Our study provides new insights on the mechanisms that shaped the substrate specificity and other peculiar properties of the modern mammalian sialidases. Moreover, we further confirm findings on the catalytic residues and identified enzyme loop portions that behave as rapidly diverging regions and may be involved in the evolution of specific properties of sialidases.

IndelFR: a database of indels in protein structures and their flanking regions

Insertion/deletion (indel) is one of the most common methods of protein sequence variation. Recent studies showed that indels could affect their flanking regions and they are important for protein function and evolution. Here, we describe the Indel Flanking Region Database (IndelFR, http://indel.bioinfo.sdu.edu.cn), which provides sequence and structure information about indels and their flanking regions in known protein domains. The indels were obtained through the pairwise alignment of homologous structures in SCOP superfamilies. The IndelFR database contains 2,925,017 indels with flanking regions extracted from 373,402 structural alignment pairs of 12,573 non-redundant domains from 1053 superfamilies. IndelFR provides access to information about indels and their flanking regions, including amino acid sequences, lengths, locations, secondary structure constitutions, hydrophilicity/hydrophobicity, domain information, 3D structures and so on. IndelFR has already been used for molecular evolution studies and may help to promote future functional studies of indels and their flanking regions.

Gallus gallus NEU3 sialidase as model to study protein evolution mechanism based on rapid evolving loops

BACKGROUND

Large surface loops contained within compact protein structures and not involved in catalytic process have been proposed as preferred regions for protein family evolution. These loops are subjected to lower sequence constraints and can evolve rapidly in novel structural variants. A good model to study this hypothesis is represented by sialidase enzymes. Indeed, the structure of sialidases is a β-propeller composed by anti-parallel β-sheets connected by loops that suit well with the rapid evolving loop hypothesis. These features prompted us to extend our studies on this protein family in birds, to get insights on the evolution of this class of glycohydrolases.

RESULTS

Gallus gallus (Gg) genome contains one NEU3 gene encoding a protein with a unique 188 amino acid sequence mainly constituted by a peptide motif repeated six times in tandem with no homology with any other known protein sequence. The repeat region is located at the same position as the roughly 80 amino acid loop characteristic of mammalian NEU4. Based on molecular modeling, all these sequences represent a connecting loop between the first two highly conserved β-strands of the fifth blade of the sialidase β-propeller. Moreover this loop is highly variable in sequence and size in NEU3 sialidases from other vertebrates. Finally, we found that the general enzymatic properties and subcellular localization of Gg NEU3 are not influenced by the deletion of the repeat sequence.

CONCLUSION

In this study we demonstrated that sialidase protein structure contains a surface loop, highly variable both in sequence and size, connecting two conserved β-sheets and emerging on the opposite site of the catalytic crevice. These data confirm that sialidase family can serve as suitable model for the study of the evolutionary process based on rapid evolving loops, which may had occurred in sialidases. Giving the peculiar organization of the loop region identified in Gg NEU3, this protein can be considered of particular interest in such evolutionary studies and to get deeper insights in sialidase evolution.

Impact of indels on the flanking regions in structural domains

Amino acid substitution and insertions/deletions (indels) are two common events in protein evolution; however, current knowledge on indels is limited. In this study, we investigated the effects of indels on the flanking regions in protein structure superfamilies. Comprehensive analysis of structural classification of proteins superfamilies revealed that indels lead to a series of changes in the flanking regions, including the following: 1) structural shift in the tertiary structure, with a first-order exponential decay relation between structural shift and the distance to indels, 2) instability of the secondary structure elements in which parts of the α helix and β sheet are destroyed, and 3) an increase in the amino acid substitution rate of the primary structure and the nonsimilar amino acid substitution rate. In general, these quality changes are due to the combined effects of the "regional-inherent effect," "indel-accompanied effect," and "indel-following effect." Furthermore, these quality changes reflect changes in selective pressure. Indels are more likely to be preserved in regions with low selective pressure, and indels can further reduce the selective pressure on the flanking regions. These findings improve our understanding of the role of indels in protein evolution.

Evolutionary pattern of protein architecture in mammal and fruit fly genomes

Mutations, which can alter amino acid constitution, contribute greatly to protein evolution. However, little is reported of their pattern during protein structural evolution. We investigated the distribution of non-synonymous single nucleotide polymorphisms (nsSNPs) and insertions/deletions (indels) along mammal and fruit fly proteins. We found the nsSNPs (and d(N)) and indels increased in protein boundary regions, and this pattern is inversely correlated with the distribution of protein domain density. Additionally, synonymous substitutions (and d(S)) are reduced in 5' and 3' regions, indicating more variable protein boundaries, compared with central interior. All evidence suggests that the inner part of coding sequences (CDSs) is comparatively conserved, whereas the 5' and 3' regions, with higher evolution rates, are more variable. We assumed that due to greater frequencies of nsSNPs and indels in adaptive regions of CDSs it could be easier to ultimately alter, gain, or lose amino acids, thus becoming the front line of protein evolution.

Evolution of specificity in the eukaryotic endomembrane system

Two hundred years after Darwin's birth, our understanding of genetic mechanisms and cell biology has advanced to a level unimaginable in the 19th century. We now know that eukaryotic cells contain a huge variety of internal compartments, each with their own function, identity and history. For the compartments that together form the membrane-trafficking system, one of the central questions is how that identity is encoded and how it evolved. Here we review the key components involved in membrane-trafficking events, including SNAREs, Rabs, vesicle coats, and tethers and what is known about their evolutionary history. Our current understanding suggests a possible common mechanism by which the membrane-trafficking organelles might have evolved. This model of increased organellar complexity by gene duplication and co-evolution of multiple, interacting, specificity-encoding proteins could well be applicable to other non-endosymbiotic organelles as well. The application of basic evolutionary principles well beyond their original scope has been exceedingly powerful not only in reconstructing the history of cellular compartments, but for medical and applied research as well, and underlines the contributions of Darwin's ideas in modern biology.

Identification and evolution of fungal mitochondrial tyrosyl-tRNA synthetases with group I intron splicing activity

The bifunctional Neurospora crassa mitochondrial tyrosyl-tRNA synthetase (CYT-18 protein) both aminoacylates mitochondrial tRNA(Tyr) and acts as a structure-stabilizing splicing cofactor for group I introns. Previous studies showed that CYT-18 has distinct tRNA(Tyr) and group I intron-binding sites, with the latter formed by three small "insertions" in the nucleotide-binding fold and other structural adaptations compared with nonsplicing bacterial tyrosyl-tRNA synthetases. Here, analysis of genomic sequences shows that mitochondrial tyrosyl-tRNA synthetases with structural adaptations similar to CYT-18's are uniquely characteristic of fungi belonging to the subphylum Pezizomycotina, and biochemical assays confirm group I intron splicing activity for the enzymes from several of these organisms, including Aspergillus nidulans and the human pathogens Coccidioides posadasii and Histoplasma capsulatum. By combining multiple sequence alignments with a previously determined cocrystal structure of a CYT-18/group I intron RNA complex, we identify conserved features of the Pezizomycotina enzymes related to group I intron and tRNA interactions. Our results suggest that mitochondrial tyrosyl-tRNA synthetases with group I intron splicing activity evolved during or after the divergence of the fungal subphyla Pezizomycotina and Saccharomycotina by a mechanism involving the concerted differentiation of preexisting protein loop regions. The unique group I intron splicing activity of these fungal enzymes may provide a new target for antifungal drugs.

Insertions and the emergence of novel protein structure: a structure-based phylogenetic study of insertions

BACKGROUND

In protein evolution, the mechanism of the emergence of novel protein domain is still an open question. The incremental growth of protein variable regions, which was produced by stochastic insertions, has the potential to generate large and complex sub-structures. In this study, a deterministic methodology is proposed to reconstruct phylogenies from protein structures, and to infer insertion events in protein evolution. The analysis was performed on a broad range of SCOP domain families.

RESULTS

Phylogenies were reconstructed from protein 3D structural data. The phylogenetic trees were used to infer ancestral structures with a consensus method. From these ancestral reconstructions, 42.7% of the observed insertions are nested insertions, which locate in previous insert regions. The average size of inserts tends to increase with the insert rank or total number of insertions in the variable regions. We found that the structures of some nested inserts show complex or even domain-like fold patterns with helices, strands and loops. Furthermore, a basal level of structural innovation was found in inserts which displayed a significant structural similarity exclusively to themselves. The beta-Lactamase/D-ala carboxypeptidase domain family is provided as an example to illustrate the inference of insertion events, and how the incremental growth of a variable region is capable to generate novel structural patterns.

CONCLUSION

Using 3D data, we proposed a method to reconstruct phylogenies. We applied the method to reconstruct the sequences of insertion events leading to the emergence of potentially novel structural elements within existing protein domains. The results suggest that structural innovation is possible via the stochastic process of insertions and rapid evolution within variable regions where inserts tend to be nested. We also demonstrate that the structure-based phylogeny enables the study of new questions relating to the evolution of protein domain and biological function.

A supersecondary structure library and search algorithm for modeling loops in protein structures

We present a fragment-search based method for predicting loop conformations in protein models. A hierarchical and multidimensional database has been set up that currently classifies 105 950 loop fragments and loop flanking secondary structures. Besides the length of the loops and types of bracing secondary structures the database is organized along four internal coordinates, a distance and three types of angles characterizing the geometry of stem regions. Candidate fragments are selected from this library by matching the length, the types of bracing secondary structures of the query and satisfying the geometrical restraints of the stems and subsequently inserted in the query protein framework where their fit is assessed by the root mean square deviation (r.m.s.d.) of stem regions and by the number of rigid body clashes with the environment. In the final step remaining candidate loops are ranked by a Z-score that combines information on sequence similarity and fit of predicted and observed ϕ/ψ main chain dihedral angle propensities. Confidence Z-score cut-offs were determined for each loop length that identify those predicted fragments that outperform a competitive ab initio method. A web server implements the method, regularly updates the fragment library and performs prediction. Predicted segments are returned, or optionally, these can be completed with side chain reconstruction and subsequently annealed in the environment of the query protein by conjugate gradient minimization. The prediction method was tested on artificially prepared search datasets where all trivial sequence similarities on the SCOP superfamily level were removed. Under these conditions it is possible to predict loops of length 4, 8 and 12 with coverage of 98, 78 and 28% with at least of 0.22, 1.38 and 2.47 Å of r.m.s.d. accuracy, respectively. In a head-to-head comparison on loops extracted from freshly deposited new protein folds the current method outperformed in a ∼5:1 ratio an earlier developed database search method.

Analysis of protein homology by assessing the (dis)similarity in protein loop regions

Two proteins are considered to have a similar fold if sufficiently many of their secondary structure elements are positioned similarly in space and are connected in the same order. Such a common structural scaffold may arise due to either divergent or convergent evolution. The intervening unaligned regions ("loops") between the superimposable helices and strands can exhibit a wide range of similarity and may offer clues to the structural evolution of folds. One might argue that more closely related proteins differ less in their nonconserved loop regions than distantly related proteins and, at the same time, the degree of variability in the loop regions in structurally similar but unrelated proteins is higher than in homologs. Here we introduce a new measure for structural (dis)similarity in loop regions that is based on the concept of the Hausdorff metric. This measure is used to gauge protein relatedness and is tested on a benchmark of homologous and analogous protein structures. It has been shown that the new measure can distinguish homologous from analogous proteins with the same or higher accuracy than the conventional measures that are based on comparing proteins in structurally aligned regions. We argue that this result can be attributed to the higher sensitivity of the Hausdorff (dis)similarity measure in detecting particularly evident dissimilarities in structures and draw some conclusions about evolutionary relatedness of proteins in the most populated protein folds.