Evolutionary comparisons of RecA-like proteins across all major kingdoms of living organisms

Protein sequences with similarities to Escherichia coli RecA were compared across the major kingdoms of eubacteria, archaebacteria, and eukaryotes. The archaeal sequences branch monophyletically and are most closely related to the eukaryotic paralogous Rad51 and Dmc1 groups. A multiple alignment of the sequences suggests a modular structure of RecA-like proteins consisting of distinct segments, some of which are conserved only within subgroups of sequences. The eukaryotic and archaeal sequences share an N-terminal domain which may play a role in interactions with other factors and nucleic acids. Several positions in the alignment blocks are highly conserved within the eubacteria as one group and within the eukaryotes and archaebacteria as a second group, but compared between the groups these positions display nonconservative amino acid substitutions. Conservation within the RecA-like core domain identifies possible key residues involved in ATP-induced conformational changes. We propose that RecA-like proteins derive evolutionarily from an assortment of independent domains and that the functional homologs of RecA in noneubacteria comprise an array of RecA-like proteins acting in series or cooperatively.     

The preference for GT-rich DNA by the yeast Rad51 protein defines a set of universal pairing sequences

The Rad51 protein of Saccharomyces cerevisiae is a eukaryotic homolog of the RecA protein, the prototypic DNA strand-exchange protein of Escherichia coli. RAD51 gene function is required for efficient genetic recombination and for DNA double-strand break repair. Recently, we demonstrated that RecA protein has a preferential affinity for GT-rich DNA sequences-several of which exhibit enhanced RecA protein-promoted homologous pairing activity. The fundamental similarity between the RecA and Rad51 proteins suggests that Rad51 might display an analogous bias. Using in vitro selection, here we show that the yeast Rad51 protein shares the same preference for GT-rich sequences as its prokaryotic counterpart. This bias is also manifest as an increased ability of Rad51 protein to promote the invasion of supercoiled DNA by homologous GT-rich single-stranded DNA, an activity not previously described for the eukaryotic pairing protein. We propose that the preferred utilization of GT-rich sequences is a conserved feature among all homologs of RecA protein, and that GT-rich regions are loci for increased genetic exchange in both prokaryotes and eukaryotes.     

Heat shock protein 70 family: multiple sequence comparisons, function, and evolution

The heat shock protein 70 kDa sequences (HSP70) are of great importance as molecular chaperones in protein folding and transport. They are abundant under conditions of cellular stress. They are highly conserved in all domains of life: Archaea, eubacteria, eukaryotes, and organelles (mitochondria, chloroplasts). A multiple alignment of a large collection of these sequences was obtained employing our symmetric-iterative ITERALIGN program (Brocchieri and Karlin 1998). Assessments of conservation are interpreted in evolutionary terms and with respect to functional implications. Many archaeal sequences (methanogens and halophiles) tend to align best with the Gram-positive sequences. These two groups also miss a signature segment [about 25 amino acids (aa) long] present in all other HSP70 species (Gupta and Golding 1993). We observed a second signature sequence of about 4 aa absent from all eukaryotic homologues, significantly aligned in all prokaryotic sequences. Consensus sequences were developed for eight groups [Archaea, Gram-positive, proteobacterial Gram-negative, singular bacteria, mitochondria, plastids, eukaryotic endoplasmic reticulum (ER) isoforms, eukaryotic cytoplasmic isoforms]. All group consensus comparisons tend to summarize better the alignments than do the individual sequence comparisons. The global individual consensus "matches" 87% with the consensus of consensuses sequence. A functional analysis of the global consensus identifies a (new) highly significant mixed charge cluster proximal to the carboxyl terminus of the sequence highlighting the hypercharge run EEDKKRRER (one-letter aa code used). The individual Archaea and Gram-positive sequences contain a corresponding significant mixed charge cluster in the location of the charge cluster of the consensus sequence. In contrast, the four Gram-negative proteobacterial sequences of the alignment do not have a charge cluster (even at the 5% significance level). All eukaryotic HSP70 sequences have the analogous charge cluster. Strikingly, several of the eukaryotic isoforms show multiple mixed charged clusters. These clusters were interpreted with supporting data related to HSP70 activity in facilitating chaperone, transport, and secretion function. We observed that the consensus contains only a single tryptophan residue and a single conserved cysteine. This is interpreted with respect to the target rule for disaggregating misfolded proteins. The mitochondrial HSP70 connections to bacterial HSP70 are analyzed, suggesting a polyphyletic split of Trypanosoma and Leishmania protist mitochondrial (Mt) homologues separated from Mt-animal/fungal/plant homologues. Moreover, the HSP70 sequences from the amitochondrial Entamoeba histolytica and Trichomonas vaginalis species were analyzed. The E. histolytica HSP70 is most similar to the higher eukaryotic cytoplasmic sequences, with significantly weaker alignments to ER sequences and much diminished matching to all eubacterial, mitochondrial, and chloroplast sequences. This appears to be at variance with the hypothesis that E. histolytica rather recently lost its mitochondrial organelle. T. vaginalis contains two HSP70 sequences, one Mt-like and the second similar to eukaryotic cytoplasmic sequences suggesting two diverse origins.     

