The archaeo-eukaryotic GINS proteins and the archaeal primase catalytic subunit PriS share a common domain

Primase and GINS are essential factors for chromosomal DNA replication in eukaryotic and archaeal cells. Here we describe a previously undetected relationship between the C-terminal domain of the catalytic subunit (PriS) of archaeal primase and the B-domains of the archaeo-eukaryotic GINS proteins in the form of a conserved structural domain comprising a three-stranded antiparallel β-sheet adjacent to an α-helix and a two-stranded β-sheet or hairpin. The presence of a shared domain in archaeal PriS and GINS proteins, the genes for which are often found adjacent on the chromosome, suggests simple mechanisms for the evolution of these proteins.     

Analysis of 3D structural differences in the IgG-binding domains based on the interresidue average-distance statistics

It is well-known that the IgG-binding domain from staphylococcal protein A folds into a 3α helix bundle structure, while the IgG-binding domain of streptococcal protein G forms an (α + β) structure. Recently, He et al. (Biochemistry 44:14055–14061, 2005) made mutants of these proteins from the wild types of protein A and protein G strains. These mutants are referred to as protein A219 and protein G311, and it was showed that these two mutants have different 3D structures, i.e., the 3α helix bundle structure and the (α + β) structure, respectively, despite the high sequence identity (59%). The purpose of our study was to clarify how such 3D structural differences are coded in the sequences with high homology. To address this problem, we introduce a predicted contact map constructed based on the interresidue average-distance statistics for prediction of folding properties of a protein. We refer to this map as an average distance map (ADM). Furthermore, the statistics of interresidue distances can be converted to an effective interresidue potential. We calculated the contact frequency of each residue of a protein in random conformations with this effective interresidue potential, and then we obtained values similar to ϕ values. We refer to this contact frequency of each residue as a p(μ) value. The comparison of the p(μ) values to the ϕ values for a protein suggests that p(μ) values reveal the information on the folding initiation site. Using these techniques, we try to extract the information on the difference in the 3D structures of protein A219 and protein G311 coded in their amino acid sequences in the present work. The results show that the ADM analyses and the p(μ) value analyses predict the information of folding initiation sites, which can be used to detect the 3D difference in both proteins.     

Physical principles of protein structure and protein folding

Physical principles determining the protein structure and protein folding are reviewed: (i) the molecular theory of protein secondary structure and the method of its prediction based on this theory; (ii) the existence of a limited set of thermodynamically favourable folding patterns of α- and β-regions in a compact globule which does not depend on the details of the amino acid sequence; (iii) the moderns approaches to the prediction of the folding patterns of α- and β-regions in concrete proteins; (iv) experimental approaches to the mechanism of protein folding. The review reflects theoretical and experimental works of the author and his collaborators as well as those of other groups.     

Pathway of Oxidative Folding of a 3-Disulfide α-Lactalbumin May Resemble Either BPTI Model or Hirudin Model

Pathways of oxidative folding of disulfide proteins display a high degree of diversity and vary among two extreme models. The BPTI model is defined by limited species of folding intermediates adopting mainly native disulfide bonds. The hirudin model is characterized by highly heterogeneous folding intermediates containing mostly non-native disulfide bonds. αLA-IIIA is a 3-disulfide variant of α-lactalbumin (αLA) with a 3-D conformation essentially identical to that of intact αLA. αLA-IIIA contains 3 native disulfide bonds of αLA, two of them are located at the calcium binding β-subdomain (Cys⁶¹–Cys⁷⁷ and Cys⁷³–Cys⁹¹) and the third bridge is located within the α-helical domain of the molecule (Cys²⁸–Cys¹¹¹). We investigate here the pathway of oxidative folding of fully reduced αLA-IIIA with and without stabilization of its β-subdomain by calcium binding. In the absence of calcium, the folding pathway of αLA-IIIA was shown to resemble that of hirudin model. Upon stabilization of β-sheet domain by calcium binding, the folding pathway of αLA-IIIA exhibits a striking similarity to that of BPTI model. Three predominant folding intermediates of αLA-IIIA containing exclusively native disulfide bonds were isolated and structurally characterized. Our results further demonstrate that stabilization of subdomains in a protein may dictate its folding pathway and represent a major cause for the existing diversity in the folding pathways of the disulfide-containing proteins.     

Structure-function relationships of heterodimeric amino acid transporters

