Evolutionary aspects of protein structure and folding

Four basic stages of evolution of protein structure are described based on recent work of the authors targeted specifically on reconstruction of the earliest events in the protein evolution. According to this reconstruction, the initial stage of short peptides of, probably, only few amino-acid residues had been followed by formation of closed loops of the size 25-30 residues, which corresponds to the polymer-statistically optimal ring closure size for mixed polypeptide chains. The next stage involved fusion of the respective small linear genes and formation of protein structures consisting of several closed loops of the nearly standard size, up to 4-6 loops (100-200 amino acid residues) in a typical protein fold. The last, modern stage began with combinatorial fusion of the presumably circular 300-600 bp DNA units and, accordingly, formation of multidomain proteins.     

Distinct Stages of Protein Evolution as Suggested by Protein Sequence Analysis

Evolution of proteins encoded in nucleotide sequences began with the advent of the triplet code. The chronological order of the appearance of amino acids on the evolution scene and the steps in the evolution of the triplet code have been recently reconstructed (Trifonov, 2000b) on the basis of 40 different ranking criteria and hypotheses. According to the consensus chronology, the pair of complementary GGC and GCC codons for the amino acids alanine and glycine appeared first. Other codons appeared as complementary pairs as well, which divided their respective amino acids into two alphabets, encoded by triplets with either central purines or central pyrimidines: G, D, S, E, N, R, K, Q, C, H, Y, and W (Glycine alphabet G) and A, V, P, S, L, T, I, F, and M (Alanine alphabet A). It is speculated that the earliest polypeptide chains were very short, presumably of uniform length, belonging to two alphabet types encoded in the two complementary strands of the earliest mRNA duplexes. After the fusion of the minigenes, a mosaic of the alphabets would form. Traces of the predicted mosaic structure have been, indeed, detected in the protein sequences of complete prokaryotic genomes in the form of weak oscillations with the period 12 residues in the form of alteration of two types of 6 residue long units. The next stage of protein evolution corresponded to the closure of the chains in the loops of the size 25–30 residues (Berezovsky et al., 2000). Autocorrelation analysis of proteins of 23 complete archaebacterial and eubacterial genomes revealed that the preferred distances between valine, alanine, glycine, leucine, and isoleucine along the sequences are in the same range of 25–30 residues, indicating that the loops are primarily closed by hydrophobic interactions between the ends of the loops. The loop closure stage is followed by the formation of typical folds of 100–200 amino acids, via end-to-end fusion of the genes encoding the loop-size chains. This size was apparently dictated by the optimal ring closure for DNA. In both cases the closure into the ring (loop) rendered evolutionarily advantageous stability to the respective structures. Further gene fusions lead to the formation of modern multidomain proteins. Recombinational gene splicing is likely to have appeared after the DNA circularization stage. Received: 21 December 2000 / Accepted: 28 February 2001     

Flowering buds of globular proteins: transpiring simplicity of protein organization

Structural and functional complexity of proteins is dramatically reduced to a simple linear picture when the laws of polymer physics are considered. A basic unit of the protein structure is a nearly standard closed loop of 25-35 amino acid residues, and every globular protein is built of consecutively connected closed loops. The physical necessity of the closed loops had been apparently imposed on the early stages of protein evolution. Indeed, the most frequent prototype sequence motifs in prokaryotic proteins have the same sequence size, and their high match representatives are found as closed loops in crystallized proteins. Thus, the linear organization of the closed loop elements is a quintessence of protein evolution, structure and folding.     

Evolutionary conservation of protein regions in the protonmotive cytochromeb and their possible roles in redox catalysis

Summary The amino acid sequences of the protonmotive cytochromeb from seven representative and phylogenetically diverse species have been compared to identify protein regions or segments that are conserved during evolution. The sequences analyzed included both prokaryotic and eukaryotic examples as well as mitochondrial cytochromeb and chloroplastb ₆ proteins. The principal conclusion from these analyses is that there are five protein regions-each comprising about 20 amino acid residues—that are consistently conserved during evolution. These domains are evident despite the low density of invariant residues. The two most highly conserved regions, spanning approximately consensus residues 130–150 and 270–290, are located in extramembrane loops and are hypothesized to constitute part of the Q_o reaction center. The intramembrane, hydrophobic protein regions containing the heme-ligating histidines are also conserved during evolution. It was found, however, that the conservation of the protein segments extramembrane to the histidine residues ligating the low potential b566 heme group showed a higher degree of sequence conservation. The location of these conserved regions suggests that these extramembrane segments are also involved in forming the Q_o reaction center. A protein segment putatively constituting a portion of the Q_i reaction center, located approximately in the region spanned by consensus residues 20–40, is conserved in species as divergent as mouse andRhodobacter. This region of the protein shows substantially less sequence conservation in the chloroplast cytochromeb ₆. The catalytic role of these conserved regions is strongly supported by locations of residues that are altered in mutants resistant to inhibitors of cytochromeb electron transport.     