Xrcc3 is required for assembly of Rad51 complexes in vivo

Rad51 is a member of a family of eukaryotic proteins related to the bacterial recombinational repair protein RecA. Rad51 protein localizes to multiple subnuclear foci in Chinese hamster ovary cells. Subnuclear Rad51 foci are induced by ionizing radiation or the DNA cross-linking agent cisplatin. Formation of these foci is likely to reflect assembly of a multimeric form of Rad51 that promotes DNA repair. Formation of damage-induced Rad51 foci does not occur in the Chinese hamster ovary cell line irs1SF, which is sensitive to DNA damaging agents. The Rad51 focus formation defect of irs1SF cells is corrected by a construct that encodes the repair protein Xrcc3. Xrcc3 is a human homolog of Rad51 previously isolated by virtue of its ability to correct the radiation sensitivity of irs1SF cells. Changes in the steady state level of Rad51 protein do not account for the irs1SF defect nor do they account for the appearance of foci following DNA damage. These results suggest that Xrcc3 is required for the assembly or stabilization of a multimeric form of Rad51 during DNA repair. Cell lines defective in two different components of DNA protein kinase formed Rad51 foci in response to damage, indicating DNA protein kinase is not required for damaged-induced mobilization of Rad51.     

Rad52 protein stimulates DNA strand exchange by Rad51 and replication protein A

The generation of a double-strand break in the Saccharomyces cerevisiae genome is a potentially catastrophic event that can induce cell-cycle arrest or ultimately result in loss of cell viability. The repair of such lesions is strongly dependent on proteins encoded by the RAD52 epistasis group of genes (RAD50-55, RAD57, MRE11, XRS2), as well as the RFA1 and RAD59 genes. rad52 mutants exhibit the most severe phenotypic defects in double-strand break repair, but almost nothing is known about the biochemical role of Rad52 protein. Rad51 protein promotes DNA strand exchange and acts similarly to RecA protein. Yeast Rad52 protein interacts with Rad51 protein, binds single-stranded DNA and stimulates annealing of complementary single-stranded DNA. We find that Rad52 protein stimulates DNA strand exchange by targeting Rad51 protein to a complex of replication protein A (RPA) with single-stranded DNA. Rad52 protein affects an early step in the reaction, presynaptic filament formation, by overcoming the inhibitory effects of the competitor, RPA. Furthermore, stimulation is dependent on the concerted action of both Rad51 protein and RPA, implying that specific protein-protein interactions between Rad52 protein, Rad51 protein and RPA are required.     

Prediction of homology and divergence in the secondary structure of polypeptides

A quantitative procedure is described for the comparison of secondary structure of homologous proteins. Standard predictive methods are used to generate probability profiles from pairs of homologous amino acid sequences; correlation coefficients (R) are then computed between each pair of amino acids for alpha-helix (R alpha), extended structure (R beta), turn (R(t)), and coil (R(c)). R values are >0.2 for correctly aligned homologous sequences. Unrelated or incorrectly aligned sequences give R values near zero. Lack of correlation for a segment of otherwise well-correlated sequences is used to identify structural divergence, which is then evaluated graphically by using difference profiles. A combination of these techniques correctly predicts secondary structural differences between melittin or beta-endorphin and their respective synthetic analogs. The method is potentially useful to describe evolutionary changes in protein secondary structure as well as in the design of peptide analogs.     