Heterodimeric amino acid transporters mediate the transfer of amino acids between organs and between different cell types. Members of this particular family of amino acid transporters are constituted by a heavy chain and an associated light chain. The heavy chain is a type II membrane protein with an intracellular amino terminus, a single transmembrane helix, and a large extracellular domain. The light chain, in contrast, is a typical helix-bundle protein with 12 putative transmembrane helices. Two different heavy chains, designated 4F2hc and rbAT, and seven different light chains have been identified to date. Deletion studies indicate that the extracellular domain of the heavy chain has two subdomains. The carboxy-terminal tip of 4F2hc is critical for recognition of certain light chains, whereas the carboxy-terminal tip of rbAT is involved in substrate transport. Sequence alignments suggest that the major part of the extracellular domain forms an α/β domain similar to bacterial α-amylases. A structural model of the rbAT extracellular domain is presented that is in agreement with experimental observations from several mutations and that aligns well with the α-amylase domain.     

A new computational model to study mass inhomogeneity and hydrophobicity inhomogeneity in proteins

We propose a simple yet reliable computational framework that characterizes the differential mass and hydrophobicity distribution within structural classes of proteins. Radial partitioning of protein interior that could successfully distinguish the mass and hydrophobicity distribution patterns in extremophilic proteins from that in their structurally aligned mesophilic counterparts. Distance-dependent mass and hydrophobicity magnitudes could retrieve vital structural insights; needed to probe the hidden connections between packing, folding and stability within different structural classes of proteins, with causality. New computational markers; one, to represent the total mass content; other, related to hydrophobic centrality of proteins, are proposed as well. Results reveal that mass and hydrophobicity packing within extremophilic proteins is indeed more compact than that in their mesophilic counterparts. Analysis of structural constraints within them vindicate it. Total mass (and hydrophobicity) content is found to be maximum in α/β thermophilic proteins and minimum for the all-α mesophilic proteins. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Evidence for myxobacterial origin of eukaryotic defensins

Antimicrobial defensins with the cysteine-stabilized α-helical and β-sheet (CSαβ) motif are a large family of ancient, evolutionarily related innate immunity effectors of multicellular organisms. Although the widespread distribution in plants, fungi, and invertebrates suggests their uniqueness to Eukarya, it is unknown whether these eukaryotic defensins originated before or posterior to the emergence of eukaryotes. In this study, we provide evidence in support of the existence of defensin-like peptides (DLPs) in myxobacteria based on structural bioinformatics analysis, which recognized two bacterial peptides with a conserved cysteine-stabilized α-helical motif, a nested structural unit of the CSαβ motif. Similarity in sequence and structure to fungal DLPs together with restricted distribution to the myxobacteria as well as central role of the myxobacteria in the origin of eukaryotes suggest that the bacterial DLPs represent the ancestor of the eukaryotic defensins and could mediate immune defense of early eukaryotes after gene transfer to the proto-eukaryotic genome. Our work thus offers a basis for further investigation of prokaryotic origin of eukaryotic immune effector molecules. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Structures closed into cycles in globular proteins

Different types of structures closed into cycles are widespread at all the levels of structural organization of proteins. β-Hairpins, triple-stranded β-sheets, and βαβ-units represent simple structural motifs closed into cycles by systems of hydrogen bonds. Secondary closing of these simple motifs into larger cycles by means of different superhelices, split β-hairpins, or SS-bridges results in formation of complex structural motifs such as abcd-units, φ-motifs, five- and seven-segment α/β-motifs, etc. At the level of tertiary structure many proteins and domains fold into structures closed into cylinders. Apparently, closing the motifs and domains into cycles and cylinders results in formation of more cooperative and stable structures as compared with open ones, and this may be the reason for high frequencies of occurrence of the motifs in proteins.     

The NMR structure of the domain II of a chloroplastic NifU-like protein OsNifU1A

NifU-like proteins are a highly conserved protein that serves as the scaffold for assembly of Fe-S clusters. Chloroplastic NifU-like proteins have tandem NifU like domains, named domain I and domain II. Although the amino acid sequences of these domains are very similar to each other, the predicted functional region for the Fe-S cluster assembly, the CXXC motif, exists only in domain I. The structure of the domain II of chloroplastic NifU-like protein OsNifU1A has an α-β sandwich structure containing two α helices located on one side of the β-sheet. The electrostatic surface potential of OsNifU1A domain II is predominantly positively charged. Chloroplastic NifU-like proteins are targeted to ferredoxin for transferring the Fe-S cluster. The ferredoxin presents an overall negatively charged surface, which may evoke an electrostatic association with OsNifU1A domain II.     

