To Investigate Protein Evolution by Detecting Suppressed Epitope Structures

Material remains of ancestor nucleotides and proteins are largely unavailable, thus sequence comparison among homologous genes in present-day organisms forms the core of current knowledge of molecular evolution. Variation in protein three-dimensional structure is a basis for functional diversity. To study the evolution of three-dimensional structures in related proteins would significantly improve our understanding of protein evolution and function. A protein may contain ancestor conformations that have been allosterically suppressed by evolutionarily additive structures. Using monoclonal antibody probes to detect such conformation in proteins after removing the suppressor structure, our study demonstrates three-dimensional structure evidence for the evolutionary relationship between troponin I and troponin T, two subunits of the troponin complex in the Ca²⁺-regulatory system of striated muscle, and among their muscle type-specific isoforms. The experimental data show the feasibility of detecting evolutionarily suppressed history-telling structural states in proteins by removing conformational modulator segments added during evolution. In addition to identifying structural modifications that were critical to the emergence of diverged proteins, investigating this novel mode of evolution will help us to understand the origin and functional potential of protein structures.     

Solution structure of Archaeglobus fulgidis peptidyl-tRNA hydrolase (Pth2) provides evidence for an extensive conserved family of Pth2 enzymes in archea, bacteria, and eukaryotes

The solution structure of protein AF2095 from the thermophilic archaea Archaeglobus fulgidis, a 123-residue (13.6-kDa) protein, has been determined by NMR methods. The structure of AF2095 is comprised of four alpha-helices and a mixed beta-sheet consisting of four parallel and anti-parallel beta-strands, where the alpha-helices sandwich the beta-sheet. Sequence and structural comparison of AF2095 with proteins from Homo sapiens, Methanocaldococcus jannaschii, and Sulfolobus solfataricus reveals that AF2095 is a peptidyl-tRNA hydrolase (Pth2). This structural comparison also identifies putative catalytic residues and a tRNA interaction region for AF2095. The structure of AF2095 is also similar to the structure of protein TA0108 from archaea Thermoplasma acidophilum, which is deposited in the Protein Data Bank but not functionally annotated. The NMR structure of AF2095 has been further leveraged to obtain good-quality structural models for 55 other proteins. Although earlier studies have proposed that the Pth2 protein family is restricted to archeal and eukaryotic organisms, the similarity of the AF2095 structure to human Pth2, the conservation of key active-site residues, and the good quality of the resulting homology models demonstrate a large family of homologous Pth2 proteins that are conserved in eukaryotic, archaeal, and bacterial organisms, providing novel insights in the evolution of the Pth and Pth2 enzyme families.     

Emergence of polyproline II-like structure at early stages of experimental evolution from random polypeptides

To examine whether a primordial functional protein at the early stages of evolution has structural features, we carried out experimental evolution consisting of 25 cycles (generations) of mutation and selection toward DNA-binding function using a random-sequence polypeptide of 139 amino acid residues with no secondary structure as the initial sequence. In each generation, 16 clones were sampled arbitrarily for sequence analysis, and a phylogenetic tree was constructed. Polypeptide evolution proceeded from the initial point on branch I in 2 main directions of branches II and III. The initial and 2 evolved polypeptides (one at the 24th generation on branch III and the other at the 25th generation on branch II) were purified to examine their functional and structural properties. Although binding of the initial polypeptide to the target DNA was not detected by surface plasmon resonance measurements, the 2 evolved polypeptides bound to the DNA with dissociation constants of 1.6 and 1.0 microM, respectively, indicating an increase in affinity during the experimental evolution. Circular dichroism spectra of the evolved polypeptides, but not of the initial polypeptide, showed features characteristic of the polyproline II (PPII)-like structure, a left-handed helical structure commonly found in natural proteins, suggesting that the structure emerged through the experimental evolution. The same structural feature was found in another experimental evolution toward catalytic activity. These results demonstrate that the PPII-like structure is one of the common features that could have appeared in the early evolutionary stages of primordial functional protein.     

The role of protein conformational fluctuations in allostery, function, and evolution