The DNA binding properties of Saccharomyces cerevisiae Rad51 protein

Saccharomyces cerevisiae Rad51 protein is the paradigm for eukaryotic ATP-dependent DNA strand exchange proteins. To explain some of the unique characteristics of DNA strand exchange promoted by Rad51 protein, when compared with its prokaryotic homologue the Escherichia coli RecA protein, we analyzed the DNA binding properties of the Rad51 protein. Rad51 protein binds both single-stranded DNA (ssDNA) and double-stranded DNA (dsDNA) in an ATP- and Mg2+-dependent manner, over a wide range of pH, with an apparent binding stoichiometry of approximately 1 protein monomer per 4 (+/-1) nucleotides or base pairs, respectively. Only dATP and adenosine 5'-gamma-(thiotriphosphate) (ATPgammaS) can substitute for ATP, but binding in the presence of ATPgammaS requires more than a 5-fold stoichiometric excess of protein. Without nucleotide cofactor, Rad51 protein binds both ssDNA and dsDNA but only at pH values lower than 6.8; in this case, the apparent binding stoichiometry covers the range of 1 protein monomer per 6-9 nucleotides or base pairs. Therefore, Rad51 protein displays two distinct modes of DNA binding. These binding modes are not inter-convertible; however, their initial selection is governed by ATP binding. On the basis of these DNA binding properties, we conclude that the main reason for the low efficiency of the DNA strand exchange promoted by Rad51 protein in vitro is its enhanced dsDNA-binding ability, which inhibits both the presynaptic and synaptic phases of the DNA strand exchange reaction as follows: during presynapsis, Rad51 protein interacts with and stabilizes secondary structures in ssDNA thereby inhibiting formation of a contiguous nucleoprotein filament; during synapsis, Rad51 protein inactivates the homologous dsDNA partner by directly binding to it.     

Structural polymorphism of the RecA protein from the thermophilic bacterium Thermus aquaticus

The Escherichia coli RecA protein has served as a model for understanding protein-catalyzed homologous recombination, both in vitro and in vivo. Although RecA proteins have now been sequenced from over 60 different bacteria, almost all of our structural knowledge about RecA has come from studies of the E. coli protein. We have used electron microscopy and image analysis to examine three different structures formed by the RecA protein from the thermophilic bacterium Thermus aquaticus. This protein has previously been shown to catalyze an in vitro strand exchange reaction at an optimal temperature of about 60 degrees C. We show that the active filament formed by the T. aquaticus RecA on DNA in the presence of a nucleotide cofactor is extremely similar to the filament formed by the E. coli protein, including the extension of DNA to a 5.1-A rise per base pair within this filament. This parameter appears highly conserved through evolution, as it has been observed for the eukaryotic RecA analogs as well. We have also characterized bundles of filaments formed by the T. aquaticus RecA in the absence of both DNA and nucleotide cofactor, as well as hexameric rings of the protein formed under all conditions examined. The bundles display a very large plasticity of mass within the RecA filament, as well as showing a polymorphism in filament-filament contacts that may be important to understanding mutations that affect surface residues on the RecA filament.     

Human Rad51 protein can form homologous joints in the absence of net strand exchange

The eukaryotic homologs of RecA protein are central enzymes of recombination and repair, and notwithstanding a high degree of conservation they differ sufficiently from RecA to offer insights into mechanisms and biological roles. The yield of DNA strand exchange reactions driven by both Escherichia coli RecA protein and its human homolog HsRad51 protein was inversely related to the GC content of oligonucleotide substrates, but at any given GC composition, HsRad51 promoted less exchange than RecA. When 40% of bases were GC pairs, the rate constant for strand exchange by HsRad51 was unmeasurable, whereas the rate constants for homologous pairing were unaltered relative to more AT-rich DNA. The ability of HsRad51 to form joints in the absence of net strand exchange was confirmed by experiments in which heterologous blocks at both ends of linear duplex oligonucleotides produced joints that instantly dissociated upon deproteinization. These findings suggest that HsRad51 acting alone on human DNA in vivo is a pairing protein that cannot form extensive heteroduplex DNA.     

Validity of the children and youth physical self-perception profile: a confirmatory factor analysis

