Modeling protein network evolution under genome duplication and domain shuffling

Background  

Successive whole genome duplications have recently been firmly established in all major eukaryote kingdoms. Such exponential evolutionary processes must have largely contributed to shape the topology of protein-protein interaction (PPI) networks by outweighing, in particular, all time-linear network growths modeled so far.     

Light signaling mediated by PAS domain-containing proteins in Xanthomonas campestris pv. campestris

Per-ARNT-Sim (PAS) domains are important signalling modules that possibly monitor changes in various stimuli such as light. For the majority of PAS domains that have been identified by sequence similarity, the biological function of the signalling pathways has not yet been experimentally investigated.Thirty-three PAS proteins were discovered in Xanthomonas campestris pv. campestris(Xcc) by genome/proteome analysis. Thirteen PAS proteins were identified as contributing to light signalling and Xcc growth, motility or virulence using molecular genetics and bioinformatics methods. The PAS domains played important roles in light signalling to regulate the growth, motility and virulence of Xcc. They might be regulated by not only light quality (wavelength)but also quantity (intensity) as potential light-signalling components. Evaluating the light wavelength, three light-signalling types of PAS proteins in Xcc were shown to be involved in blue light signalling, tricolour (blue, red and far red)signalling or red/far-red signalling. This showed that Xcc had evolved a complicated light-signalling system to adapt to a complex environment.     

Directed evolution of proteins by exon shuffling

Evolution of eukaryotes is mediated by sexual recombination of parental genomes. Crossovers occur in random, but homologous, positions at a frequency that depends on DNA length. As exons occupy only 1% of the human genome and introns about 24%, by far most of the crossovers occur between exons, rather than inside. The natural process of creating new combinations of exons by intronic recombination is called exon shuffling. Our group is developing in vitro formats for exon shuffling and applying these to the directed evolution of proteins. Based on the splice frame junctions, nine classes of exons and three classes of introns can be distinguished. Splice frame diagrams of natural genes show how the splice frame rules govern exon shuffling. Here, we review various approaches to constructing libraries of exon-shuffled genes. For example, exon shuffling of human pharmaceutical proteins can generate libraries in which all of the sequences are fully human, without the point mutations that raise concerns about immunogenicity.     

Molecular cloning and characterization of hetR genes from filamentous cyanobacteria

HetR, a serine type protease, plays an important role in heterocyst differentiation in filamentous cyanobacteria. We isolated and sequenced the hetR genes from different heterocystous and filamentous nonheterocystous cyanobacteria. The hetR gene in the heterocyst forming Anabaena variabilis ATCC 29413 FD was interrupted by interposon mutagenesis (mutant strain WSIII8). This mutant does not form heterocysts and shows no diazotrophic growth under aerobic conditions. However, under anaerobic N(2)-fixing conditions, the WSIII8 cells are able to grow, and high nitrogenase (Nif2) activity is detectable. Nif2 expression was demonstrated in each vegetative cell of the filament by immunolocalization 4 h after nitrogen step-down.     

The hypoxic testis and post-meiotic expression of PAS domain proteins

Spermatogenesis in the seminiferous tubuli of the testis occurs under a high proliferation rate, suggesting considerable oxygen consumption. Because of the lack of blood vessels, the oxygen partial pressure in the lumen of the tubuli is very low. However, the consequences of these environmental conditions on spermatogenesis are unknown. The PAS domain is found in environmental protein sensors involved in the perception of oxygen partial pressure, light intensity, redox potentials, voltage and certain ligands. We previously identified two PAS proteins highly expressed in the testis: a novel isoform of the hypoxia-inducible factor (HIF)-1alpha and PASKIN, a PAS-Ser/Thr kinase related to bacterial oxygen sensing PAS-domain proteins.     

Conservation of telomere protein complexes: shuffling through evolution

The rapid evolution of telomere proteins has hindered identification of orthologs from diverse species and created the impression that certain groups of eukaryotes have largely non-overlapping sets of telomere proteins. However, the recent identification of additional telomere proteins from various model organisms has dispelled this notion by expanding our understanding of the composition, architecture and range of telomere protein complexes present in individual species. It is now apparent that versions of the budding yeast CST complex and mammalian shelterin are present in multiple phyla. While the precise subunit composition and architecture of these complexes vary between species, the general function is often conserved. Despite the overall conservation of telomere protein complexes, there is still considerable species-specific variation, with some organisms having lost a particular subunit or even an entire complex. In some cases, complex components appear to have migrated between the telomere and the telomerase RNP. Finally, gene duplication has created telomere protein paralogs with novel functions. While one paralog may be part of a conserved telomere protein complex and have the expected function, the other paralog may serve in a completely different aspect of telomere biology.     

