Side-chain rotamers in protein alpha-alpha-hairpins and a mechanism of their selection

The observed features of side-chain rotamer distributions in protein alpha-alpha-hairpins are described. It was found that in left-turned alpha-alpha-hairpins most side chains occupying d-positions have t-rotamers and those in g-positions g- -rotamers. In right-turned alpha-alpha-hairpins, most side chains in a-positions adopt g- -rotamers and those in e-positions t-rotamers. Analysis of these features enables us to conclude that selection of side-chain rotamers in alpha-alpha-hairpins depends on both the type of the alpha-helix packing and the residue position. The observed features can be explained taking into account the squeezing mechanism according to which interhelical interactions bring alpha-helices closer to each other and this effect squeezes side chains out of the helix-helix interface and as a result they adopt unique conformations.     

Structure of alpha-spiral hairpins with short connections in globular proteins

The analysis of conformations of more than 100 alpha-alpha-hairpins with closely packed helical segments and connections up to four amino acid residues in length was carried out. Five types of the connections were revealed and their phi and psi values on the Ramachandran map were found. Each type of alpha-alpha-hairpins was shown to have a unique sequence pattern for hydrophobic and hydrophilic residues.     

Mechanism of selection of side-chain rotamers in α-helices

In this study, a possible mechanism of selection of side-chain rotamers based on the rotamer distributions in known coiled-coil proteins is suggested. According to this mechanism, interhelical hydrophobic, polar, and packing interactions bring alpha-helices closer to each other and this effect squeezes side chains out of the helix-helix interface. As a result, in dimeric coiled coils and long alpha-alpha-hairpins where alpha-helices are packed in a face-to-face manner, most side chains occupying the a-positions have t-rotamers and those in the d-positions g(-)-rotamers. In tetramers, where alpha-helices are packed side-by-side, most side chains in the a-positions adopt g(-)-rotamers and those in the d-positions t-rotamers.     

From structure and dynamics of protein L7/L12 to molecular switching in ribosome

Based on the (1)H-(15)N NMR spectroscopy data, the three-dimensional structure and internal dynamic properties of ribosomal protein L7 from Escherichia coli were derived. The structure of L7 dimer in solution can be described as a set of three distinct domains, tumbling rather independently and linked via flexible hinge regions. The dimeric N-terminal domain (residues 1-32) consists of two antiparallel alpha-alpha-hairpins forming a symmetrical four-helical bundle, whereas the two identical C-terminal domains (residues 52-120) adopt a compact alpha/beta-fold. There is an indirect evidence of the existence of transitory helical structures at least in the first part (residues 33-43) of the hinge region. Combining structural data for the ribosomal protein L7/L12 from NMR spectroscopy and x-ray crystallography, it was suggested that its hinge region acts as a molecular switch, initiating "ratchet-like" motions of the L7/L12 stalk with respect to the ribosomal surface in response to elongation factor binding and GTP hydrolysis. This hypothesis allows an explanation of events observed during the translation cycle and provides useful insights into the role of protein L7/L12 in the functioning of the ribosome.     

5S rRNA sequences of myxobacteria and radioresistant bacteria and implications for eubacterial evolution

5S rRNA sequences were determined for the myxobacteria Cystobacter fuscus, Myxococcus coralloides, Sorangium cellulosum, and Nannocystis exedens and for the radioresistant bacteria Deinococcus radiodurans and Deinococcus radiophilus. A dendrogram was constructed by using weighted pairwise grouping based on these and all other previously known eubacterial 5S rRNA sequences, and this dendrogram showed differences as well as similarities compared with results derived from 16S rRNA analyses. In the dendrogram, Deinococcus 5S rRNA sequences clustered with 5S rRNA sequences of the genus Thermus, as suggested by the results of 16S rRNA analyses. However, in contrast to the 16S rRNA results, the Deinococcus-Thermus cluster divided the 5S rRNA sequences of the alpha subdivision of the class Proteobacteria from the 5S rRNA sequences of the beta and gamma subgroups of the Proteobacteria. The myxobacterial 5S rRNA sequence data failed to confirm the existence of a delta subgroup of the class Proteobacteria, which was suggested by the results of 16S rRNA analyses.     