Is protein folding rate dependent on number of folding stages? Modeling of protein folding with ferredoxin-like fold

Statistical analysis of protein folding rates has been done for 84 proteins with available experimental data. A surprising result is that the proteins with multi-state kinetics from the size range of 50–100 amino acid residues (a.a.) fold as fast as proteins with two-state kinetics from the same size range. At the same time, the proteins with two-state kinetics from the size range 101–151 a.a. fold faster than those from the size range 50–100 a.a. Moreover, it turns out unexpectedly that usually in the group of structural homologs from the size range 50–100 a.a., proteins with multi-state kinetics fold faster than those with two-state kinetics. The protein folding for six proteins with a ferredoxin-like fold and with a similar size has been modeled using Monte Carlo simulations and dynamic programming. Good correlation between experimental folding rates, some structural parameters, and the number of Monte Carlo steps has been obtained. It is shown that a protein with multi-state kinetics actually folds three times faster than its structural homologs.     

Domain Evolution in the α-Amylase Family

The available amino acid sequences of the α-amylase family (glycosyl hydrolase family 13) were searched to identify their domain B, a distinct domain that protrudes from the regular catalytic (β/α)₈-barrel between the strand β3 and the helix α3. The isolated domain B sequences were inspected visually and also analyzed by Hydrophobic Cluster Analysis (HCA) to find common features. Sequence analyses and inspection of the few available three-dimensional structures suggest that the secondary structure of domain B varies with the enzyme specificity. Domain B in these different forms, however, may still have evolved from a common ancestor. The largest number of different specificities was found in the group with structural similarity to domain B from Bacillus cereus oligo-1,6-glucosidase that contains an α-helix succeeded by a three-stranded antiparallel β-sheet. These enzymes are α-glucosidase, cyclomaltodextrinase, dextran glucosidase, trehalose-6-phosphate hydrolase, neopullulanase, and a few α-amylases. Domain B of this type was observed also in some mammalian proteins involved in the transport of amino acids. These proteins show remarkable similarity with (β/α)₈-barrel elements throughout the entire sequence of enzymes from the oligo-1,6-glucosidase group. The transport proteins, in turn, resemble the animal 4F2 heavy-chain cell surface antigens, for which the sequences either lack domain B or contain only parts thereof. The similarities are compiled to indicate a possible route of domain evolution in the α-amylase family. Received: 4 December 1996 / Accepted: 13 March 1997     

Domain III function of Mu transposase analysed by directed placement of subunits within the transpososome

Assembly of the functional tetrameric form of Mu transposase (MuA protein) at the two att ends of Mu depends on interaction of MuA with multiple att and enhancer sites on supercoiled DNA, and is stimulated by MuB protein. The N-terminal domain I of MuA harbours distinct regions for interaction with the att ends and enhancer; the C-terminal domain III contains separate regions essential for tetramer assembly and interaction with MuB protein (IIIα and IIIβ, respectively). Although the central domain II (the ‘DDE’ domain) of MuA harbours the known catalytic DDE residues, a 26 amino acid peptide within IIIα also has a non-specific DNA binding and nuclease activity which has been implicated in catalysis. One model proposes that active sites for Mu transposition are assembled by sharing structural/catalytic residues between domains II and III present on separate MuA monomers within the MuA tetramer. We have used substrates with altered att sites and mixtures of MuA proteins with either wild-type or altered att DNA binding specificities, to create tetrameric arrangements wherein specific MuA subunits are nonfunctional in II, IIIα or IIIβ domains. From the ability of these oriented tetramers to carry out DNA cleavage and strand transfer we conclude that domain IIIα or IIIβ function is not unique to a specific subunit within the tetramer, indicative of a structural rather than a catalytic function for domain III in Mu transposition.     

Solution NMR structures reveal a distinct architecture and provide first structures for protein domain family PF04536