It is now well-known that proteins exist at equilibrium as ensembles of conformational states rather than as unique static structures. Here we review from an ensemble perspective important biological effects of such spontaneous fluctuations on protein allostery, function, and evolution. However, rather than present a thorough literature review on each subject, we focus instead on connecting these phenomena through the ensemble-based experimental, theoretical, and computational investigations from our laboratory over the past decade. Special emphasis is given to insights that run counter to some of the prevailing ideas that have emerged over the past 40 years of structural biology research. For instance, when proteins are viewed as conformational ensembles rather than as single structures, the commonly held notion of an allosteric pathway as an obligate series of individual structural distortions loses its meaning. Instead, allostery can result from energetic linkage between distal sites as one Boltzmann distribution of states transitions to another. Additionally, the emerging principles from this ensemble view of proteins have proven surprisingly useful in describing the role of intrinsic disorder in inter-domain communication, functional adaptation mediated by mutational control of fluctuations, and evolutionary conservation of the energetics of protein stability.     

Origins of globular structure in proteins

Since natural proteins are the products of a long evolutionary process, the structural properties of present-day proteins should depend not only on physico-chemical constraints, but also on evolutionary constraints. Here we propose a model for protein evolution, in which membranes play a key role as a scaffold for supporting the gradual evolution from flexible polypeptides to well-folded proteins. We suggest that the folding process of present-day globular proteins is a relic of this putative evolutionary process. To test the hypothesis that membranes once acted as a cradle for the folding of globular proteins, extensive research on membrane proteins and the interactions of globular proteins with membranes will be required.     

Evolutionary computer programming of protein folding and structure predictions

In order to understand the mechanism of protein folding and to assist the rational de-novo design of fast-folding, non-aggregating and stable artificial enzymes it is very helpful to be able to simulate protein folding reactions and to predict the structures of proteins and other biomacromolecules. Here, we use a method of computer programming called "evolutionary computer programming" in which a program evolves depending on the evolutionary pressure exerted on the program. In the case of the presented application of this method on a computer program for folding simulations, the evolutionary pressure exerted was towards faster finding deep minima in the energy landscape of protein folding. Already after 20 evolution steps, the evolved program was able to find deep minima in the energy landscape more than 10 times faster than the original program prior to the evolution process.     

Coronavirus escape from heptad repeat 2 (HR2)-derived peptide entry inhibition as a result of mutations in the HR1 domain of the spike fusion protein

Peptides based on heptad repeat (HR) domains of class I viral fusion proteins are considered promising antiviral drugs targeting virus cell entry. We have analyzed the evolution of the mouse hepatitis coronavirus during multiple passaging in the presence of an HR2-based fusion inhibitor. Drug-resistant variants emerged as a result of multiple substitutions in the spike fusion protein, notably within a 19-residue segment of the HR1 region. Strikingly, one mutation, an A1006V substitution, which consistently appeared first in four independently passaged viruses, was the main determinant of the resistance phenotype, suggesting that only limited options exist for escape from the inhibitory effect of the HR2 peptide.     

Are Proposed Early Genetic Codes Capable of Encoding Viable Proteins?

Proteins are elaborate biopolymers balancing between contradicting intrinsic propensities to fold, aggregate, or remain disordered. Assessing their primary structural preferences observable without evolutionary optimization has been reinforced by the recent identification of de novo proteins that have emerged from previously non-coding sequences. In this paper we investigate structural preferences of hypothetical proteins translated from random DNA segments using the standard genetic code and three of its proposed evolutionarily predecessor models encoding 10, 6, and 4 amino acids, respectively. Our only main assumption is that the disorder, aggregation, and transmembrane helix predictions used are able to reflect the differences in the trends of the protein sets investigated. We found that the 10-residue code encodes proteins that resemble modern proteins in their predicted structural properties. All of the investigated early genetic codes give rise to proteins with enhanced disorder and diminished aggregation propensities. Our results suggest that an ancestral genetic code similar to the proposed 10-residue one is capable of encoding functionally diverse proteins but these might have existed under conditions different from today’s common physiological ones. The existence of a protein functional repertoire for the investigated earlier stages which is quite distinct as it is today can be deduced from the presented results.     

