Structural convergence during protein evolution

The three-dimensional structure of the Vint domain from Tetrahymena thermophila suggests a ligand-regulated cleavage mechanism by the HINT fold

Vint proteins have been identified in unicellular metazoans as a novel hedgehog-related gene family, merging the von Willebrand factor type A domain and the Hedgehog/INTein (HINT) domains. We present the first three-dimensional structure of the Vint domain from Tetrahymena thermophila corresponding to the auto-processing domain of hedgehog proteins, shedding light on the unique features, including an adduct recognition region (ARR). Our results suggest a potential binding between the ARR and sulfated glycosaminoglycans like heparin sulfate. Moreover, we uncover a possible regulatory role of the ARR in the auto-processing by Vint domains, expanding our understanding of the HINT domain evolution and their use in biotechnological applications. Vint domains might have played a crucial role in the transition from unicellular to multicellular organisms.

Using AlphaFold Predictions in Viral Research

Elucidation of the tertiary structure of proteins is an important task for biological and medical studies. AlphaFold, a modern deep-learning algorithm, enables the prediction of protein structure to a high level of accuracy. It has been applied in numerous studies in various areas of biology and medicine. Viruses are biological entities infecting eukaryotic and procaryotic organisms. They can pose a danger for humans and economically significant animals and plants, but they can also be useful for biological control, suppressing populations of pests and pathogens. AlphaFold can be used for studies of molecular mechanisms of viral infection to facilitate several activities, including drug design. Computational prediction and analysis of the structure of bacteriophage receptor-binding proteins can contribute to more efficient phage therapy. In addition, AlphaFold predictions can be used for the discovery of enzymes of bacteriophage origin that are able to degrade the cell wall of bacterial pathogens. The use of AlphaFold can assist fundamental viral research, including evolutionary studies. The ongoing development and improvement of AlphaFold can ensure that its contribution to the study of viral proteins will be significant in the future.

Use of an Integrated Approach Involving AlphaFold Predictions for the Evolutionary Taxonomy of Duplodnaviria Viruses

The evaluation of the evolutionary relationships is exceptionally important for the taxonomy of viruses, which is a rapidly expanding area of research. The classification of viral groups belonging to the realm Duplodnaviria, which include tailed bacteriophages, head-tailed archaeal viruses and herpesviruses, has undergone many changes in recent years and continues to improve. One of the challenging tasks of Duplodnaviria taxonomy is the classification of high-ranked taxa, including families and orders. At the moment, only 17 of 50 families have been assigned to orders. The evaluation of the evolutionary relationships between viruses is complicated by the high level of divergence of viral proteins. However, the development of structure prediction algorithms, including the award-winning AlphaFold, encourages the use of the results of structural predictions to clarify the evolutionary history of viral proteins. In this study, the evolutionary relationships of two conserved viral proteins, the major capsid protein and terminase, representing different viruses, including all classified Duplodnaviria families, have been analysed using AlphaFold modelling. This analysis has been undertaken using structural comparisons and different phylogenetic methods. The results of the analyses mainly indicated the high quality of AlphaFold modelling and the possibility of using the AlphaFold predictions, together with other methods, for the reconstruction of the evolutionary relationships between distant viral groups. Based on the results of this integrated approach, assumptions have been made about refining the taxonomic classification of bacterial and archaeal Duplodnaviria groups, and problems relating to the taxonomic classification of Duplodnaviria have been discussed.

PR10/Bet v1-like Proteins as Novel Contributors to Plant Biochemical Diversity

Pathogenesis-related (PR) proteins constitute a broad class of plant proteins with analogues found throughout nature from bacteria to higher eukaryotes. PR proteins were first noted in plants as part of the hypersensitive response, but have since been assigned an array of biological roles. The PR10/Bet v1-like proteins are a subset of PR proteins characterized by an ability to bind a wide range of lipophilic ligands, uniquely positioning them as contributors to specialized biosynthetic pathways. PR10/Bet v1-like proteins participate in the production of plant alkaloids and phenolics including flavonoids, both as general binding proteins and in special cases as catalysts. Owing initially to the perceived allergenic properties of PR10/Bet v1-like proteins, many were studied at the structural level to elucidate the basis for ligand binding. These studies provided a foundation for more recent efforts to understand higher-level structural order and how PR10/Bet v1-like proteins catalyse key reactions in plant pathways. Synthetic biology aimed at reconstituting plant-specialized metabolism in microorganisms uses knowledge of these proteins to fine-tune performance in new systems.

