Code of codon roots, determining the intra- and intermolecular interaction of amino acids in peptide chains]

To study the recognition processes and interaction of peptides and proteins, a model has been suggested according to which the first steps of complex formation of molecules are defined by G/C and A/U complementarity of codon roots of amino acid forming the molecules contact sites or surfaces. In contrast to amino acid--antiamino acid interaction code (L. B. Mekler, 1969), the code of codon roots involves the interaction of amino acids independently of the base structures in nucleotide triplets in positions 1 and 3. The analysis of the spectra of point mutations homologous proteins confirms the possible role of the root code.     

Hidden symmetry of peptide and protein primary structures]

The internal symmetry of peptide chains was considered. To identify symmetrically located equivalent amino acids, the signatures method and the code of amino acid codon roots were applied. There was revealed the hidden symmetry of amino acid sequences of peptides and proteins as well as of their active centres. Amino acids having common codon roots in primary (and supposedly in the spatial "biologically active") molecular structures, are located symmetrically. Definition of local symmetry of peptide chains was proposed to use as one of the elements of complex analysis to determine location of molecular active centres.     

Hidden symmetry of the genetic code and laws of amino acid interaction]

Natural amino acids having common antiamino acids are divided into families and groups according to the algorithm of the genetic code (a-n-n-a, amino acid-codon-anticodon-antiamino acid). Members of these groups are placed symmetrically in the structure of the genetic code. In the course of evolution, those point mutations are predominantly accepted retained. In homologous proteins of phylogenetically related organisms which lend to amino acids belonging to one family or group and having common antiamino acids. This assumption is in agreement with L. B. Mekler's theory (1969) of the amino acid interaction code a-a.     

Interaction code for polar and nonpolar amino acids: "ice-breaker" model]

A novel model for the study of recognition and interaction code of amino acids in peptides, proteins and their complexes has been proposed. The model is designed on the modern notions on the structure and properties of water and hydrophobic bonds. It is assumed that the polar side chains of amino acids during the formation of the hydrophobic bonds act as "ice-breaker", thus destroying the organized structure of water (clusters or "icebergs") around the hydrophobic radicals of amino acids.     

The optimality of the standard genetic code assessed by an eight-objective evolutionary algorithm

Background

The standard genetic code (SGC) is a unique set of rules which assign amino acids to codons. Similar amino acids tend to have similar codons indicating that the code evolved to minimize the costs of amino acid replacements in proteins, caused by mutations or translational errors. However, if such optimization in fact occurred, many different properties of amino acids must have been taken into account during the code evolution. Therefore, this problem can be reformulated as a multi-objective optimization task, in which the selection constraints are represented by measures based on various amino acid properties.

Results

To study the optimality of the SGC we applied a multi-objective evolutionary algorithm and we used the representatives of eight clusters, which grouped over 500 indices describing various physicochemical properties of amino acids. Thanks to that we avoided an arbitrary choice of amino acid features as optimization criteria. As a consequence, we were able to conduct a more general study on the properties of the SGC than the ones presented so far in other papers on this topic. We considered two models of the genetic code, one preserving the characteristic codon blocks structure of the SGC and the other without this restriction. The results revealed that the SGC could be significantly improved in terms of error minimization, hereby it is not fully optimized. Its structure differs significantly from the structure of the codes optimized to minimize the costs of amino acid replacements. On the other hand, using newly defined quality measures that placed the SGC in the global space of theoretical genetic codes, we showed that the SGC is definitely closer to the codes that minimize the costs of amino acids replacements than those maximizing them.

Conclusions

The standard genetic code represents most likely only partially optimized systems, which emerged under the influence of many different factors. Our findings can be useful to researchers involved in modifying the genetic code of the living organisms and designing artificial ones.

     

Genetic code redundancy and its differential influence on the evolution of protein interiors versus exteriors

The distinctive amino acid compositions of protein exteriors and interiors were compared to the composition bias imposed by genetic code redundancy. It transpired that the synonym allocation is biased more in favour of those residues which are preferred in interiors, and this leads to an average interior residue being more probable and less mutable compared to an exterior residue. The general implications for protein evolution are discussed in association with the known evolutionary behavior of particular protein families. It is suggested that some proteins may have their structural history "fossilised" in their interiors and that the "amino acid" code is in reality a "protein" code.     

Linear and inverted repetitions in protein sequences