Recombination as a Point Process along Sequences

Histories of sequences in the coalescent model with recombination can be simulated using an algorithm that takes as input a sample of extant sequences. The algorithm traces the history of the sequences going back in time, encountering recombinations and coalescence (duplications) until the ancestral material is located on one sequence for homologous positions in the present sequences. Here an alternative algorithm is formulated not as going back in time and operating on sequences, but by moving spatially along the sequences, updating the history of the sequences as recombination points are encountered. This algorithm focuses on spatial aspects of the coalescent with recombination rather than on temporal aspects as is the case of familiar algorithms. Mathematical results related to spatial aspects of the coalescent with recombination are derived.     

Comparison of mealybug (Planococcus lilacinus) and fruit fly genomes: isolation and analysis of conserved sequences and their utility in studying synteny in the mealybug

By using ligation-mediated PCR products from mealybug DNA as tester and biotinylated fly DNA as driver, we recovered a fraction of the tester that remains hybridized to driver following high-stringency washing conditions. This fraction is expected to contain mealybug sequences conserved in the fly (MCF). Reciprocal experiments enabled the isolation of fly sequences conserved in the mealybug (FCM). Coding sequences among MCF show amino acid identities >40% with fly proteins, allowing a reliable identification of orthologs. Three sequences from the fly cytogenetic positions 98-99 were hybridized onto mealybug chromosomes and the results identified differences in synteny between the two species. Taken together, our results present a method for direct isolation of sequences conserved between an 'orphan' (mealybug) genome and a 'reference' (fly) genome and showed that these sequences can be used to study chromosome synteny in the mealybug.     

Monophyly of beta-tubulin and H+-ATPase gene variants in Glomus intraradices: consequences for molecular evolutionary studies of AM fungal genes

Arbuscular mycorrhizal fungi (AMF) are an ecologically important group of fungi. Previous studies showed the presence of divergent copies of beta-tubulin and V-type vacuolar H+-ATPase genes in AMF genomes and suggested horizontal gene transfer from host plants or mycoparasites to AMF. We sequenced these genes from DNA isolated from an in vitro cultured isolate of Glomus intraradices that was free of any obvious contaminants. We found two highly variable beta-tubulin sequences and variable H+-ATPase sequences. Despite this high variation, comparison of the sequences with those in gene banks supported a glomeromycotan origin of G. intraradices beta-tubulin and H+-ATPase sequences. Thus, our results are in sharp contrast with the previously reported polyphyletic origin of those genes. We present evidence that some highly divergent sequences of beta-tubulin and H+-ATPase deposited in the databases are likely to be contaminants. We therefore reject the prediction of horizontal transfer to AMF genomes. High differences in GC content between glomeromycotan sequences and sequences grouping in other lineages are shown and we suggest they can be used as an indicator to detect such contaminants. H+-ATPase phylogeny gave unexpected results and failed to resolve fungi as a natural group. beta-Tubulin phylogeny supported Glomeromeromycota as sister group of the Chytridiomycota. Contrasts between our results and trees previously generated using rDNA sequences are discussed.     

Soybean Recombination Sites are Present as Dispersed Segments in Arabidopsis and Liverwort Mitochondrial DNA

In previous studies we have shown that recombination across a 299-bp interspersed sequence accounts for the diversity of the mitochondrial genome in wild and cultivated soybeans. In this study, a computer-assisted survey of databases was performed using sequences of the repeat and its neighboring regions as query sequences. The sequences of soybean were found to be present as many short segments that include repeated sequences in the mitochondrial genomes of Arabidopsis and liverwort. Taken together with the results of a DNA gel-blot analysis, this suggests that the soybean sequences were found to have originated during land plant evolution and are present as small-interspersed segments in many taxa of land plants.     

Enhancer-like structures in moderately repetitive sequences of eukaryotic genomes

The results of contextual analysis of 25 different middle repetitive DNA sequences are presented. It was shown that each of these repetitive DNA sequences contains at least one enhancer-like structure homologous to real enhancers, as well as to their consensus. The enhancer-like structures have been also revealed in the replication origin of some prokaryote genomes. The results are discussed in the light of a possible role of middle repetitive DNA sequences in the modulation of gene expression. Some aspects of genomes' evolution, in relation to enhancers, are also considered.     