Maintenance of a Protein Structure in the Dynamic Evolution of TIMPs over 600 Million Years

Deciphering the events leading to protein evolution represents a challenge, especially for protein families showing complex evolutionary history. Among them, TIMPs represent an ancient eukaryotic protein family widely distributed in the animal kingdom. They are known to control the turnover of the extracellular matrix and are considered to arise early during metazoan evolution, arguably tuning essential features of tissue and epithelial organization. To probe the structure and molecular evolution of TIMPs within metazoans, we report the mining and structural characterization of a large data set of TIMPs over approximately 600 Myr. The TIMPs repertoire was explored starting from the Cnidaria phylum, coeval with the origins of connective tissue, to great apes and humans. Despite dramatic sequence differences compared with highest metazoans, the ancestral proteins displayed the canonical TIMP fold. Only small structural changes, represented by an α-helix located in the N-domain, have occurred over the evolution. Both the occurrence of such secondary structure elements and the relative solvent accessibility of the corresponding residues in the three-dimensional structures raises the possibility that these sites represent unconserved element prone to accept variations.

True arrestins and arrestin-fold proteins: a structure-based appraisal

Arrestin-clan proteins are folded alike, a feature responsible for their recent grouping in a single clan. In human, it includes the well-characterized visual and β-arrestins, the arrestin domain-containing proteins (ARRDCs), isoforms of the retromer subunit VPS26, and DSCR3, a protein involved in Down syndrome. A new arrestin-fold-predicted protein, RGP1, described here may join the clan. Unicellular organisms like the yeast Saccharomyces cerevisiae or the amoeba Dictyostelium discoideum harbor VPS26, DSCR3, and RGP1 isoforms as well as arrestin-related trafficking adaptors or ADCs, but true arrestins are missing. Functionally, members of the arrestin clan have generally a scaffolding role in various membrane protein trafficking events. Despite their similar structure, the mechanism of cargo recognition and internalization and the nature of recruited partners differ for the different members. Based on the recent literature, true arrestins (visual and β-arrestins), ARRDCs, and yeast ARTS are the closest from a functional point of view.

Local structural differences in homologous proteins: specificities in different SCOP classes

The constant increase in the number of solved protein structures is of great help in understanding the basic principles behind protein folding and evolution. 3-D structural knowledge is valuable in designing and developing methods for comparison, modelling and prediction of protein structures. These approaches for structure analysis can be directly implicated in studying protein function and for drug design. The backbone of a protein structure favours certain local conformations which include α-helices, β-strands and turns. Libraries of limited number of local conformations (Structural Alphabets) were developed in the past to obtain a useful categorization of backbone conformation. Protein Block (PB) is one such Structural Alphabet that gave a reasonable structure approximation of 0.42 Å. In this study, we use PB description of local structures to analyse conformations that are preferred sites for structural variations and insertions, among group of related folds. This knowledge can be utilized in improving tools for structure comparison that work by analysing local structure similarities. Conformational differences between homologous proteins are known to occur often in the regions comprising turns and loops. Interestingly, these differences are found to have specific preferences depending upon the structural classes of proteins. Such class-specific preferences are mainly seen in the all-β class with changes involving short helical conformations and hairpin turns. A test carried out on a benchmark dataset also indicates that the use of knowledge on the class specific variations can improve the performance of a PB based structure comparison approach. The preference for the indel sites also seem to be confined to a few backbone conformations involving β-turns and helix C-caps. These are mainly associated with short loops joining the regular secondary structures that mediate a reversal in the chain direction. Rare β-turns of type I’ and II’ are also identified as preferred sites for insertions.

SitEx: a computer system for analysis of projections of protein functional sites on eukaryotic genes