An extensive search for internal regularities in amino acid sequences has been made, using both the genetic code and the relative frequencies of amino acid alternatives in homologous proteins. The two methods give very similar results and strongly suggest the occurrence of significant linear and inverted repetitions (similar sequences of opposite polarity) in several proteins. A hypothesis is developed to explain the occurrence of such internal regularities in proteins. This hypothesis is based on a process of duplication of an ancestral loop in which a symmetrical arrangement of amino acid allows stabilization by interaction between the amino acid side chains.     

Amino acid substitutions in VH CDR2 change the idiotype but not the antigen-binding of monoclonal antibodies to alpha(1----6)dextrans

An idiotype defined by mAb and polyclonal antibodies to 10.16.1, an anti-alpha(1----6) dextran was previously reported to be expressed on most BALB/c anti-alpha(1----6)dextrans with groove-type sites and to involved CDR3 and probably CDR2. By comparing amino acid sequences of VH and VL derived from cDNA of idiotype+ and idiotype- anti-alpha(1----6)dextran hybridoma proteins, an idiotope was assigned to VH CDR2. Substitution of phenylalanine for leucine at residue 52 in CDR2 coupled with amino acid changes at either residue 58 or residues 57 and 60 abolished expression of this idiotype without affecting Ag binding.     

Regrouping of peptide chains in molecules and molecular complexes: recognition and spatial conjugation codes of amino acids

A two-step model for reactions between peptide and protein molecules in aqueous medium is considered. The first stage of the reaction involves specific recognition and primary complex formation. This process is governed by the amino acid interaction code a-a as a part of genetic code (algorithm a-n-n-a, amino acid-codon-anticodon-anti-amino acid). According to the a-a code, the primary complex formation is determined by amino acid pairs of opposite polarity. During the second stage of the reaction, when the contacting ligand and receptor surfaces undergo dehydration, the primary complex becomes rearranged. The new structure is mainly determined by pairwise contacts of amino acids having similar polarity and belonging to the same amino acid family.     

Folding of peptide chains during formation of the spatial structure of protein molecules is determined by the code of amino acid codon roots]

Basing on the analysis of a large number of protein sequences (Cserzo M., Simon I., 1989), the structure of the amino acid nearest neighbour pair whose occurrence has a maximal positive deviation from the mean statistical value, is shown to correspond in most cases to the code of the amino acid codon roots. It reveals particularly amino acid pairs in n and n+5 positions of polypeptide chains. Amino acids belonging to A/U family contribute mostly to the folding of peptide chains.     

Neutral adaptation of the genetic code to double-strand coding

Summary We lay new foundations to the hypothesis that the genetic code is adapted to evolutionary retention of information in the antisense strands of natural DNA/RNA sequences. In particular, we show that the genetic code exhibits, beyond the neutral replacement patterns of amino acid substitutions, optimal properties by favoring simultaneous evolution of proteins encoded in DNA/RNA sense-antisense strands. This is borne out in the sense-antisense transformations of the codons of every amino acid which target amino acids physicochemically similar to each other. Moreover, silent mutations in the sense strand generate conservative ones in its antisense counterpart and vice versa. Coevolution of proteins coded by complementary strands is shown to be a definite possibility, a result which does not depend on any physical interaction between the coevolving proteins. Likewise, the degree to which the present genetic code is dedicated to evolutionary sense-antisense tolerance is demonstrated by comparison with many randomized codes. Double-strand coding is quantified from an information-theoretical point of view.     

Protein evolution drives the evolution of the genetic code and vice versa

A model for the developmental pathway of the genetic code, grounded on group theory and the thermodynamics of codon-anticodon interaction is presented. At variance with previous models, it takes into account not only the optimization with respect to amino acid attributes but, also physicochemical constraints and initial conditions. A 'simple-first' rule is introduced after ranking the amino acids with respect to two current measures of chemical complexity. It is shown that a primeval code of only seven amino acids is enough to build functional proteins. It is assumed that these proteins drive the further expansion of the code. The proposed primeval code is compared with surrogate codes randomly generated and with another proposal for primeval code found in the literature. The departures from the 'universal' code, observed in many organisms and cellular compartments, fit naturally in the proposed evolutionary scheme. A strong correlation is found between, on one side, the two classes of aminoacyl-tRNA synthetases, and on the other, the amino acids grouped by end-atom-type and by codon type. An inverse of Davydov's rules, to associate the amino acid end atoms (O/N and non-O/non-N) of 18 amino acids with codons containing a weak base (A/U), extended to the 20 amino acids, is derived.     

Natural expansion of the genetic code