Internal gene duplication in the evolution of prokaryotic transmembrane proteins

We investigated the evolution of transmembrane (TM) topology by detecting partial sequence repeats in TM protein sequences and analyzing them in detail. A total of 377 sequences that seem to have evolved by internal gene duplication events were found among 38,124 predicted TM protein sequences (except for single-spannings) from 87 prokaryotic genomes. Various types of internal duplication patterns were identified in these sequences. The majority of them are diploid-type (including quasi-diploid-type) duplication in which a primordial protein sequence was duplicated internally to become an extant TM protein with twice as many TM segments as the primordial one, and the remaining ones are partial duplications including triploid-type. The diploid-type repeats are recognized in many 8-tms, 10-tms and 12-tms TM protein sequences, suggesting the diploid-type duplication was a principle mechanism in the evolutionary development of these types of TM proteins. The "positive-inside" rule is satisfied in whole sequences of both 10-tms and 8-tms TM proteins and in both halves of 10-tms proteins while not necessarily in the second half of 8-tms proteins, providing fit examples of "internal divergent topology evolution" likely occurred after a diploid-type internal duplication event. From analyzing the partial duplication patterns, several evolutionary pathways were recognized for 6-tms TM proteins, i.e. from primordial 2-tms, 3-tms and 4-tms TM proteins to extant 6-tms proteins. Similarly, the duplication pattern analysis revealed plausible evolution scenarios that 7-tms TM proteins have arisen from 3-tms, 4-tms and 5-tms TM protein precursors via partial internal gene duplications.     

A Specific interaction between SecA2 and a region of the preprotein adjacent to the signal peptide occurs during transport via the accessory Sec system

The accessory Sec systems of streptococci and staphylococci mediate the transport of a family of large, serine-rich glycoproteins to the bacterial cell surface. These systems are comprised of SecA2, SecY2, and three core accessory Sec proteins (Asp1-3). In Streptococcus gordonii, transport of the serine-rich glycoprotein GspB requires both a unique 90-residue N-terminal signal peptide and an adjacent 24-residue segment (the AST domain). We used in vivo site-specific photo-cross-linking to identify proteins that interact with the AST domain during transport. To facilitate this analysis, the entire accessory Sec system of S. gordonii was expressed in Escherichia coli. The determinants of GspB trafficking to the accessory Sec system in E. coli matched those in S. gordonii, establishing the validity of this approach. When the photo-cross-linker was placed within the AST domain, the preprotein was found to cross-link to SecA2. Importantly, no cross-linking to SecA was detected. Cross-linking of the N-terminal end of the AST domain to SecA2 occurred regardless of whether Asp1-3 were present. However, cross-linking to the C-terminal end was dependent on the Asps. The combined results indicate that full engagement of the AST domain by SecA2 is modulated by one or more of the Asps, and suggest that this process is important for initiating transport.     

Ancestral sequence reconstruction for protein engineers

In addition to its value in the study of molecular evolution, ancestral sequence reconstruction (ASR) has emerged as a useful methodology for engineering proteins with enhanced properties. Proteins generated by ASR often exhibit unique or improved activity, stability, and/or promiscuity, all of which are properties that are valued by protein engineers. Comparison between extant proteins and evolutionary intermediates generated by ASR also allows protein engineers to identify substitutions that have contributed to functional innovation or diversification within protein families. As ASR becomes more widely adopted as a protein engineering approach, it is important to understand the applications, limitations, and recent developments of this technique. This review highlights recent exemplifications of ASR, as well as technical aspects of the reconstruction process that are relevant to protein engineering.     

Use of structural phylogenetic networks for classification of the ferritin-like superfamily