Optimized procedures for producing biologically active chemokines

We describe here two strategies to produce biologically active chemokines with authentic N-terminal amino acid residues. The first involves producing the target chemokine with an N-terminal 6×His-SUMO tag in Escherichia coli as inclusion bodies. The fusion protein is solubilized and purified with Ni–NTA–agarose in denaturing reagents. This is further followed by tag removal and refolding in a redox refolding buffer. The second approach involves expressing the target chemokine with an N-terminal 6×His-Trx-SUMO tag in an engineered E. coli strain that facilitates formation of disulfide bonds in the cytoplasm. Following purification of the fusion protein via Ni–NTA and tag removal, the target chemokine is refolded without redox buffer and purified by reverse phase chromatography. Using the procedures, we have produced more than 15 biologically active chemokines, with a yield of up to 15 mg/L.     

Molecular characterization of Alg8, a putative glycosyltransferase, involved in alginate polymerisation

The topology of Alg8, the proposed catalytic subunit of the alginate polymerase, was assessed using PhoA and LacZ fusion protein analysis. This analysis suggested that the periplasmic loop comprises only three amino acid residues with the adjacent transmembrane helices at positions 361–387 and 393–416. Accordingly, the extended cytosolic loop could be located at positions 71–361 and was proposed to contain important catalytic residues. Further experimental evidence for this cytosolic domain was obtained by independently demonstrating this protein region as purified soluble protein domain. The soluble protein domain was identified by MALDI-TOF/MS and presumably represents the cytosolic catalytic domain of Alg8. Site-directed mutagenesis of 11 conserved residues in the cytosolic loop showed that D-188/D-190 (DXD motif), D-295/D-296 (acid–base catalysts) and K-297 were each essential for in vivo polymerase activity, whereas D-179/D-181 (DXD motif), C-244, R-263, D-279, and E-282 were not directly involved in the polymerisation reaction. The role of these amino acid residues with respect to the catalysed alginate polymerisation reaction was discussed with the aid of the recently developed structural model of Alg8.     

Loop fold nature of globular proteins

Protein chains make numerous returns in globules, thus forming loops, closed by tight residue-to-residue contacts-closed loops. Previous statistical analysis of the sizes and locations of the closed loops in all major protein folds revealed that the loops have an almost standard contour length of 25-30 amino acid residues and follow one after another along the chain. In this work the closed loops of the major folds are presented in three dimensions. A special image filtering procedure is introduced that allows one to visualize the standard size closed loops for the first time. The loop positions along the sequences are verified by detection of loop-end clusters.     

Comparison of the computed three-dimensional structures of oncogenic forms (bound to GDP) of theras-Gene-Encoded p21 protein with the structure of the normal (non-transforming) wild-type protein

Theras-oncogene-encoded p21 protein becomes oncogenic if amino acid substitutions occur at critical positions in the polypeptide chain. The most commonly found oncogenic forms contain Val in place of Gly 12 or Leu in place of Gln 61. To determine the effects of these substitutions on the three-dimensional structure of the whole p21 protein, we have performed molecular dynamics calculations on each of these three proteins bound to GDP and magnesium ion to compute the average structures of each of the three forms. Comparisons of the computed average structures shows that both oncogenic forms with Val 12 and Leu 61 differ substantially in structure from that of the wild type (containing Gly 12 and Gln 61) in discrete regions: residues 10–16, 32–47, 55–74, 85–89, 100–110, and 119–134. All of these regions occur in exposed loops, and several of them have already been found to be involved in the cellular functioning of the p21 protein. These regions have also previously been identified as the most flexible domains of the wild-type protein and have been bound to be the same ones that differ in conformation between transforming and nontransforming p21 mutant proteins neither of which binds nucleotide. The two oncogenic forms have similar conformations in their carboxyl-terminal domains, but differ in conformation at residues 32–47 and 55–74. The former region is known to be involved in the interaction with at least three downstream effector target proteins. Thus, differences in structure between the two oncogenic proteins may reflect different relative affinities of each oncogenic protein for each of these effector targets. The latter region, 55–74, is known to be a highly mobile segment of the protein. The results strongly suggest that critical oncogenic amino acid substitutions in the p21 protein cause changes in the structures of vital domains of this protein.     