At the time of its discovery four decades ago, the genetic code was viewed as the result of a "frozen accident." Our current knowledge of the translation process and of the detailed structure of its components highlights the roles of RNA structure (in mRNA and tRNA), RNA modification (in tRNA), and aminoacyl-tRNA synthetase diversity in the evolution of the genetic code. The diverse assortment of codon reassignments present in subcellular organelles and organisms of distinct lineages has 'thawed' the concept of a universal immutable code; it may not be accidental that out of more than 140 amino acids found in natural proteins, only two (selenocysteine and pyrrolysine) are known to have been added to the standard 20-member amino acid alphabet. The existence of phosphoseryl-tRNA (in the form of tRNACys and tRNASec) may presage the discovery of other cotranslationally inserted modified amino acids.     

Relationships between genomic base content and distribution of mass in coded proteins

The aim of this research was to examine the possible significance of genome/protein relationships in terms of effects on distribution of mass, especially in proteins. Amino acid residues in proteins have side-chains and polypeptide segments. We use "SCM" (side-chain mass), "MCM" (main-chain mass), and "deltaM" (SCM-MCM) as the deviation from "mass balance." Total MCM of the 61 amino acids in the standard code, 3412, equals total SCM: they form a mass balanced set (mean deltaM = 0). Of 14 natural variants of the code, seven have slightly positive mean deltaM values and seven have slightly negative values. Codes with the standard amino acids assigned randomly to the 20 codon sets of the standard code have about one chance in 3,300 of producing a mass balanced set. In natural proteins, as %A + T increases, the proportion of the mass in the side-chains also increases, by about half the amount calculated for standard genes with various AT/GC ratios, partly due to selection of codons with greater variability in composition at synonymous sites. For 203 representative species (including organelles), the total protein mass is distributed approximately equally between SCM and MCM (overall mean deltaM/amino acid residue, -0.06). The attainment of some overall macromolecular mass balance may have been a criterion for selecting the codon/amino acid pairs. When both structural and dynamic requirements are considered, a genetic code based on hydrophobicity and mass balance as key properties seems likely.     

A model of intra- and intermolecular interactions of peptide chains]

Results of attempts to determine the code of interaction of amino acids in peptide chains proceeding from their coding nucleotide sequences have been summarized. According to the model suggested the G/C and A/U complementarity of codon roots determines the mutual binding of coded amino acid residues. Structures of analogs of the immunoactive peptide, a fragment of IgG1 (336-370) EPQVY have been constructed on the basis of the model.     

Complete heavy and light chain variable region sequence of anti-arsonate monoclonal antibodies from BALB/c and A/J mice sharing the 36-60 idiotype are highly homologous

Structural and serologic studies on murine A/J monoclonal anti-arsonate antibodies resulted in the identification of a second idiotype family (Id36-60) in addition to the predominant idiotype family (IdCR). Id36-60, unlike IdCR, is a dominant idiotype in the BALB/c strain but is a "minor" idiotype in the A/J strain. The complete heavy and light chain variable region (VH and VL) amino acid sequences of a representative Id36-60 hybridoma protein from both the A/J and BALB/c strains have been determined. There are only four amino acid sequence differences between the VH of antibody 36-60 (A/J) and antibody 1210.7 (BALB/c). Two of these differences arise from single nucleotide changes in which the A/J and BALB/c Id36-60 VH germline gene sequences differ. The two other differences are the result of somatic mutation in hybridoma protein 36-60. In addition, Id36-60 heavy chains employ the same D and JH3 segments in both strains. The entire Vk2 VL of 36-60 and 1210.7 differ by only two amino acids, suggesting that like the heavy chains, they are derived from highly homologous VL genes. The same Jk segment is used in both antibodies. A comparison of the amino acid sequence data from Id36-60-bearing hybridomas suggests that a heavy chain amino acid difference accounts for the diminished arsonate binding by the 1210.7 hybridoma protein. Because the 1210.7 heavy chain is the unmutated product of the BALB/c VH gene, somatic mutation in VH may be required to enhance Ars affinity in this system.     

Genetic code redundancy and the evolutionary stability of protein secondary structure

The genetic code has an inherent bias towards some amino acids because of the variable number of synonymous codons per amino acid. The extent to which these biases are expressed in protein secondary structure is described through the analysis of the overall amino acid compositions of the alpha-helix, beta-sheet, beta-turn and random coil segments elucidated by X-ray crystallography. Given the concept of neutral mutation in proteins, the allocation of synonyms in the genetic code appears to protect secondary structures from amino acid changes and discourages the appearance of chemically complex residues. The level of protection is similar for each structural form, despite their clear preferences for certain amino acids. The organization of the code is therefore relevant to the preservation of conformation seen in the evolution of many protein families.     