Search of interrelationships between the structural-functional protein organization and exon structure of encoding gene provides insights into issues concerned with the function, origin and evolution of genes and proteins. The functions of proteins and their domains are defined mostly by functional sites. The relation of the exon-intron structure of the gene to the protein functional sites has been little studied. Development of resources containing data on projections of protein functional sites on eukaryotic genes is needed. We have developed SitEx, a database that contains information on functional site amino acid positions in the exon structure of encoding gene. SitEx is integrated with the BLAST and 3DExonScan programs. BLAST is used for searching sequence similarity between the query protein and polypeptides encoded by single exons stored in SitEx. The 3DExonScan program is used for searching for structural similarity of the given protein with these polypeptides using superimpositions. The developed computer system allows users to analyze the coding features of functional sites by taking into account the exon structure of the gene, to detect the exons involved in shuffling in protein evolution, also to design protein-engineering experiments. SitEx is accessible at http://www-bionet.sscc.ru/sitex/. Currently, it contains information about 9994 functional sites presented in 2021 proteins described in proteomes of 17 organisms.

Homology versus analogy: possible evolutionary relationship of immunoglobulins, cupredoxins, and Cu,Zn-superoxide dismutase

The 'immunoglobulin-like' fold is one of most common structural motifs observed in proteins. This topology is found in more than 80 superfamilies of proteins, including Cu,Zn-superoxide dismutase (SOD) and cupredoxin. Evolutionary relationships have not been identified, but may exist. The challenge remains, therefore, of resolving the issue of whether the diverse distribution of the fold is accounted for by divergent evolution of function or convergent evolution of structure following multiple independent origins of function. Since the early studies that revealed conformational similarity of immunoglobulins and other proteins, the number of primary structures available for comparison has dramatically increased and new computational approaches for analysis of sequences have been developed. It now appears that a hypothesis of a common evolutionary origin for cupredoxins, Cu,Zn-SOD, and immunoglobulins may be credible. The distinction between protein homology and protein analogy is fundamental. The immunoglobulin-like fold may represent a robust system within which to examine again the issue of protein homology versus analogy.

Optimality of the genetic code with respect to protein stability and amino-acid frequencies

BACKGROUND

The genetic code is known to be efficient in limiting the effect of mistranslation errors. A misread codon often codes for the same amino acid or one with similar biochemical properties, so the structure and function of the coded protein remain relatively unaltered. Previous studies have attempted to address this question quantitatively, by estimating the fraction of randomly generated codes that do better than the genetic code in respect of overall robustness. We extended these results by investigating the role of amino-acid frequencies in the optimality of the genetic code.

RESULTS

We found that taking the amino-acid frequency into account decreases the fraction of random codes that beat the natural code. This effect is particularly pronounced when more refined measures of the amino-acid substitution cost are used than hydrophobicity. To show this, we devised a new cost function by evaluating in silico the change in folding free energy caused by all possible point mutations in a set of protein structures. With this function, which measures protein stability while being unrelated to the code's structure, we estimated that around two random codes in a billion (109) are fitter than the natural code. When alternative codes are restricted to those that interchange biosynthetically related amino acids, the genetic code appears even more optimal.

CONCLUSIONS

These results lead us to discuss the role of amino-acid frequencies and other parameters in the genetic code's evolution, in an attempt to propose a tentative picture of primitive life.

Second codon positions of genes and the secondary structures of proteins. Relationships and implications for the origin of the genetic code

The nucleotide frequencies in the second codon positions of genes are remarkably different for the coding regions that correspond to different secondary structures in the encoded proteins, namely, helix, beta-strand and aperiodic structures. Indeed, hydrophobic and hydrophilic amino acids are encoded by codons having U or A, respectively, in their second position. Moreover, the beta-strand structure is strongly hydrophobic, while aperiodic structures contain more hydrophilic amino acids. The relationship between nucleotide frequencies and protein secondary structures is associated not only with the physico-chemical properties of these structures but also with the organisation of the genetic code. In fact, this organisation seems to have evolved so as to preserve the secondary structures of proteins by preventing deleterious amino acid substitutions that could modify the physico-chemical properties required for an optimal structure.

On the classes of aminoacyl-tRNA synthetases and the error minimization in the genetic code

As a consequence of the existence of two classes of aminoacyl-tRNA synthetases (aaRSs), we defined two types of mutations: g (mutations that do not change the class of the involved amino acids) and u (those which change the class). We have found that the mean chemical distance resulting from g mutations is smaller than that corresponding to u mutations, indicating that g mutations are responsible for most of the known minimization of the genetic code. This supports models for the origin and evolution of the code, in which new amino acids were added after duplications or modification of existing aaRSs.

Evolution of the genetic code