Thermostable malate synthase of<Emphasis Type="Italic"> Streptomyces thermovulgaris</Emphasis>

The gene, encoding malate synthase (MS), aceB, was cloned from the thermophilic bacterium Streptomyces thermovulgaris by homology-based PCR. The 1,626-bp cloned fragment encodes a protein consisting of 541 amino acids. S. thermovulgaris malate synthase (stMS) gene was over-expressed in Escherichia coli using a glutathione-S transferase (GST) fusion vector (pGEX-6P-1), purified by affinity chromatography, and subsequently cleaved from its GST fusion partner. The purified stMS was characterized and compared to a mesophilic malate synthase (scMS) from Streptomyces coelicolor. stMS exhibited higher temperature optima (40–60 °C) than those of scMS (28–37 °C). It was more thermostable and very resistant to the chemical denaturant urea. Amino acid sequence comparison of stMS with four mesophilic streptomycete MSs indicated that they share 70.9–91.4% amino acid identities, with stMS possessing slightly more charged residues (~31%) than its mesophilic counterparts (~28–29%). Seven charged residues (E85, R187, R209, H239, H364, R382 and K520) that were unique to stMS may be selectively and strategically placed to support its peculiar characteristics.     

Gene assembly consisting of small units with similar amino acid composition in the <Emphasis Type="Italic">Saccharomyces cerevisiae</Emphasis> genome

Amino acid compositions of all genes in Saccharomyces cerevisiae were determined using a computer analysis of the complete genome. The amino acid composition of an assembly of several genes or a single gene consisted of 3000–7000 amino acid residues forming a certain pattern of amino acid composition. This rule was independent not only of the gene size, but also of the gene position. Thus, the small units, consisting of 3000–7000 amino acid residues, showed a similar amino acid composition, and they formed all the genes in the complete genome.     

Protein Modules Conserved Since LUCA

Universal scale of the sequence conservation has been recently introduced based on omnipresence of the protein sequence motifs across species. A large spectrum of short sequences, up to eight residues has been found to reside in all or almost all prokaryotic organisms. By this discovery a principally novel quantitative approach is introduced to the problem of reconstruction of the last universal common ancestor (LUCA). The most conserved elements (protein modules) with defined structures and sequences harboring the omnipresent motifs are outlined in this work, by combining the sequence and protein crystal structure data. The structurally conserved modules involve 25–30 amino acid residues and have appearance of closed loops, loop-n-lock structures. This confirms earlier conclusions on the loop-fold structure of globular proteins. Many of the topmost conserved modules represent the primary closed loop prototypes, that have been derived by whole genome sequence searches. The data presented, thus, make a basis for further developments toward the earliest stages of protein evolution. [Reviewing Editor: Dr. Martin Kreitman]     

Protein folding: looping from hydrophobic nuclei.

Protein structure can be viewed as a compact linear array of nearly standard size closed loops of 25-30 amino acid residues (Berezovsky et al., FEBS Letters 2000; 466: 283-286) irrespective of details of secondary structure. The end-to-end contacts in the loops are likely to be hydrophobic, which is a testable hypothesis. This notion could be verified by direct comparison of the loop maps with Kyte and Doolittle hydropathicity plots. This analysis reveals that most of the ends of the loops are hydrophobic, indeed. The same conclusion is reached on the basis of positional autocorrelation analysis of protein sequences of 23 fully sequenced bacterial genomes. Hydrophobic residues valine, alanine, glycine, leucine, and isoleucine appear preferentially at the 25-30 residues distance one from another. These observations open a new perspective in the understanding of protein structure and folding: a consecutive looping of the polypeptide chain with the loops ending primarily at hydrophobic nuclei.     

Collapsin response mediator protein-2 regulates neurite formation by modulating tubulin GTPase activity