Use of Chou-Fasman amino acid conformational parameters to analyze the organization of the genetic code and to construct protein genealogies

Summary Chou-Fasman parameters, measuring preferences of each amino acid for different conformational regions in proteins, were used to obtain an amino acid difference index of conformational parameter distance (CPD) values. CPD values were found to be significantly lower for amino acid exchanges representing in the genetic code transitions of purines, GA than for exchanges representing either transitions of pyrimidines, CU, or transversions of purines and pyrimidines. Inasmuch as the distribution of CPD values in these non GA exchanges resembles that obtained for amino acid pairs with double or triple base differences in their underlying codons, we conclude that the genetic code was not particularly designed to minimize effects of mutation on protein conformation. That natural selection minimizes these changes, however, was shown by tabulating results obtained by the maximum parsimony method for eight protein genealogies with a total occurrence of 4574 base substitutions. At the beginning position of the codons GA transitions were in very great excess over other base substitutions, and, conversely, CU transitions were deficient. At the middle position of the codons only fast evolving proteins showed an excess of GA transitions, as though selection mainly preserved conformation in these proteins while weeding out mutations affecting chemical properties of functional sites in slow evolving proteins. In both fast and slow evolving proteins the net direction of transitions and transversions was found to be from G beginning codons to non-G beginning codons resulting in more commonly occurring amino acids, especially alanine with its generalized conformational properties, being replaced at suitable sites by amino acids with more specialized conformational and chemical properties. Historical circumstances pertaining to the origin of the genetic code and the nature of primordial proteins could account for such directional changes leading to increases in the functional density of proteins.In order to further explore the course of protein evolution, a modified parsimony algorithm was developed for constructing protein genealogies on the basis of minimum CPD length. The algorithm's ability to judge with finer discrimination that in protein evolution certain pathways of amino acid substitution should occur more readily than others was considered a potential advantage over strict maximum parsimony. In developing this CPD algorithm, the path of minimum CPD length through intermediate amino acids allowed by the genetic code for each pair of amino acids was determined. It was found that amino acid exchanges representing two base changes have a considerably lower average CPD value per base substitution than the amino acid exchanges representing single base changes. Amino acid exchanges representing three base changes have yet a further marked reduction in CPD per base change. This shows how extreme constraining effects of stabilizing selection can be circumvented, for by way of intermediate amino acids almost any amino acid can ultimately be substituted for another without damage to an evolving protein's conformation during the process.     

Structure-based prediction of DNA target sites by regulatory proteins

Regulatory proteins play a critical role in controlling complex spatial and temporal patterns of gene expression in higher organism, by recognizing multiple DNA sequences and regulating multiple target genes. Increasing amounts of structural data on the protein-DNA complex provides clues for the mechanism of target recognition by regulatory proteins. The analyses of the propensities of base-amino acid interactions observed in those structural data show that there is no one-to-one correspondence in the interaction, but clear preferences exist. On the other hand, the analysis of spatial distribution of amino acids around bases shows that even those amino acids with strong base preference such as Arg with G are distributed in a wide space around bases. Thus, amino acids with many different geometries can form a similar type of interaction with bases. The redundancy and structural flexibility in the interaction suggest that there are no simple rules in the sequence recognition, and its prediction is not straightforward. However, the spatial distributions of amino acids around bases indicate a possibility that the structural data can be used to derive empirical interaction potentials between amino acids and bases. Such information extracted from structural databases has been successfully used to predict amino acid sequences that fold into particular protein structures. We surmised that the structures of protein-DNA complexes could be used to predict DNA target sites for regulatory proteins, because determining DNA sequences that bind to a particular protein structure should be similar to finding amino acid sequences that fold into a particular structure. Here we demonstrate that the structural data can be used to predict DNA target sequences for regulatory proteins. Pairwise potentials that determine the interaction between bases and amino acids were empirically derived from the structural data. These potentials were then used to examine the compatibility between DNA sequences and the protein-DNA complex structure in a combinatorial "threading" procedure. We applied this strategy to the structures of protein-DNA complexes to predict DNA binding sites recognized by regulatory proteins. To test the applicability of this method in target-site prediction, we examined the effects of cognate and noncognate binding, cooperative binding, and DNA deformation on the binding specificity, and predicted binding sites in real promoters and compared with experimental data. These results show that target binding sites for several regulatory proteins are successfully predicted, and our data suggest that this method can serve as a powerful tool for predicting multiple target sites and target genes for regulatory proteins.