In the postgenomic era, bioinformatic analysis of sequence similarity is an immensely powerful tool to gain insight into evolution and protein function. Over long evolutionary distances, however, sequence-based methods fail as the similarities become too low for phylogenetic analysis. Macromolecular structure generally appears better conserved than sequence, but clear models for how structure evolves over time are lacking. The exponential growth of three-dimensional structural information may allow novel structure-based methods to drastically extend the evolutionary time scales amenable to phylogenetics and functional classification of proteins. To this end, we analyzed 80 structures from the functionally diverse ferritin-like superfamily. Using evolutionary networks, we demonstrate that structural comparisons can delineate and discover groups of proteins beyond the "twilight zone" where sequence similarity does not allow evolutionary analysis, suggesting that considerable and useful evolutionary signal is preserved in three-dimensional structures.     

The Phylogenomic Roots of Modern Biochemistry: Origins of Proteins, Cofactors and Protein Biosynthesis

The complexity of modern biochemistry developed gradually on early Earth as new molecules and structures populated the emerging cellular systems. Here, we generate a historical account of the gradual discovery of primordial proteins, cofactors, and molecular functions using phylogenomic information in the sequence of 420 genomes. We focus on structural and functional annotations of the 54 most ancient protein domains. We show how primordial functions are linked to folded structures and how their interaction with cofactors expanded the functional repertoire. We also reveal protocell membranes played a crucial role in early protein evolution and show translation started with RNA and thioester cofactor-mediated aminoacylation. Our findings allow elaboration of an evolutionary model of early biochemistry that is firmly grounded in phylogenomic information and biochemical, biophysical, and structural knowledge. The model describes how primordial α-helical bundles stabilized membranes, how these were decorated by layered arrangements of β-sheets and α-helices, and how these arrangements became globular. Ancient forms of aminoacyl-tRNA synthetase (aaRS) catalytic domains and ancient non-ribosomal protein synthetase (NRPS) modules gave rise to primordial protein synthesis and the ability to generate a code for specificity in their active sites. These structures diversified producing cofactor-binding molecular switches and barrel structures. Accretion of domains and molecules gave rise to modern aaRSs, NRPS, and ribosomal ensembles, first organized around novel emerging cofactors (tRNA and carrier proteins) and then more complex cofactor structures (rRNA). The model explains how the generation of protein structures acted as scaffold for nucleic acids and resulted in crystallization of modern translation.     

Intra-chain 3D segment swapping spawns the evolution of new multidomain protein architectures

Multidomain proteins form in evolution through the concatenation of domains, but structural domains may comprise multiple segments of the chain. In this work, we demonstrate that new multidomain architectures can evolve by an apparent three-dimensional swap of segments between structurally similar domains within a single-chain monomer. By a comprehensive structural search of the current Protein Data Bank (PDB), we identified 32 well-defined segment-swapped proteins (SSPs) belonging to 18 structural families. Nearly 13% of all multidomain proteins in the PDB may have a segment-swapped evolutionary precursor as estimated by more permissive searching criteria. The formation of SSPs can be explained by two principal evolutionary mechanisms: (i) domain swapping and fusion (DSF) and (ii) circular permutation (CP). By large-scale comparative analyses using structural alignment and hidden Markov model methods, it was found that the majority of SSPs have evolved via the DSF mechanism, and a much smaller fraction, via CP. Functional analyses further revealed that segment swapping, which results in two linkers connecting the domains, may impart directed flexibility to multidomain proteins and contributes to the development of new functions. Thus, inter-domain segment swapping represents a novel general mechanism by which new protein folds and multidomain architectures arise in evolution, and SSPs have structural and functional properties that make them worth defining as a separate group.     

An amino acid sequence which directs membrane insertion causes loss of membrane potential

A 55-amino acid segment, normally present between residues 241 and 295 of the 348-residue gene I protein of the filamentous bacteriophage f1, acts as an internal signal sequence for gene I protein or, when present in fusion proteins, for EcoRI endonuclease or alkaline phosphatase. The resulting proteins are inserted so that they span the membrane with sequences on the amino side of the 55-residue segment in the cytoplasm and those near the carboxy side outside the cytoplasmic membrane. The presence of these proteins in the membrane results in the rapid inhibition of cell growth, probably from a loss of the membrane potential. We describe some of the elements in this 55-residue segment that appear to be crucial for its interaction with the membrane.     