The nuclease sensitivity of active genes.

Brief micrococcal nuclease digestion of chick embryonic red blood cells results in preferential excision and solubilization of monomer nucleosomes associated with beta-globin sequences and also 5'-sequences flanking the beta-globin gene. Both regions are DNAse-I sensitive in nuclei. Such salt-soluble nucleosomes are enriched in all four major HMG proteins but HMG1 and 2 are only weakly associated. These nucleosomes appear to have lost much of the DNAse-I sensitivity of active genes. The HMG14 and 17-containing salt-soluble nucleosomes separated by electrophoresis are not DNAse-I sensitive and contain inactive gene sequences as well as active sequences. Reconstitution of HMG proteins onto bulk nucleosomes or chromatin failed to reveal an HMG-dependent sensitivity of active genes as assayed by dot-blot hybridization and it was found that the DNAse-I sensitivity of ASV proviral sequences as assayed by dot-blot hybridization was not HMG-dependent. These results indicate that higher order chromatin structures might be responsible for nuclease sensitivity of active genes.     

Nucleotide sequences of three different isoforms of beta-tubulin cDNA from Chinese hamster ovary cells.

The complete cDNA sequences of two clones encoding beta-tubulin isotypes and the partial sequence of a third isoform from Chinese hamster ovary cells have been determined. The deduced amino acid sequences of the three isoforms show extensive homology to each other as well as with other alpha and beta-tubulin sequences from various species. These results provide evidence for the expression of three different isoforms of beta-tubulin in Chinese hamster ovary cells.     

Effect of training datasets on support vector machine prediction of protein-protein interactions

Knowledge of protein-protein interaction is useful for elucidating protein function via the concept of 'guilt-by-association'. A statistical learning method, Support Vector Machine (SVM), has recently been explored for the prediction of protein-protein interactions using artificial shuffled sequences as hypothetical noninteracting proteins and it has shown promising results (Bock, J. R., Gough, D. A., Bioinformatics 2001, 17, 455-460). It remains unclear however, how the prediction accuracy is affected if real protein sequences are used to represent noninteracting proteins. In this work, this effect is assessed by comparison of the results derived from the use of real protein sequences with that derived from the use of shuffled sequences. The real protein sequences of hypothetical noninteracting proteins are generated from an exclusion analysis in combination with subcellular localization information of interacting proteins found in the Database of Interacting Proteins. Prediction accuracy using real protein sequences is 76.9% compared to 94.1% using artificial shuffled sequences. The discrepancy likely arises from the expected higher level of difficulty for separating two sets of real protein sequences than that for separating a set of real protein sequences from a set of artificial sequences. The use of real protein sequences for training a SVM classification system is expected to give better prediction results in practical cases. This is tested by using both SVM systems for predicting putative protein partners of a set of thioredoxin related proteins. The prediction results are consistent with observations, suggesting that real sequence is more practically useful in development of SVM classification system for facilitating protein-protein interaction prediction.     

The evolutionary stability of cytochrome c-551 in Pseudomonas aeruginosa and Pseudomonas fluorescens biotype C

Cytochrome c-551 was prepared from nine different strains of Pseudomonas aeruginosa and six of Pseudomonas fluorescens biotype C, and their amino acid sequences were compared with the sequences previously determined for the cytochromes of type strains of each species. The standard of sequence examination was such that all single amino acid substitutions, delections or insertions ought to have been detected. Balanced double changes in sites in the same part of the sequence might have escaped detection. The standard of some of the quantitative amino acid analyses was not as high as would be required for the investigation of completely unknown sequences. Eight of the Ps. aeruginosa sequences could not be distinguished from the type sequence, whereas the ninth had a single amino acid substitution. The sequences from Ps. fluorescens biotype C were more varied, differing in from zero to four substitutions from the type sequence, with the most diverse sequences differing in seven positions. The results for Ps. aeruginosa are interpreted as evidence that neutral mutations are not responsible for much molecular evolution. The superficially paradoxical differences in the results for the two species are discussed.     

