Backbones of folded proteins reveal novel invariant amino acid neighborhoods

Folding of naturally occurring proteins has eluded a universal molecular level explanation till date. Rather, there is an abundance of diverse views on dominant factors governing protein folding. Through rigorous analyses of several thousand crystal structures, we observe that backbones of folded proteins display some remarkable invariant features. Folded proteins are characterized by spatially well-defined, distance dependent, and universal, neighborhoods of amino acids which defy any of the conventionally prevalent views. These findings present a compelling case for a newer view of protein folding which takes into account solvent mediated and amino acid shape and size assisted optimization of the tertiary structure of the polypeptide chain to make a functional protein.     

Backbones of Folded Proteins Reveal Novel Invariant Amino Acid Neighborhoods

Abstract

Folding of naturally occurring proteins has eluded a universal molecular level explanation till date. Rather, there is an abundance of diverse views on dominant factors governing protein folding. Through rigorous analyses of several thousand crystal structures, we observe that backbones of folded proteins display some remarkable invariant features. Folded proteins are characterized by spatially well-defined, distance dependent, and universal, neighborhoods of amino acids which defy any of the conventionally prevalent views. These findings present a compelling case for a newer view of protein folding which takes into account solvent mediated and amino acid shape and size assisted optimization of the tertiary structure of the polypeptide chain to make a functional protein.     

A Model of Proto-Anti-Codon RNA Enzymes Requiring l-Amino Acid Homochirality

All living organisms encode the 20 natural amino acid units of polypeptides using a universal scheme of triplet nucleotide "codons". Disparate features of this codon scheme are potentially informative of early molecular evolution: (i) the absence of any codons for D-amino acids; (ii) the odd combination of alternate codon patterns for some amino acids; (iii) the confinement of synonymous positions to a codon's third nucleotide; (iv) the use of 20 specific amino acids rather than a number closer to the full coding potential of 64; and (v) the evolutionary relationship of patterns in stop codons to amino acid codons. Here I propose a model for an ancestral proto-anti-codon RNA (pacRNA) auto-aminoacylation system and show that pacRNAs would naturally manifest features of the codon table. I show that pacRNAs could implement all the steps for auto-aminoacylation: amino acid coordination, intermediate activation of the amino acid by the 5'-end of the pacRNA, and 3'-aminoacylation of the pacRNA. The anti-codon cradles of pacRNAs would have been able to recognize and coordinate only a small number of L-amino acids via hydrogen bonding. A need for proper spatial coordination would have limited the number of chargeable amino acids for all anti-codon sequences, in addition to making some anti-codon sequences unsuitable. Thus, the pacRNA model implies that the idiosyncrasies of the anti-codon table and L-amino acid homochirality co-evolved during a single evolutionary period. These results further imply that early life consisted of an aminoacylated RNA world with a richer enzymatic potential than ribonucleotides alone.     

Concentration and Extinction Coefficient Determination for Oligonucleotides and Analogs Using a General Phosphate Analysis

A general phosphate analysis (GPA) is developed which assays the concentration of nucleic acid oligomers and their analogs based on stoichiometric phosphorus in the sequence. The method involves complete digestion of the oligomer sample to orthophosphate using acid at high temperature and subsequent colorimetric analysis by phosphomolybdate complex formation. GPA is applied to oligomers having phosphodiester, methylphosphonate, and phosphorothioate backbone linkages. Given the absorption spectra of oligomers having these backbones, extinction coefficients are obtained and compared to other quantitative and predictive methods. In addition to sequences having the usual nucleoside residues found in naturally occurring nucleic acids, oligomers having base analog residues can be readily quantified by GPA.     

A large scale test of computational protein design: folding and stability of nine completely redesigned globular proteins

A previously developed computer program for protein design, RosettaDesign, was used to predict low free energy sequences for nine naturally occurring protein backbones. RosettaDesign had no knowledge of the naturally occurring sequences and on average 65% of the residues in the designed sequences differ from wild-type. Synthetic genes for ten completely redesigned proteins were generated, and the proteins were expressed, purified, and then characterized using circular dichroism, chemical and temperature denaturation and NMR experiments. Although high-resolution structures have not yet been determined, eight of these proteins appear to be folded and their circular dichroism spectra are similar to those of their wild-type counterparts. Six of the proteins have stabilities equal to or up to 7kcal/mol greater than their wild-type counterparts, and four of the proteins have NMR spectra consistent with a well-packed, rigid structure. These encouraging results indicate that the computational protein design methods can, with significant reliability, identify amino acid sequences compatible with a target protein backbone.     