Envelope structure of four gliding filamentous cyanobacteria.

The cell walls of four gliding filamentous Oscillatoriaceae species comprising three different genera were studied by freeze substitution, freeze fracturing, and negative staining. In all species, the multilayered gram-negative cell wall is covered with a complex external double layer. The first layer is a tetragonal crystalline S-layer anchored on the outer membrane. The second array is formed by parallel, helically arranged surface fibrils with diameters of 8 to 12 nm. These fibrils have a serrated appearance in cross sections. In all cases, the orientation of the surface fibrils correlates with the sense of revolution of the filaments during gliding, i.e., clockwise in both Phormidium strains and counterclockwise in Oscillatoria princeps and Lyngbya aeruginosa. The lack of longitudinal corrugations or contractions of the surface fibrils and the identical appearances of motile and nonmotile filaments suggest that this structure plays a passive screw thread role in gliding. It is hypothesized that the necessary propulsive force is generated by shear forces between the surface fibrils and the continuing flow of secreted extracellular slime. Furthermore, the so-called junctional pores seem to be the extrusion sites of the slime. In motile cells, these pores exhibit a different staining behavior than that seen in nonmotile ones. In the former, the channels of the pores are filled with electron-dense material, whereas in the latter, the channels appear comparatively empty, highly contrasting the peptidoglycan. Finally, the presence of regular surface structures in other gliding prokaryotes is considered an indication that comparable structures are general features of the cell walls of gliding microbes.     

Use of filamentous cyanobacteria for biodegradation of organic pollutants.

Biodegradation is increasingly being considered as a less expensive alternative to physical and chemical means of decomposing organic pollutants. Pathways of biodegradation have been characterized for a number of heterotrophic microorganisms, mostly soil isolates, some of which have been used for remediation of water. Because cyanobacteria are photoautotrophic and some can fix atmospheric nitrogen, their use for bioremediation of surface waters would circumvent the need to supply biodegradative heterotrophs with organic nutrients. This paper demonstrates that two filamentous cyanobacteria have a natural ability to degrade a highly chlorinated aliphatic pesticide, lindane (gamma-hexachlorocyclohexane); presents quantitative evidence that this ability can be enhanced by genetic engineering; and provides qualitative evidence that those two strains can be genetically engineered to degrade another chlorinated pollutant, 4-chlorobenzoate.     

Occurrence of facultative anoxygenic photosynthesis among filamentous and unicellular cyanobacteria.

Eleven of 21 cyanobacteria strains examined are capable of facultative anoxygenic photosynthesis, as shown by their ability to photoassimilate CO2 in the presence of Na2S, 3-(3,4-dichlorophenyl)-1,1-dimethylurea and 703-nm light. These include different cyanobacterial types (filamentous and unicellular) of different growth histories (aerobic, anaerobic, and marine and freshwater). Oscillatoria limnetica, Aphanothece halophytica (7418), and Lyngbya (7104) have different optimal concentrations of Na2S permitting CO2 photoassimilation, above which the rate decreases: 3.5, 0.7, and 0.1 mM, respectively. In A. halophytica, for each CO2 molecule photoassimilated two sulfide molecules are oxidized to elemental sulfur, which is excreted from the cells.The ecological and evolutionary significance of anoxygenic photosynthesis in the cyanobacteria is discussed.     

Response regulators of bacterial signal transduction systems: Selective domain shuffling during evolution

Response regulators of bacterial sensory transduction systems generally consist of receiver module domains covalently linked to effector domains. The effector domains include DNA binding and/or catalytic units that are regulated by sensor kinase-catalyzed aspartyl phosphorylation within their receiver modules. Most receiver modules are associated with three distinct families of DNA binding domains, but some are associated with other types of DNA binding domains, with methylated chemotaxis protein (MCP) demethylases, or with sensor kinases. A few exist as independent entities which regulate their target systems by noncovalent interactions.In this study the molecular phylogenies of the receiver modules and effector domains of 49 fully sequenced response regulators and their homologues were determined. The three major, evolutionarily distinct, DNA binding domains found in response regulators were evaluated for their phylogenetic relatedness, and the phylogenetic trees obtained for these domains were compared with those for the receiver modules. Members of one family (family 1) of DNA binding domains are linked to large ATPase domains which usually function cooperatively in the activation of E. Coli ⁵⁴-dependent promoters or their equivalents in other bacteria. Members of a second family (family 2) always function in conjunction with the E. Coli ⁷⁰ or its equivalent in other bacteria. A third family of DNA binding domains (family 3) functions by an uncharacterized mechanism involving more than one a factor. These three domain families utilize distinct helix-turn-helix motifs for DNA binding.The phylogenetic tree of the receiver modules revealed three major and several minor clusters of these domains. The three major receiver module clusters (clusters 1, 2, and 3) generally function with the three major families of DNA binding domains (families 1, 2, and 3, respectively) to comprise three classes of response regulators (classes 1, 2, and 3), although several exceptions exist. The minor clusters of receiver modules were usually, but not always, associated with other types of effector domains. Finally, several receiver modules did not fit into a cluster. It was concluded that receiver modules usually diverged from common ancestral protein domains together with the corresponding effector domains, although domain shuffling, due to intragenic splicing and fusion, must have occurred during the evolution of some of these proteins.Multiple sequence alignments of the 49 receiver modules and their various types of effector domains, together with other homologous domains, allowed definition of regions of striking sequence similarity and degrees of conservation of specific residues. Sequence data were correlated with structure/function when such information was available. These studies should provide guides for extrapolation of results obtained with one response regulator to others as well as for the design of future structure/function analyses. Correspondence to: M.H. Saier, Jr.     