Comparative path lengths in amino acid biosynthesis and other molecular indicators of the timing of codon assignment were examined to reconstruct the main stages of code evolution. The codon tree obtained was rooted in the 4 N-fixing amino acids (Asp, Glu, Asn, Gln) and 16 triplets of the NAN set. This small, locally phased (commaless) code evidently arose from ambiguous translation on a poly(A) collector strand, in a surface reaction network. Copolymerisation of these amino acids yields polyanionic peptide chains, which could anchor uncharged amide residues to a positively charged mineral surface. From RNA virus structure and replication in vitro, the first genes seemed to be RNA segments spliced into tRNA. Expansion of the code reduced the risk of mutation to an unreadable codon. This step was conditional on initiation at the 5'-codon of a translated sequence. Incorporation of increasingly hydrophobic amino acids accompanied expansion. As codons of the NUN set were assigned most slowly, they received the most nonpolar amino acids. The origin of ferredoxin and Gln synthetase was traced to mid-expansion phase. Surface metabolism ceased by the end of code expansion, as cells bounded by a proteo-phospholipid membrane, with a protoATPase, had emerged. Incorporation of positively charged and aromatic amino acids followed. They entered the post-expansion code by codon capture. Synthesis of efficient enzymes with acid-base catalysis was then possible. Both types of aminoacyl-tRNA synthetases were attributed to this stage. tRNA sequence diversity and error rates in RNA replication indicate the code evolved within 20 million yr in the preIsuan era. These findings on the genetic code provide empirical evidence, from a contemporaneous source, that a surface reaction network, centred on C-fixing autocatalytic cycles, rapidly led to cellular life on Earth.

The beta-sheets of proteins, the biosynthetic relationships between amino acids, and the origin of the genetic code

Two forces are generally hypothesised as being responsible for conditioning the origin of the organization of the genetic code: the physicochemical properties of amino acids and their biosynthetic relationships (relationships between precursor and product amino acids). If we assume that the biosynthetic relationships between amino acids were fundamental in defining the genetic code, then it is reasonable to expect that the distribution of physicochemical properties among the amino acids in precursor-product relationships cannot be random but must, rather, be affected by some selective constraints imposed by the structure of primitive proteins. Analysis shows that measurements representing the 'size' of amino acids, e.g. bulkiness, are specifically associated to the pairs of amino acids in precurso-product relationships. However, the size of amino acids cannot have been selected per se but, rather, because it reflects the beta-sheets of proteins which are, therefore, identified as the main adaptive theme promoting the origin of genetic code organization. Whereas there are no traces of the alpha-helix in the genetic code table. The above considerations make it necessary to re-examine the relationship linking the hydrophilicity of the dinucleoside monophosphates of anticodons and the polarity and bulkiness of amino acids. It can be concluded that this relationship seems to be meaningful only between the hydrophilicity of anticodons and the polarity of amino acids. The latter relationship is supposed to have been operative on hairpin structures, ancestors of the tRNA molecule. Moreover, it is on these very structures that the biosynthetic links between precursor and product amino acids might have been achieved, and the interaction between the hydrophilicity of anticodons and the polarity of amino acids might have had a role in the concession of codons (anticodons) from precursors to products.

A quantitative measure of error minimization in the genetic code

We have calculated the average effect of changing a codon by a single base for all possible single-base changes in the genetic code and for changes in the first, second, and third codon positions separately. Such values were calculated for an amino acid's polar requirement, hydropathy, molecular volume, and isoelectric point. For each attribute the average effect of single-base changes was also calculated for a large number of randomly generated codes that retained the same level of redundancy as the natural code. Amino acids whose codons differed by a single base in the first and third codon positions were very similar with respect to polar requirement and hydropathy. The major differences between amino acids were specified by the second codon position. Codons with U in the second position are hydrophobic, whereas most codons with A in the second position are hydrophilic. This accounts for the observation of complementary hydropathy. Single-base changes in the natural code had a smaller average effect on polar requirement than all but 0.02% of random codes. This result is most easily explained by selection to minimize deleterious effects of translation errors during the early evolution of the code.

Antisense peptides: tools for receptor isolation? Lack of antisense MSH and ACTH to interact with their sense peptides and to induce receptor-specific antibodies