In vitro selection for sense codon suppression

The universal genetic code links the 20 naturally occurring amino acids to the 61 sense codons. Previously, the UAG amber stop codon (a nonsense codon) has been used as a blank in the code to insert natural and unnatural amino acids via nonsense suppression. We have developed a selection methodology to investigate whether the unnatural amino acid biocytin could be incorporated into an mRNA display library at sense codons. In these experiments we probed a single randomized NNN codon with a library of 16 orthogonal, biocytin-acylated tRNAs. In vitro selection for efficient incorporation of the unnatural amino acid resulted in templates containing the GUA codon at the randomized position. This sense suppression occurs via Watson-Crick pairing with similar efficiency to UAG-mediated nonsense suppression. These experiments suggest that sense codon suppression is a viable means to expand the chemical and functional diversity of the genetic code.     

ENVIRON: a software package to compare protein three-dimensional structures with homologous sequences using local structural motifs

This work presents a method to compare local clusters of interactingresidues as observed in a known three-dimensional protein structurewith corresponding clusters inferred from homologous proteinsequences, assuming conserved protein folding. For this purposethe local environment of a selected residue in a known proteinstructure is defined as the ensemble of amino acids in contactwith it in the folded state. Using a multiple sequence alignmentto identify corresponding residues in homologous proteins, adetailed comparison can be performed between the local environmentof a selected amino acid in the template protein structure andthe expected local environments at the sets of equivalent residues,derived from the aligned protein sequences. The comparison makesit possible to detect conserved local features such as hydrogenbonding or complementarity in residue substitution. A globalmeasure of environmental similarity is also defined, to searchfor conserved amino acid clusters subject to functional or structural constraints. The proposed approach is useful for investigatingprotein function as well as for site-directed mutagenesis experiments,where appropriate amino acid substitutions can be suggestedby observing naturally occurring protein variants.     

Systematic variation of amino acid substitutions for stringent assessment of pairwise covariation

During protein evolution, amino acids change due to a combination of functional constraints and genetic drift. Proteins frequently contain pairs of amino acids that appear to change together (covariation). Analysis of covariation from naturally occurring sets of orthologs cannot distinguish between residue pairs retained by functional requirements of the protein and those pairs existing due to changes along a common evolutionary path. Here, we have separated the two types of covariation by independently recombining every naturally occurring amino acid variant within a set of 15 subtilisin orthologs. Our analysis shows that in this family of subtilisin orthologs, almost all possible pairwise combinations of amino acids can coexist. This suggests that amino acid covariation found in the subtilisin orthologs is almost entirely due to common ancestral origin of the changes rather than functional constraints. We conclude that naturally occurring sequence diversity can be used to identify positions that can vary independently without destroying protein function.     

Reduction of protein sequence complexity by residue grouping

It is well known that there are some similarities among various naturally occurring amino acids. Thus, the complexity in protein systems could be reduced by sorting these amino acids with similarities into groups and then protein sequences can be simplified by reduced alphabets. This paper discusses how to group similar amino acids and whether there is a minimal amino acid alphabet by which proteins can be folded. Various reduced alphabets are obtained by reserving the maximal information for the simplified protein sequence compared with the parent sequence using global sequence alignment. With these reduced alphabets and simplified similarity matrices, we achieve recognition of the protein fold based on the similarity score of the sequence alignment. The coverage in dataset SCOP40 for various levels of reduction on the amino acid types is obtained, which is the number of homologous pairs detected by program BLAST to the number marked by SCOP40. For the reduced alphabets containing 10 types of amino acids, the ability to detect distantly related folds remains almost at the same level as that by the alphabet of 20 types of amino acids, which implies that 10 types of amino acids may be the degree of freedom for characterizing the complexity in proteins.     

Comparative characterization of random‐sequence proteins consisting of 5, 12, and 20 kinds of amino acids