The protein family (Pfam) PF04536 is a broadly conserved domain family of unknown function (DUF477), with more than 1,350 members in prokaryotic and eukaryotic proteins. High-quality NMR structures of the N-terminal domain comprising residues 41–180 of the 684-residue protein CG2496 from Corynebacterium glutamicum and the N-terminal domain comprising residues 35–182 of the 435-residue protein PG0361 from Porphyromonas gingivalis both exhibit an α/β fold comprised of a four-stranded β-sheet, three α-helices packed against one side of the sheet, and a fourth α-helix attached to the other side. In spite of low sequence similarity (18%) assessed by structure-based sequence alignment, the two structures are globally quite similar. However, moderate structural differences are observed for the relative orientation of two of the four helices. Comparison with known protein structures reveals that the α/β architecture of CG2496(41–180) and PG0361(35–182) has previously not been characterized. Moreover, calculation of surface charge potential and identification of surface clefts indicate that the two domains very likely have different functions.     

Non-erythroid beta spectrin interacting proteins and their effects on spectrin tetramerization

With yeast two-hybrid methods, we used a C-terminal fragment (residues 1697–2145) of non-erythroid beta spectrin (βII-C), including the region involved in the association with alpha spectrin to form tetramers, as the bait to screen a human brain cDNA library to identify proteins interacting with βII-C. We applied stringent selection steps to eliminate false positives and identified 17 proteins that interacted with βII-C (IP_βII-C s). The proteins include a fragment (residues 38–284) of “THAP domain containing, apoptosis associated protein 3, isoform CRA g”, “glioma tumor suppressor candidate region gene 2” (residues 1-478), a fragment (residues 74–442) of septin 8 isoform c, a fragment (residues 704–953) of “coatomer protein complex, subunit beta 1, a fragment (residues 146–614) of zinc-finger protein 251, and a fragment (residues 284–435) of syntaxin binding protein 1. We used yeast three-hybrid system to determine the effects of these βII-C interacting proteins as well as of 7 proteins previously identified to interact with the tetramerization region of non-erythroid alpha spectrin (IP_αII-N s) [1] on spectrin tetramer formation. The results showed that 3 IP_βII-C s were able to bind βII-C even in the presence of αII-N, and 4 IP_αII-N s were able to bind αII-N in the presence of βII-C. We also found that the syntaxin binding protein 1 fragment abolished αII-N and βII-C interaction, suggesting that this protein may inhibit or regulate non-erythroid spectrin tetramer formation.     

Acyl-CoA Dehydrogenases: Dynamic History of Protein Family Evolution

The acyl-CoA dehydrogenases (ACADs) are enzymes that catalyze the α,β-dehydrogenation of acyl-CoA esters in fatty acid and amino acid catabolism. Eleven ACADs are now recognized in the sequenced human genome, and several homologs have been reported from bacteria, fungi, plants, and nematodes. We performed a systematic comparative genomic study, integrating homology searches with methods of phylogenetic reconstruction, to investigate the evolutionary history of this family. Sequence analyses indicate origin of the family in the common ancestor of Archaea, Bacteria, and Eukaryota, illustrating its essential role in the metabolism of early life. At least three ACADs were already present at that time: ancestral glutaryl-CoA dehydrogenase (GCD), isovaleryl-CoA dehydrogenase (IVD), and ACAD10/11. Two gene duplications were unique to the eukaryotic domain: one resulted in the VLCAD and ACAD9 paralogs and another in the ACAD10 and ACAD11 paralogs. The overall patchy distribution of specific ACADs across the tree of life is the result of dynamic evolution that includes numerous rounds of gene duplication and secondary losses, interdomain lateral gene transfer events, alteration of cellular localization, and evolution of novel proteins by domain acquisition. Our finding that eukaryotic ACAD species are more closely related to bacterial ACADs is consistent with endosymbiotic origin of ACADs in eukaryotes and further supported by the localization of all nine previously studied ACADs in mitochondria. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Close Evolutionary Relatedness of α-Amylases from Archaea and Plants