MOTIVATION: The Match-Box software comprises protein sequence alignment tools based on strict statistical thresholds of similarity between protein segments. The method circumvents the gap penalty requirement: gaps being the result of the alignment and not a governing parameter of the procedure. The reliable conserved regions outlined by Match-Box are particularly relevant for homology modelling of protein structures, prediction of essential residues for site-directed mutagenesis and oligonucleotide design for cloning homologous genes by polymerase chain reaction (PCR). RESULTS: The method produces reliable results, as assessed by tests performed on protein families of known structures and of low sequence similarity. A reliability score is computed in relation to a threshold of similarity progressively raised to extend the aligned regions to their maximal length, up to the significance limit of matching segments. The score obtained at each position is printed below the sequences and allows a discriminant reading of each aligned region. AVAILABILITY: Sequences may be submitted to a Web server at http://www.fundp.ac.be/sciences/biologie/bms/+ ++matchbox_submit.html or sent by e-mail to matchbox/biq.fundp.ac.be (help available by just mailing help).     

Evolutionary relationships among the serpins

The serpins are a widely distributed group of serine proteinase inhibitors found in plants, birds, mammals and viruses. Despite the great evolutionary divergence of these organisms, their serpins are highly conserved, both in sequence and structurally. Amino acid sequences were aligned by a combination of automatic algorithms and by consideration of conserved structural elements in those serpins for which crystal structures exist. The program HOMED was used which allowed the alignment of amino acids to be simultaneously converted into the equivalently aligned nucleotide sequences. The aligned amino acids were used as the basis for superposition of the four known three-dimensional structures for which coordinates are available and compared with an optimal three-dimensional superposition in order to estimate the reliability of the sequence alignment. Phylogenetic relationships implied by these nucleotide sequence alignments were determined by the method of maximum parsimony. The proposed gene tree suggested that as much diversity existed between the plant serpin and mammalian serpins as was present among mammalian serpins and provided further evidence that the architecture of serpin molecules is highly constrained.     

Regulation of Rad51 function by c-Abl in response to DNA damage

The Rad51 protein, a homolog of bacterial RecA, functions in DNA double-strand break repair and genetic recombination. Whereas Rad51 catalyzes ATP-dependent pairing and strand exchange between homologous DNA molecules, regulation of this function is unknown. The c-Abl tyrosine kinase is activated by ionizing radiation and certain other DNA-damaging agents. Here we demonstrate that c-Abl interacts constitutively with Rad51. We show that c-Abl phosphorylates Rad51 on Tyr-54 in vitro. The results also show that treatment of cells with ionizing radiation induces c-Abl-dependent phosphorylation of Rad51. Phosphorylation of Rad51 by c-Abl inhibits the binding of Rad51 to DNA and the function of Rad51 in ATP-dependent DNA strand exchange reactions. These findings represent the first demonstration that Rad51 is regulated by phosphorylation and support a functional role for c-Abl in regulating Rad51-dependent recombination in the response to DNA damage.     

The human Rad51 protein: polarity of strand transfer and stimulation by hRP-A

The human Rad51 protein is homologous to the RecA protein and catalyses homologous pairing and strand transfer reactions in vitro. Using single-stranded circular and homologous linear duplex DNA, we show that hRad51 forms stable joint molecules by transfer of the 5' end of the complementary strand of the linear duplex to the ssDNA. The polarity of strand transfer is therefore 3' to 5', defined relative to the ssDNA on which hRad51 initiates filament formation. This polarity is opposite to that observed with RecA. Homologous pairing and strand transfer require stoichiometric amounts of hRad51, corresponding to one hRad51 monomer per three nucleotides of ssDNA. Joint molecules are not observed when the protein is present in limiting or excessive amounts. The human ssDNA binding-protein, hRP-A, stimulates hRad51-mediated reactions. Its effect is consistent with a role in the removal of secondary structures from ssDNA, thereby facilitating the formation of continuous Rad51 filaments.     

No sliding during homology search by RecA protein