Screening of functional proteins from a random‐sequence library has been used to evolve novel proteins in the field of evolutionary protein engineering. However, random‐sequence proteins consisting of the 20 natural amino acids tend to aggregate, and the occurrence rate of functional proteins in a random‐sequence library is low. From the viewpoint of the origin of life, it has been proposed that primordial proteins consisted of a limited set of amino acids that could have been abundantly formed early during chemical evolution. We have previously found that members of a random‐sequence protein library constructed with five primitive amino acids show high solubility (Doi et al., Protein Eng Des Sel 2005;18:279–284). Although such a library is expected to be appropriate for finding functional proteins, the functionality may be limited, because they have no positively charged amino acid. Here, we constructed three libraries of 120‐amino acid, random‐sequence proteins using alphabets of 5, 12, and 20 amino acids by preselection using mRNA display (to eliminate sequences containing stop codons and frameshifts) and characterized and compared the structural properties of random‐sequence proteins arbitrarily chosen from these libraries. We found that random‐sequence proteins constructed with the 12‐member alphabet (including five primitive amino acids and positively charged amino acids) have higher solubility than those constructed with the 20‐member alphabet, though other biophysical properties are very similar in the two libraries. Thus, a library of moderate complexity constructed from 12 amino acids may be a more appropriate resource for functional screening than one constructed from 20 amino acids.     

The effects of bioregulators upon amino acid transport and protein synthesis in isolated rat hepatocytes

Isolated rat hepatocytes prepared by an enzyme perfusion technique possess a functional amino acid transport system and retain the capacity to synthesize protein. Amino acid transport was studied using the non-metabolizable amino acid analog alpha-aminoisobutyric acid. The transport process was time, temperature and concentration dependent. Similarly, leucine incorporation into protein was time and temperature dependent being optimal at 3m degrees C. Amino acid, fetal calf serum, growth hormone and glucose all produced small, reproducible increases in protein synthesis rates. Bovine serum albumin diminished the uptake of alpha-aminoisobutyric acid and leucine incorporation into protein. The amino acid content on either side of the cell membrane was found to affect transport into or out of the cellular compartment (transconcentration effects). High cell concentrations decreased transport and protein synthesis as a result of isotopic dilution of labelled amino acids with those released by the hepatocytes. This was consistent with the capacity of naturally occurring amino aicds to compete with alpha-aminoisobutyric acid for uptake into the hepatocyte. In order to define more precisely the effects of bioregulators on transport and protein synthesis it will be necessary to define and subfractionate cellular compartments and proteins which are the specific targets of cellular regulation.     

Evolution of a pathway to novel long-chain carotenoids

Using methods of laboratory evolution to force the C(30) carotenoid synthase CrtM to function as a C(40) synthase, followed by further mutagenesis at functionally important amino acid residues, we have discovered that synthase specificity is controlled at the second (rearrangement) step of the two-step reaction. We used this information to engineer CrtM variants that can synthesize previously unknown C(45) and C(50) carotenoid backbones (mono- and diisopentenylphytoenes) from the appropriate isoprenyldiphosphate precursors. With this ability to produce new backbones in Escherichia coli comes the potential to generate whole series of novel carotenoids by using carotenoid-modifying enzymes, including desaturases, cyclases, hydroxylases, and dioxygenases, from naturally occurring pathways.     

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.     

The triplet genetic code had a doublet predecessor

Information theoretic analysis of genetic languages indicates that the naturally occurring 20 amino acids and the triplet genetic code arose by duplication of 10 amino acids of class-II and a doublet genetic code having codons NNY and anticodons GNN. Evidence for this scenario is presented based on the properties of aminoacyl-tRNA synthetases, amino acids and nucleotide bases.     

Recapitulation of protein family divergence using flexible backbone protein design

We use flexible backbone protein design to explore the sequence and structure neighborhoods of naturally occurring proteins. The method samples sequence and structure space in the vicinity of a known sequence and structure by alternately optimizing the sequence for a fixed protein backbone using rotamer based sequence search, and optimizing the backbone for a fixed amino acid sequence using atomic-resolution structure prediction. We find that such a flexible backbone design method better recapitulates protein family sequence variation than sequence optimization on fixed backbones or randomly perturbed backbone ensembles for ten diverse protein structures. For the SH3 domain, the backbone structure variation in the family is also better recapitulated than in randomly perturbed backbones. The potential application of this method as a model of protein family evolution is highlighted by a concerted transition to the amino acid sequence in the structural core of one SH3 domain starting from the backbone coordinates of an homologous structure.     

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.     