A computer method for finding common base paired helices in aligned sequences: application to the analysis of random sequences.

We describe a new computer program that identifies conserved secondary structures in aligned nucleotide sequences of related single-stranded RNAs. The program employs a series of hash tables to identify and sort common base paired helices that are located in identical positions in more than one sequence. The program gives information on the total number of base paired helices that are conserved between related sequences and provides detailed information about common helices that have a minimum of one or more compensating base changes. The program is useful in the analysis of large biological sequences. We have used it to examine the number and type of complementary segments (potential base paired helices) that can be found in common among related random sequences similar in base composition to 16S rRNA from Escherichia coli. Two types of random sequences were analyzed. One set consisted of sequences that were independent but they had the same mononucleotide composition as the 16S rRNA. The second set contained sequences that were 80% similar to one another. Different results were obtained in the analysis of these two types of random sequences. When 5 sequences that were 80% similar to one another were analyzed, significant numbers of potential helices with two or more independent base changes were observed. When 5 independent sequences were analyzed, no potential helices were found in common. The results of the analyses with random sequences were compared with the number and type of helices found in the phylogenetic model of the secondary structure of 16S ribosomal RNA. Many more helices are conserved among the ribosomal sequences than are found in common among similar random sequences. In addition, conserved helices in the 16S rRNAs are, on the average, longer than the complementary segments that are found in comparable random sequences. The significance of these results and their application in the analysis of long non-ribosomal nucleotide sequences is discussed.     

Distribution pattern and enzymic hypermethylation of inverted repetitive DNA sequences in P815 mastocytoma cells.

A specific class of DNA sequences, the inverted repetitive sequences, forms a double-stranded structure within a single linear polynucleotide chain in denatured DNA. The reassociation process is unimolecular and occurs very fast. Quantitative analyses have shown that in mouse P815 cells these sequences comprise about 4% of the nuclear DNA and are interspersed within sequences of other degrees of repetitiveness. After labeling the cells with L-[Me-3H]methionine and [14C]deoxycytidine, relative rates of enzymic DNA methylation were computed on the basis of radioactivities found in pyrimidine residues of the nuclear DNA. The results indicate that in P815 cells, DNA of inverted repetitive sequences is methylated to a level about 50% higher than the normal repetitive DNA sequences and to about 300% higher than the unique and intermediary intermediatry sequences. The biological function of the inverted repetitive sequences, as well as of the role of enzymic methylation of DNA remains unknown.     

Rarity associated with specific ecological niches in the bacterial world: the 'Synergistes' example

The 'Synergistes' group, which apparently represents an as yet unnamed division of the bacteria, was explored in 93 anaerobic environments (guts, soils, digestors, etc.). From 16S rDNA gene-targeted polymerase chain reaction (PCR) assays, this group appeared to be present in 90% of the anaerobic microbial ecosystems analysed. The phylogeny of 103 16S rDNA sequences from 30 ecosystems showed a strong link between 16S rDNA sequences and given ecosystems. 'Synergistes' 16S rDNA sequences from animal sources (termites, guinea pigs, pigs, birds, etc.) formed clustered phylogenetical groups. 'Synergistes' groups were also associated either with anaerobic digestors and soils or with thermophilic conditions. Sequences available from the DNA database were consistent with the results. These results show the wide diversity of the 'Synergistes' division as well as the specific ecological niche of each 16S rDNA sequences.     

Structural organization of the genome of the cellular slime mold Dictyostelium discoideum: interspersion of repetitive and single-copy DNA sequences.

The length and interspersion of reiterated and single-copy DNA sequences in Dictyostelium have been examined. The results indicate that approximately 50-60% of the single-copy sequences in DNA fragments 1500 nucleotides long and 75% of the single-copy sequences in fragments 3000 nucleotides long are linked to short interspersed repeat DNA sequences. The average length of these single-copy sequences is 1500 nucleotides. The length of the reiterated DNA has also been analyzed and shows a bimodal distribution. One half is present in sequences greater than 2000 nucleotides long, while the remainder is present as short fragments 250-450 nucleotides long. These shorter fragments are interspersed with the bulk of the single-copy DNA.