The amino acid sequences of 22 α-amylases from family 13 of glycosyl hydrolases were analyzed with the aim of revealing the evolutionary relationships between the archaeal α-amylases and their eubacterial and eukaryotic counterparts. Two evolutionary distance trees were constructed: (i) the first one based on the alignment of extracted best-conserved sequence regions (58 residues) comprising β2, β3, β4, β5, β7, and β8 strand segments of the catalytic (α/β)₈-barrel and a short conserved stretch in domain B protruding out of the barrel in the β3 →α3 loop, and (ii) the second one based on the alignment of the substantial continuous part of the (α/β)₈-barrel involving the entire domain B (consensus length: 386 residues). With regard to archaeal α-amylases, both trees compared brought, in fact, the same results; i.e., all family 13 α-amylases from domain Archaea were clustered with barley pI isozymes, which represent all plant α-amylases. The enzymes from Bacillus licheniformis and Escherichia coli, representing liquefying and cytoplasmic α-amylases, respectively, seem to be the further closest relatives to archaeal α-amylases. This evolutionary relatedness clearly reflects the discussed similarities in the amino acid sequences of these α-amylases, especially in the best-conserved sequence regions. Since the results for α-amylases belonging to all three domains (Eucarya, Eubacteria, Archaea) offered by both evolutionary trees are very similar, it is proposed that the investigated conserved sequence regions may indeed constitute the ``sequence fingerprints' of a given α-amylase. Received: 3 June 1998 / Accepted: 20 August 1998     

Refolding of the hyperthermophilic protein Ssh10b involves a kinetic dimeric intermediate

The α/β-mixed dimeric protein Ssh10b from the hyperthermophile Sulfolobus shibatae is a member of the Sac10b family that is thought to be involved in chromosomal organization or DNA repair/recombination. The equilibrium unfolding/refolding of Ssh10b induced by denaturants and heat was fully reversible, suggesting that Ssh10b could serve as a good model for folding/unfolding studies of protein dimers. Here, we investigate the folding/unfolding kinetics of Ssh10b in detail by stopped-flow circular dichroism (SF-CD) and using GdnHCl as denaturant. In unfolding reactions, the native Ssh10b turned rapidly into fully unfolded monomers within the stopped-flow dead time with no detectable kinetic intermediate, agreeing well with the results of equilibrium unfolding experiments. In refolding reactions, two unfolded monomers associate in the burst phase to form a dimeric intermediate that undergoes a further, slower, first-order folding process to form the native dimer. Our results demonstrate that the dimerization is essential for maintaining the native tertiary interactions of the protein Ssh10b. In addition, folding mechanisms of Ssh10b and several other α/β-mixed or pure β-sheet proteins are compared.     

Bioinformatic analysis of the TonB protein family

TonB is a protein prevalent in a large number of Gram-negative bacteria that is believed to be responsible for the energy transduction component in the import of ferric iron complexes and vitamin B₁₂ across the outer membrane. We have analyzed all the TonB proteins that are currently contained in the Entrez database and have identified nine different clusters based on its conserved 90-residue C-terminal domain amino acid sequence. The vast majority of the proteins contained a single predicted cytoplasmic transmembrane domain; however, nine of the TonB proteins encompass a ∼290 amino acid N-terminal extension homologous to the MecR1 protein, which is composed of three additional predicted transmembrane helices. The periplasmic linker region, which is located between the N-terminal domain and the C-terminal domain, is extremely variable both in length (22–283 amino acids) and in proline content, indicating that a Pro-rich domain is not a required feature for all TonB proteins. The secondary structure of the C-terminal domain is found to be well preserved across all families, with the most variable region being between the second α-helix and the third β-strand of the antiparallel β-sheet. The fourth β-strand found in the solution structure of the Escherichia coli TonB C-terminal domain is not a well conserved feature in TonB proteins in most of the clusters. Interestingly, several of the TonB proteins contained two C-terminal domains in series. This analysis provides a framework for future structure-function studies of TonB, and it draws attention to the unusual features of several TonB proteins. Byron C. H. Chu and R. Sean Peacock contributed equally to this work.     

Topology, stability, sequence, and length: defining the determinants of two-state protein folding kinetics

The fastest simple, single domain proteins fold a million times more rapidly than the slowest. Ultimately this broad kinetic spectrum is determined by the amino acid sequences that define these proteins, suggesting that the mechanisms that underlie folding may be almost as complex as the sequences that encode them. Here, however, we summarize recent experimental results which suggest that (1) despite a vast diversity of structures and functions, there are fundamental similarities in the folding mechanisms of single domain proteins and (2) rather than being highly sensitive to the finest details of sequence, their folding kinetics are determined primarily by the large-scale, redundant features of sequence that determine a protein's gross structural properties. That folding kinetics can be predicted using simple, empirical, structure-based rules suggests that the fundamental physics underlying folding may be quite straightforward and that a general and quantitative theory of protein folding rates and mechanisms (as opposed to unfolding rates and thus protein stability) may be near on the horizon.