Survey of extrachromosomal DNA found in the filamentous cyanobacteria.

Cleared lysates of 13 species of filamentous cyanobacteria were examined for the presence of extrachromosomal DNA by using agarose gel electrophoresis and ethidium bromide staining. Seven of the 13 species contained extrachromosomal covalently closed circular DNA, and all but 1 species contained multiple elements. There was no correlation between the presence of extrachomosomal DNA and either the range of metabolic activities found in the cyanobacteria or the differentiated cell types or structures elaborated by the morphologically complex filamentous cyanobacteria.     

In vitro rapid evolution of fungal immunomodulatory proteins by DNA family shuffling

Fungal immunomodulatory proteins (FIPs) found in a wide variety of mushrooms hold significant therapeutic potential. Despite much research, the structural determinants for their immunomodulatory functions remain unknown. In this study, a DNA shuffling technique was used to create two shuffled FIP protein libraries: an intrageneric group containing products of shuffling between FIP-glu (FIP gene isolated from Ganoderma lucidum) and FIP-gsi (FIP gene isolated from Ganoderma sinense) genes and an intergeneric group containing the products of shuffling between FIP-glu, FIP-fve (FIP gene isolated from Flammulina velutipes), and FIP-vvo (FIP gene isolated from Volvariella volvacea) genes. The gene shuffling generated 426 and 412 recombinant clones, respectively. Using colony blot analysis, we selected clones that expressed relatively high levels of shuffled gene products recognized by specific polyclonal antibodies. We analyzed the DNA sequences of the selected shuffled genes, and testing of their protein products revealed that they maintained functional abilities to agglutinate blood cells and induce cytokine production by splenocytes from Kunming mice in vitro. Meanwhile, the relationships between protein structure and the hemagglutination activity and between the changed nucleotide sites and expression levels were explored by bioinformatic analysis. These combined analyses identified the nucleotide changes involved in regulating the expression levels and hemagglutination activities of the FIPs. Therefore, we were able to generate recombinant FIPs with improved biological activities and expression levels by using DNA shuffling, a powerful tool for the generation of novel therapeutic proteins and for their structural and functional studies.     

The PAS fold. A redefinition of the PAS domain based upon structural prediction.

In the postgenomic era it is essential that protein sequences are annotated correctly in order to help in the assignment of their putative functions. Over 1300 proteins in current protein sequence databases are predicted to contain a PAS domain based upon amino acid sequence alignments. One of the problems with the current annotation of the PAS domain is that this domain exhibits limited similarity at the amino acid sequence level. It is therefore essential, when using proteins with low-sequence similarities, to apply profile hidden Markov model searches for the PAS domain-containing proteins, as for the PFAM database. From recent 3D X-ray and NMR structures, however, PAS domains appear to have a conserved 3D fold as shown here by structural alignment of the six representative 3D-structures from the PDB database. Large-scale modelling of the PAS sequences from the PFAM database against the 3D-structures of these six structural prototypes was performed. All 3D models generated (> 5700) were evaluated using prosaii. We conclude from our large-scale modelling studies that the PAS and PAC motifs (which are separately defined in the PFAM database) are directly linked and that these two motifs form the PAS fold. The existing subdivision in PAS and PAC motifs, as used by the PFAM and SMART databases, appears to be caused by major differences in sequences in the region connecting these two motifs. This region, as has been shown by Gardner and coworkers for human PAS kinase (Amezcua, C.A., Harper, S.M., Rutter, J. & Gardner, K.H. (2002) Structure 10, 1349-1361, [1]), is very flexible and adopts different conformations depending on the bound ligand. Some PAS sequences present in the PFAM database did not produce a good structural model, even after realignment using a structure-based alignment method, suggesting that these representatives are unlikely to have a fold resembling any of the structural prototypes of the PAS domain superfamily.     