Diversity of human U2AF splicing factors

U2 snRNP auxiliary factor (U2AF) is an essential heterodimeric splicing factor composed of two subunits, U2AF(65) and U2AF(35). During the past few years, a number of proteins related to both U2AF(65) and U2AF(35) have been discovered. Here, we review the conserved structural features that characterize the U2AF protein families and their evolutionary emergence. We perform a comprehensive database search designed to identify U2AF protein isoforms produced by alternative splicing, and we discuss the potential implications of U2AF protein diversity for splicing regulation.     

Artificial proteins from combinatorial approaches

How do we create new artificial proteins? In this review, we present a range of experimental approaches based on combinatorial and directed evolution methods used to explore sequence space and recreate structured or active proteins. These approaches can help to understand constraints of natural evolution and can lead to new useful proteins. Strategies such as binary patterning or modular assembly can efficiently speed structural and functional innovation. Many natural protein architectures are symmetric or repeated and presumably have emerged by coalescence of simpler fragments. This process can be experimentally reproduced; a range of artificial proteins obtained from idealized fragments has recently been described and some of these have already found direct applications.     

Proteins, exons and molecular evolution

The discovery of the eukaryotic gene structure has prompted research into the potential relationship between protein structure and function and the corresponding exon/intron patterns. The exon shuffling hypothesis put forward by Gilbert and Blake suggests the encodement of structural and functional protein elements by exons which can recombine to create novel proteins. This provides an explanation for the relatively rapid evolution of proteins from a few primordial molecules. As the number of gene and protein structures increases, evidence of exon shuffling is becoming more apparent and examples are presented both from modern multi-domain proteins and ancient proteins. Recent work into the chemical properties and catalytic functions of RNA have led to hypotheses based upon the early existence of RNA. These theories suggest that the split gene structure originated in the primordial soup as a result of random RNA synthesis. Stable regions of RNA, or exons, were utilised as primitive enzymes. In response to selective pressures for information storage, the activity was directly transferred from the RNA enzymes or ribozymes, to proteins. These short polypeptides fused together to create larger proteins with a wide range of functions. Recent research into RNA processing and exon size, discussed in this review, provides a clearer insight into the evolutionary development of the gene and protein structure.     

The NMR solution structure of the 30S ribosomal protein S27e encoded in gene RS27_ARCFU of Archaeoglobus fulgidis reveals a novel protein fold

The Archaeoglobus fulgidis gene RS27_ARCFU encodes the 30S ribosomal protein S27e. Here, we present the high-quality NMR solution structure of this archaeal protein, which comprises a C4 zinc finger motif of the CX(2)CX(14-16)CX(2)C class. S27e was selected as a target of the Northeast Structural Genomics Consortium (target ID: GR2), and its three-dimensional structure is the first representative of a family of more than 116 homologous proteins occurring in eukaryotic and archaeal cells. As a salient feature of its molecular architecture, S27e exhibits a beta-sandwich consisting of two three-stranded sheets with topology B(decreasing), A(increasing), F(decreasing), and C(increasing), D(decreasing), E(increasing). Due to the uniqueness of the arrangement of the strands, the resulting fold was found to be novel. Residues that are highly conserved among the S27 proteins allowed identification of a structural motif of putative functional importance; a conserved hydrophobic patch may well play a pivotal role for functioning of S27 proteins, be it in archaeal or eukaryotic cells. The structure of human S27, which possesses a 26-residue amino-terminal extension when compared with the archaeal S27e, was modeled on the basis of two structural templates, S27e for the carboxy-terminal core and the amino-terminal segment of the archaeal ribosomal protein L37Ae for the extension. Remarkably, the electrostatic surface properties of archaeal and human proteins are predicted to be entirely different, pointing at either functional variations among archaeal and eukaryotic S27 proteins, or, assuming that the function remained invariant, to a concerted evolutionary change of the surface potential of proteins interacting with S27.