Collapsin response mediator protein-2 (CRMP-2) plays a key role in axonal development by regulating microtubule dynamics. However, the molecular mechanisms underlying this function have not been clearly elucidated. In this study, we demonstrated that hCRMP-2, specifically amino acid residues 480–509, is essential for stimulating tubulin GTPase activity. We also found that the GTPase-activating protein (GAP) activity of hCRMP-2 was important for microtubule assembly and neurite formation in differentiated PC12 pheochromocytoma cell lines. Mutant hCRMP-2, lacking arginine residues responsible for GAP activity, inhibited microtubule assembly and neurite formation. Interestingly, we found that the N-terminal region (amino acids150–299) of hCRMP-2 had an inhibitory role on GAP activity via a direct interaction with the C-terminal region (amino acids 480–509). Our results suggest that CRMP-2 as a tubulin direct binder may be a GAP of tubulin in neurite formation and that its GAP activity may be regulated by an intramolecular interaction with an N-terminal inhibitory region.     

DNA binding domains in Tn3 transposase

Summary Various segments of Tn3 transposase were fused individually to -galactosidase, and the resulting fusion proteins were examined for their DNA binding ability by a nitrocellulose filter binding assay. Analyses of a series of the fusion proteins revealed that the N-terminal segment of the transposase (amino acid positions 1–242; the transposase gene encodes 1004 residues in all) had specific DNA binding ability for the 38 bp terminal inverted repeat (IR) sequence, and the central segment (amino acid positions 243–632) had non-specific DNA binding ability. Further analyses of each of the two regions revealed that the N-terminal segment could be divided into at least two subsegments (amino acid positions 1–86 and 87–242), neither of which had specific DNA binding ability, but which both possessed nonspecific DNA binding ability. The central segment included two subsegments (amino acid positions 243–289 and 439–505) with non-specific DNA binding ability. These results and other observations suggest that Tn3 transposase has several domains including those responsible for non-specific DNA binding, and a combination of two or more domains gives rise to specific DNA binding activity. The C-terminal segment of the transposase (amino acid positions 633-1004), which is very well conserved among transposases encoded by Tn3 family transposons, had no DNA binding ability. This segment may represent the main part of the catalytic domain responsible for the initiation step of transposition.     

Sequence Structure of van der Waals Locks in Proteins

Abstract

Recent works has suggested that proteins in early evolution have gone through a stage of closed loop elements with a typical contour size of 25–35 residues. These closed loops are still the elementary protein units to these days, and can be used to spell out protein sequence/structure relationship through a relatively small number of protein prototypes. In this study we aimed to identify the sequences that are used to lock the loop ends to one another, and to show how an extensive dictionary of such locking pairs can be created using positional correlation data from a large proteome database, and structural data from PDB databases. Such a dictionary can be used in reconstructing the evolutionary pathway the modern proteins have gone through, and in identifying closed loop elements in modern proteins with yet unknown 3D structure.     

Features of a Spatially Constrained Cystine Loop in the p10 FAST Protein Ectodomain Define a New Class of Viral Fusion Peptides

The reovirus fusion-associated small transmembrane (FAST) proteins are the smallest known viral membrane fusion proteins. With ectodomains of only ∼20–40 residues, it is unclear how such diminutive fusion proteins can mediate cell-cell fusion and syncytium formation. Contained within the 40-residue ectodomain of the p10 FAST protein resides an 11-residue sequence of moderately apolar residues, termed the hydrophobic patch (HP). Previous studies indicate the p10 HP shares operational features with the fusion peptide motifs found within the enveloped virus membrane fusion proteins. Using biotinylation assays, we now report that two highly conserved cysteine residues flanking the p10 HP form an essential intramolecular disulfide bond to create a cystine loop. Mutagenic analyses revealed that both formation of the cystine loop and p10 membrane fusion activity are highly sensitive to changes in the size and spatial arrangement of amino acids within the loop. The p10 cystine loop may therefore function as a cystine noose, where fusion peptide activity is dependent on structural constraints within the noose that force solvent exposure of key hydrophobic residues. Moreover, inhibitors of cell surface thioreductase activity indicate that disruption of the disulfide bridge is important for p10-mediated membrane fusion. This is the first example of a viral fusion peptide composed of a small, spatially constrained cystine loop whose function is dependent on altered loop formation, and it suggests the p10 cystine loop represents a new class of viral fusion peptides.     

Crystal Structure of OXA-58 with the Substrate-Binding Cleft in a Closed State: Insights into the Mobility and Stability of the OXA-58 Structure