The use of antisense peptides for receptor isolation as proposed by Blalock and his colleagues (e.g. TIBTECH 8, 140-144, 1990) was tested for human ACTH as well as alpha- and beta-MSH. We synthesized the corresponding antisense peptides HTCAh, HSM-alpha and HSM-beta and analyzed them for specific interaction with the sense peptides using several types of binding assay and bioassay. Similarly HTCAh antibodies were tested for binding to ACTH receptors and ACTH antibodies. All these experiments were negative, i.e. there was no specific interaction between sense and antisense peptides nor between the corresponding antibodies. Receptor binding of the sense peptides was not affected by the antisense peptides or HTCAh antibodies. Unexpectedly, HTCAh but not HSM-alpha or HSM-beta was a weak MSH agonist acting through a site independent of the MSH receptor. A detailed analysis of the concept of antisense peptides revealed that the theoretical background of the hypothesis of the 'molecular recognition theory' is rather weak, explaining the failure of various attempts to obtain specific receptor antibodies.

An estimate on the effect of point mutation and natural selection on the rate of amino acid replacement in proteins

We outline a method for estimating quantitatively the influence of point mutations and selection on the frequencies of codons and amino acids. We show how the mutation rate, i.e., the rate of amino acid replacement due to point mutation, can be affected by the codon usage as well as by the rates of the involved base exchanges. A comparison of the mutation rates calculated from reliable values of codon usage and base exchange probabilities with those that would be expected on the basis of chance reveals a notable suppression of replacements leading to tryptophan, glutamate, lysine, and methionine, and particularly of those leading to the termination codons. If selection constraints are neglected and only mutations are taken into account, the best agreement between expected and observed frequencies of both codons and amino acids is obtained for alpha = 1.13-1.15, where (Formula: see text). The "selection values" of codons and amino acids derived by our method show a pattern that partially deviates from others in the literature. For example, the selection pressure on methionine and cysteine turns out to be much more pronounced than expected if only the discrepancies between their observed and expected occurrences in proteins are considered. To estimate to what extent randomly occurring amino acid replacements are accepted by selection, we constructed an "acceptability matrix" from the well-established matrix of accepted point mutations. On the basis of this matrix "acceptability values" of the amino acids can be defined that correlate with their selection values. We also examine the significance of mutations and selection of amino acids with respect to their physicochemical properties and functions in proteins. The conservatism of amino acid replacements with respect to certain properties such as polarity can be brought about by the mutational process alone, whereas the conservatism with respect to other relevant properties--among them all measures of bulkiness--obviously is the result of additional selectional constraints on the evolution of protein structures.

First approximation of a stereochemical rationale for the genetic code based on the topography and physicochemical properties of "cavities" constructed from models of DNA

To examine the question of whether or not the genetic code has a stereochemical basis, we used artificial constructs of the topography and physicochemical features of unique "cavities" formed by removal of the second codon base in B-DNA. The effects of base changes on the stereochemistry of the cavities are consistent with the pattern of the genetic code. Fits into the cavities of the side chains of the 20 L amino acids involved in protein synthesis can be demonstrated by using conventional physicochemical principles of hydrogen bonding and steric constraints. The specificity of the fits is remarkably consistent with the genetic code.

Are there structural analogies between amino acids and nucleic acids?

Space-filling molecular models have been used to examine structural analogies between amino acids and nucleic acids. The three-dimensional structures of amino acid R groups appear to be stereochemically related to cavities formed by removal of single bases in double helical nucleic acids. The common L amino acids may thus be complementary to their codons.

Role of minimization of chemical distances between amino acids in the evolution of the genetic code

The allocation of codons in the genetic code makes possible a moderate minimization of the chemical distances between pairs of neighboring amino acids in the code. However, the code is neither a global nor a local optimum with respect to distance minimization. These findings do not support the physicochemical postulate that distance minimization was a major factor shaping the evolution of the genetic code. They agree with the coevolution theory, which proposes that genetic code evolution was predominantly determined by the concession of codons from precursor to product amino acids in an expansion of the code to accommodate new varieties of amino acids, with distance minimization playing a subsidiary role in deciding the choice of codons to be acquired by the product amino acids from the codon domains of the precursor amino acids.

Inadequacy of prebiotic synthesis as origin of proteinous amino acids

The production of some nonproteinous, and lack of production of other proteinous, amino acids in model prebiotic synthesis, along with the instability of glutamine and asparagine, suggest that not all of the 20 present day proteinous amino acids gained entry into proteins directly from the primordial soup. Instead, a process of active co-evolution of the genetic code and its constituent amino acids would have to precede the final selection of these proteinous amono acids.