Molecular evolution of the thermosensitive PAb1620 epitope of human p53 by DNA shuffling.

Conformational stability of the p53 protein is an absolute necessity for its physiological function as a tumor suppressor. Recent in vitro studies have shown that wild-type p53 is a highly temperature-sensitive protein at the structural and functional levels. Upon heat treatment at 37 degrees C, p53 loses its wild-type (PAb1620(+)) conformation and its ability to bind DNA, but can be stabilized by different classes of ligands. To further investigate the thermal instability of p53, we isolated p53 mutants resistant to heat denaturation. For this purpose, we applied a recently developed random mutagenesis technique called DNA shuffling and screened for p53 variants that could retain reactivity to the native conformation-specific anti-p53 antibody PAb1620 upon thermal treatment. After three rounds of mutagenesis and screening, mutants were isolated with the desired phenotype. The isolated mutants were translated in vitro in either Escherichia coli or rabbit reticulocyte lysate and characterized biochemically. Mutational analysis identified 20 amino acid residues in the core domain of p53 (amino acids 101-120) responsible for the thermostable phenotype. Furthermore, the thermostable mutants could partially protect the PAb1620(+) conformation of tumor-derived p53 mutants from thermal unfolding, providing a novel approach for restoration of wild-type structure and possibly function to a subset of p53 mutants in tumor cells.     

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.     

Evolution of EF-hand calcium-modulated proteins. IV. Exon shuffling did not determine the domain compositions of EF-hand Proteins

Summary In the previous three reports in this series we demonstrated that the EF-hand family of proteins evolved by a complex pattern of gene duplication, transposition, and splicing. The dendrograms based on exon sequences are nearly identical to those based on protein sequences for troponin C, the essential light chain myosin, the regulatory light chain, and calpain. This validates both the computational methods and the dendrograms for these subfamilies. The proposal of congruence for calmodulin, troponin, C, essential light chain, and regulatory light chain was confirmed. There are, however, significant differences in the calmodulin dendrograms computed from DNA and from protein sequences. In this study we find that introns are distributed throughout the EF-hand domain and the interdomain regions. Further, dendrograms based on intron type and distribution bear little resemblance to those based on protein or on DNA sequences. We conclude that introns are inserted, and probably deleted, with relatively high frequency. Further, in the EF-hand family exons do not correspond to structural domains and exon shuffling played little if any role in the evolution of this widely distributed homolog family. Calmodulin has had a turbulent evolution. Its dendrograms based on protein sequence, exon sequence, 3′-tail sequence, intron sequences, and intron positions all show significant differences.     

Characterization and cloning of extrachromosomal DNA from filamentous cyanobacteria

Abstract Physical maps of cryptic plasmids from the filamentous cyanobacteria Anabaena variabilis PCC7118 (pGL1: 3.6 MDa), Nostoc PCC6705 (pGL2: 2.6 MDa) and Plectonema PCC6306 (pGL3: 0.95 MDa) were generated. Selectable markers were introduced onto pGL2 and pGL3 by fusing them to the vector pBR328, using their single restriction sites for Cla I. The recombinant plasmids generated were characterised with respect to the orientation of the insert and the single sites for restriction endonucleases which they possess. The stability of pGL1 and of the two recombinant plasmids in culture was investigated and a method for isolating larger cyanobacterial plasmids (> 20 MDa) was devised.     

RNA association or phosphorylation of the RS domain prevents aggregation of RS domain-containing proteins

Domains rich in alternating arginine and serine residues (RS domains) are found in a large number of eukaryotic proteins involved in several cellular processes. According to the prevailing view RS domains function as protein interaction domains, thereby promoting the assembly of higher-order cellular structures. Furthermore, recent data demonstrated that the RS regions of several SR splicing factors directly contact the pre-mRNA in a nonsequence specific but functionally important fashion. Using a variety of biochemical approaches, we now demonstrate that the RS domains of three proteins, not directly associated with the splicing reaction, such as lamin b receptor, acinus and peroxisome proliferator-activated receptor gamma coactivator-1 alpha, associate mainly with nuclear RNA and that this association is conducive in retaining the proteins in a soluble form. Phosphorylation by SRPK1 prevents RNA association, yet it greatly increases the fraction of the proteins recovered in soluble form, thereby mimicking the RNA effect. Based on these results we propose that the tendency to self-associate and form aggregates is a general property of RS domain-containing proteins and could be attributed to their disordered structure. RNA binding or SRPK1-mediated phosphorylation prevents aggregation and may serve to modulate the RS domain interaction modes.