OXA-58 is a class D β-lactamase from the multi-drug resistant Acinetobacter baumannii. We determined the crystal structure of OXA-58 in a novel crystal, and revealed the structure of the substrate-binding cleft in a closed state, distinct from a previously reported OXA-58 crystal structure with the binding cleft in an open state. In the closed state, the movement of three loops (α3–α4, β6–β7, and β8–α10) forms an arch-like architecture over the binding cleft through interaction between the Phe113 residues of α3–α4 and Met225 of β6–β7. This structure suggests the involvement of these flexible loops in OXA-58 substrate binding. In contrast to the mobile loops, the Ω-loop appeared static, including the conserved loop residues and their hydrogen bonds; the pivotal residue Trp169 within the Ω-loop, ζ-carbamic acid of the modified base catalyst residue Lys86, and nucleophilic residue Ser83. The stability of OXA-58 was enhanced concomitant with an increase in the hydrolytic activity catalyzed by NaHCO₃-dependent ζ-carbamic acid formation, with an EC₅₀ of 0.34 mM. The W169A mutant enzyme was significantly thermally unstable even in the presence of 100 mM NaHCO₃, whereas the S83A mutant was stabilized with NaHCO₃-dependent activation. The ζ-carbamic acid was shown to increase not only OXA-58 hydrolytic activity but also OXA-58 stability through the formation of a hydrogen bond network connected to the Ω-loop with Ser83 and Trp169. Thus, the static Ω-loop is important for OXA-58 stability, whereas the mobile loops of the substrate-binding cleft form the basis for accommodation of the various substituents of β-lactam backbone.     

A study of extracellular matrix-cell adhesion peptidic epitopes related to human serum amyloid A (SAA)

Serum amyloid A (SAA), an acute-phase protein, exists normally in the serum while complexed with high-density lipoprotein 3 (SAA- HDL3). Its levels increase markedly during inflammatory diseases. The pentapeptide Tyr-Ile-Gly-Ser-Asp (YIGSR-like) and the tripeptide Arg-Gly-Asn (RGD-like), related to the cell adhesion domains of laminin and fibronectin, respectively, exist in SAA within close proximity (YIGSDKYFHARGNY; amino acid residues 29–42). A structure-function study of linear and head- to-tail cyclic peptides, related to the amino acid residues 29–42 and 70–76 (GRGAEDS) of human SAA, was performed in order to evaluate their ability to inhibit adhesion of human T-lymphocytes to surfaces coated with extracellular matrix purified from bovine corneal cells.     

Segmented structure of protein sequences and early evolution of genome by combinatorial fusion of DNA elements

A theory of an early stage of genome evolution by combinatorial fusion of circular DNA units is suggested, based on protein sequence fossil evidence. The evidence includes preference of protein sequence lengths for certain sizes—multiples of 123 as for eukaryotes and multiples of 152 as for prokaryotes. At the DNA level these sizes correspond to 350–450 base pairs—the known optimal range for DNA ring closure. The methionine residues repeatedly appear along the sequences with the same period of about 120 as (in eukaryotes), presumably marking the sites of insertion of the early genes—rings of protein-coding DNA. No torsional constraint in this DNA results in very sharp estimate of the helical periodicity of the early DNA, indistinguishable from the experimental mean value for extant DNA. According to the combinatorial fusion theory, based on the above evidence, in the pregenomic, prerecombinational stage the genes and the noncoding sequences existed in form of autonomously replicating DNA rings of close to standard size, randomly segregating between dividing cells, like modern plasmids do. In the recombinational early genomic stage the rings started to fuse, forming larger DNA molecules consisting of several unit genes connected in various combinations and forming long protein-coding sequences (combinatorial fusion). This process, which involved, perhaps, noncoding sequences as well, eventually resulted in the formation of large genomes. The dispersed circular DNA—or, rather, evolutionarily advanced derivatives thereof—may still exist in the form of various mobile DNA elements.     

Sequence structure of van der Waals locks in proteins

Recent works has suggested that proteins in early evolution have gone through a stage of closed loop elements with a typical contour size of 25-35 residues. These closed loops are still the elementary protein units to these days, and can be used to spell out protein sequence/structure relationship through a relatively small number of protein prototypes. In this study we aimed to identify the sequences that are used to lock the loop ends to one another, and to show how an extensive dictionary of such locking pairs can be created using positional correlation data from a large proteome database, and structural data from PDB databases. Such a dictionary can be used in reconstructing the evolutionary pathway the modern proteins have gone through, and in identifying closed loop elements in modern proteins with yet unknown 3D structure.