Removing the redundancy from randomised gene libraries

Amino acid substitution plays a vital role in both the molecular engineering of proteins and analysis of structure-activity relationships. High-throughput substitution is achieved by codon randomisation, which generates a library of mutants (a randomised gene library) in a single experiment. For full randomisation, key codons are typically replaced with NNN (64 sequences) or NN(G)(CorT) (32 sequences). This obligates cloning of redundant codons alongside those required to encode the 20 amino acids. As the number of randomised codons increases, there is therefore a progressive loss of randomisation efficiency; the number of genes required per protein rises exponentially. The redundant codons cause amino acids to be represented unevenly; for example, methionine is encoded just once within NNN, whilst arginine is encoded six times. Finally, the organisation of the genetic code makes it impossible to encode functional subsets of amino acids (e.g. polar residues only) in a single experiment. Here, we present a novel solution to randomisation where genetic redundancy is eliminated; the number of different genes equals the number of encoded proteins, regardless of codon number. There is no inherent amino acid bias and any required subset of amino acids may be encoded in one experiment. This generic approach should be widely applicable in studies involving randomisation of proteins.     

Amino acid pairing

A set of amino acid pairings are presented which may allow protein-to-protein information transfers. Amino acid pairing is only possible on a parallel β ribbon and involves both the polypeptide backbones and the side chains. Model building revealed that of the 210 possible amino acid pairs of the standard 20 amino acids, no more than 26 could be built to meet standard criteria for bonding. Of these 26, 14 were found to be genetically encoded when the codons are read as if they paired in a parallel manner (i.e. in a manner reflecting the structural parallelism of the amino acid pairings); the other 12 pairings were derivatives of the coded pairings in which a single base of the codon triplet had been varied in accordance with Crick's (1966) “wobble hypothesis.” Evidence for the pairings is presented from colligative studies of polyamino acids. Ways of testing the hypothesis further are suggested. Its implications for the Central Dogma and theories of the origin of life are discussed.     

What limits the primary sequence space of natural proteins?

Abstract

Number of naturally occurring primary sequences of proteins is an infinitesimally small subset of the possible number of primary sequences that can be synthesized using 20 amino acids. Prevailing views ascribe this to slow and incremental mutational/selection evolutionary mechanisms. However, considering the large number of avenues available in form of diversity of emerging/evolving and/or disappearing living systems for exploring the primary sequence space over the evolutionary time scale of ～3.5 billion years, this remains a conjecture. Therefore, to investigate primary sequence space limitations, we carried out a systematic study for finding primary sequences absent in nature. We report the discovery of the smallest peptide sequence “Cysteine-Glutamine-Tryptophan-Tryptophan” that is not found in over half-a-million curated protein sequences in the Uniprot (Swiss-Prot) database. Additionally, we report a library of 83605 pentapeptides that are not found in any of the known protein sequences. Compositional analyses of these absent primary sequences yield a remarkably strong power relationship between the percentage occurrence of individual amino acids in all known protein sequences and their respective frequency of occurrence in the absent peptides, regardless of their specific position in the sequences. If random evolutionary mechanisms were responsible for limitations to the primary sequence space, then one would not expect any relationship between compositions of available and absent primary sequences. Thus, we conclusively show that stoichiometric constraints on amino acids limit the primary sequence space of proteins in nature. We discuss the possibly profound implications of our findings in both evolutionary and synthetic biology.

Communicated by Ramaswamy H. Sarma     

Regulation of aromatic amino acid biosynthesis in phenazine-producing strains of Pseudomonas

Regulation of 3-deoxy-d-arabino-heptulosonate-7-phosphate (DAHP) synthetase was studied in eight strains of Pseudomonas which synthesize phenazine compounds. Repression studies with individual aromatic amino acids led to the finding that enzyme synthesis was repressed in only one strain, P. aureofaciens B1543p, and by only one amino acid, l-tyrosine. Feedback inhibition by the aromatic amino acids varied from strain to strain in terms of the type of inhibitory control, and the particular acid or acids which inhibited. Prephenate and chorismate, as well as a number of naturally occurring phenazine compounds, inhibited the DAHP synthetase activity to varying degrees.