The RecA protein of Escherichia coli is a prototype of the RecA/Rad51 family of proteins that exist in virtually all the organisms. In a process called DNA synapsis, RecA first polymerizes onto a single-stranded DNA (ssDNA) molecule; the resulting RecA-ssDNA complex then searches for and binds to a double-stranded DNA (dsDNA) molecule containing the almost identical, or "homologous, " sequence. The RecA-ssDNA complex thus can be envisioned as a sequence-specific binding entity. How does the complex search for its target buried within nonspecific sequences? One possible mechanism is the sliding mechanism, in which the complex first binds to a dsDNA molecule nonspecifically and then linearly diffuses, or slides, along the dsDNA. To understand the mechanism of homology search by RecA, this sliding model was tested. A plasmid containing four homologous targets in tandem was constructed and used as the dsDNA substrate in the synapsis reaction. If the sliding is the predominant search mode, the two outermost targets should act as more efficient targets than the inner targets. No such positional preference was observed, indicating that a long range sliding of the RecA-ssDNA complex does not occur. These and other available data can be adequately explained by a simple three-dimensional random collision mechanism.     

Bacterioopsin, haloopsin, and sensory opsin I of the halobacterial isolate Halobacterium sp. strain SG1: three new members of a growing family

The genes coding for bacterioopsin, haloopsin, and sensory opsin I of a halobacterial isolate from the Red Sea called Halobacterium sp. strain SG1 have been cloned and sequenced. The deduced protein sequences were aligned to the previously known halobacterial retinal proteins. The addition of these new sequences lowered the number of conserved residues to only 23 amino acids, or 8% of the alignment. Data base searches with two highly conserved peptides as well as with an alignment profile yielded no significant similarity to any other protein, so the halobacterial retinal proteins should be regarded as a distinct protein family. The protein alignment was used to make predictions about the structure of the retinal proteins as well as about the amino acids in contact with retinal proteins. These results were in excellent agreement with the structural model of bacteriorhodopsin of Halobacterium halobium as well as with mutant studies, indicating that (i) structure predictions based on the sequences of a membrane protein family can be quite accurate; (ii) halorhodopsin and sensory rhodopsin I have tertiary structures similar to that of bacteriorhodopsin; (iii) conserved amino acids do not take part in reactions specific for one group of proteins, e.g., proton translocation for bacteriorhodopsins, but have a crucial role in determining the conformation and reactions of the chromophore; and (iv) the general mode of action (light-induced chromophore and protein movements) is the same for all halobacterial retinal proteins, ion pumps as well as sensors.     

Domains required for dimerization of yeast Rad6 ubiquitin-conjugating enzyme and Rad18 DNA binding protein

The RAD6 gene of Saccharomyces cerevisiae encodes a ubiquitin-conjugating enzyme required for postreplicational repair of UV-damaged DNA and for damage-induced mutagenesis. In addition, Rad6 functions in the N end rule pathway of protein degradation. Rad6 mediates its DNA repair role via its association with Rad18, whose DNA binding activity may target the Rad6-Rad18 complex to damaged sites in DNA. In its role in N end-dependent protein degradation, Rad6 interacts with the UBR1-encoded ubiquitin protein ligase (E3) enzyme. Previous studies have indicated the involvement of N-terminal and C-terminal regions of Rad6 in interactions with Ubr1. Here, we identify the regions of Rad6 and Rad18 that are involved in the dimerization of these two proteins. We show that a region of 40 amino acids towards the C terminus of Rad18 (residues 371 to 410) is sufficient for interaction with Rad6. This region of Rad18 contains a number of nonpolar residues that have been conserved in helix-loop-helix motifs of other proteins. Our studies indicate the requirement for residues 141 to 149 at the C terminus, and suggest the involvement of residues 10 to 22 at the N terminus of Rad6, in the interaction with Rad18. Each of these regions of Rad6 is indicated to form an amphipathic helix.     

Automatic detection of conserved RNA structure elements in complete RNA virus genomes

We propose a new method for detecting conserved RNA secondary structures in a family of related RNA sequences. Our method is based on a combination of thermodynamic structure prediction and phylogenetic comparison. In contrast to purely phylogenetic methods, our algorithm can be used for small data sets of approximately 10 sequences, efficiently exploiting the information contained in the sequence variability. The procedure constructs a prediction only for those parts of sequences that are consistent with a single conserved structure. Our implementation produces reasonable consensus structures without user interference. As an example we have analysed the complete HIV-1 and hepatitis C virus (HCV) genomes as well as the small segment of hantavirus. Our method confirms the known structures in HIV-1 and predicts previously unknown conserved RNA secondary structures in HCV.