Biomedical subjects

M S Gelfand

Publications and source records attributed to M S Gelfand.

At least 19 recordsLinked to original sources

PubMed2006

Comparative genomics of regulation of heavy metal resistance in Eubacteria.

BACKGROUND: Heavy metal resistance (HMR) in Eubacteria is regulated by a variety of systems including transcription factors from the MerR family (COG0789). The HMR systems are characterized by the complex signal structure (strong palindrome within a 19 or 20 bp promoter spacer), and usually consist of transporter and regulator genes. Some HMR regulons also include detoxification systems. The number of sequenced bacterial genomes is constantly increasing and even though HMR resistance regulons of the COG0789 type usually consist of few genes per genome, the computational analysis may contribute to the understanding of the cellular systems of metal detoxification. RESULTS: We studied the mercury (MerR), copper (CueR and HmrR), cadmium (CadR), lead (PbrR), and zinc (ZntR) resistance systems and demonstrated that combining protein sequence analysis and analysis of DNA regulatory signals it was possible to distinguish metal-dependent members of COG0789, assign specificity towards particular metals to uncharacterized loci, and find new genes involved in the metal resistance, in particular, multicopper oxidase and copper chaperones, candidate cytochromes from the copper regulon, new cadmium transporters and, possibly, glutathione-S-transferases. CONCLUSION: Our data indicate that the specificity of the COG0789 systems can be determined combining phylogenetic analysis and identification of DNA regulatory sites. Taking into account signal structure, we can adequately identify genes that are activated using the DNA bending-unbending mechanism. In the case of regulon members that do not reside in single loci, analysis of potential regulatory sites could be crucial for the correct annotation and prediction of the specificity.

Bacteria↗

PubMed2006

Use of the flux model of amino acid metabolism of Escherichia coli.

A program implementing a flux model of Escherichia coli metabolism was used to analyze the effects of the addition of amino acids (tryptophan, tyrosine, phenylalanine, leucine, isoleucine, valine, histidine, lysine, threonine, cysteine, methionine, arginine, proline) to minimal medium or media lacking nitrogen, carbon, or both. The overall response of the metabolic system to the addition of various amino acids to the minimal medium is similar. Glycolysis and the synthesis of pyruvate with its subsequent degradation to acetate via acetyl-CoA become more efficient, whereas the fluxes through the pentose phosphate pathway and the TCA cycle decrease. If amino acids are used as the sole source of carbon, nitrogen, or both, the changes in the flux distribution are determined mainly by the carbon limitation. The phosphoenolpyruvate to glucose-6-phosphate flux increases; the flux through the pentose phosphate path is directed towards ribulose-5-phosphate. Other changes are determined by the compounds that are the primary products of catabolism of the added amino acid.

Amino Acids↗

PubMed2005

Alternative splicing and protein function.

BACKGROUND: Alternative splicing is a major mechanism of generating protein diversity in higher eukaryotes. Although at least half, and probably more, of mammalian genes are alternatively spliced, it was not clear, whether the frequency of alternative splicing is the same in different functional categories. The problem is obscured by uneven coverage of genes by ESTs and a large number of artifacts in the EST data. RESULTS: We have developed a method that generates possible mRNA isoforms for human genes contained in the EDAS database, taking into account the effects of nonsense-mediated decay and translation initiation rules, and a procedure for offsetting the effects of uneven EST coverage. Then we computed the number of mRNA isoforms for genes from different functional categories. Genes encoding ribosomal proteins and genes in the category "Small GTPase-mediated signal transduction" tend to have fewer isoforms than the average, whereas the genes in the category "DNA replication and chromosome cycle" have more isoforms than the average. Genes encoding proteins involved in protein-protein interactions tend to be alternatively spliced more often than genes encoding non-interacting proteins, although there is no significant difference in the number of isoforms of alternatively spliced genes. CONCLUSION: Filtering for functional isoforms satisfying biological constraints and accounting for uneven EST coverage allowed us to describe differences in alternative splicing of genes from different functional categories. The observations seem to be consistent with expectations based on current biological knowledge: less isoforms for ribosomal and signal transduction proteins, and more alternative splicing of interacting and cell cycle proteins.

Algorithms↗

PubMed2005

A Gibbs sampler for identification of symmetrically structured, spaced DNA motifs with improved estimation of the signal length.

MOTIVATION: Transcription regulatory protein factors often bind DNA as homo-dimers or hetero-dimers. Thus they recognize structured DNA motifs that are inverted or direct repeats or spaced motif pairs. However, these motifs are often difficult to identify owing to their high divergence. The motif structure included explicitly into the motif recognition algorithm improves recognition efficiency for highly divergent motifs as well as estimation of motif geometric parameters. RESULT: We present a modification of the Gibbs sampling motif extraction algorithm, SeSiMCMC (Sequence Similarities by Markov Chain Monte Carlo), which finds structured motifs of these types, as well as non-structured motifs, in a set of unaligned DNA sequences. It employs improved estimators of motif and spacer lengths. The probability that a sequence does not contain any motif is accounted for in a rigorous Bayesian manner. We have applied the algorithm to a set of upstream regions of genes from two Escherichia coli regulons involved in respiration. We have demonstrated that accounting for a symmetric motif structure allows the algorithm to identify weak motifs more accurately. In the examples studied, ArcA binding sites were demonstrated to have the structure of a direct spaced repeat, whereas NarP binding sites exhibited the palindromic structure. AVAILABILITY: The WWW interface of the program, its FreeBSD (4.0) and Windows 32 console executables are available at http://bioinform.genetika.ru/SeSiMCMC

Algorithms↗

PubMed2003

Experimental determination and system level analysis of essential genes in Escherichia coli MG1655.

Defining the gene products that play an essential role in an organism's functional repertoire is vital to understanding the system level organization of living cells. We used a genetic footprinting technique for a genome-wide assessment of genes required for robust aerobic growth of Escherichia coli in rich media. We identified 620 genes as essential and 3,126 genes as dispensable for growth under these conditions. Functional context analysis of these data allows individual functional assignments to be refined. Evolutionary context analysis demonstrates a significant tendency of essential E. coli genes to be preserved throughout the bacterial kingdom. Projection of these data over metabolic subsystems reveals topologic modules with essential and evolutionarily preserved enzymes with reduced capacity for error tolerance.

Aerobiosis↗

PubMed2003

Heat shock (sigma32 and HrcA/CIRCE) regulons in beta-, gamma- and epsilon-proteobacteria.

During heat shock, the main strategy of an organism is defense from denatured proteins. This is performed by chaperones that refold and proteases that cut abnormal proteins. In studying the sigma(32) and HrcA regulons in beta- and gamma-proteobacteria, we have found some new potential participants in the heat shock response and proposed the protein disulfide isomerase function for one of them. We describe the connection between the two regulons through cross-regulation of the HrcA repressor and sigma(32) in some beta-proteobacteria. Finally, we predict the binding signal for HrcA in epsilon-proteobacteria.

Amino Acid Sequence↗

PubMed2002

Finding weak motifs in DNA sequences.

Recognition of regulatory sites in unaligned DNA sequences is an old and well-studied problem in computational molecular biology. Recently, large-scale expression studies and comparative genomics brought this problem into a spotlight by generating a large number of samples with unknown regulatory signals. Here we develop algorithms for recognition of signals in corrupted samples (where only a fraction of sequences contain sites) with biased nucleotide composition. We further benchmark these and other algorithms on several bacterial and archaeal sites in a setting specifically designed to imitate the situations arising in comparative genomics studies.

Base Sequence↗

PubMed2001

Transcriptional regulation of pentose utilisation systems in the Bacillus/Clostridium group of bacteria.

In Bacillus subtilis, utilisation of xylose, arabinose and ribose is controlled by the transcriptional factors XylR, AraR and RbsR, respectively. Here we apply the comparative approach to the analysis of these regulons in the Bacillus/Clostridium group. Evolutionary variability of operon structures is demonstrated and operator sites for the main transcription factors are predicted. The consensus sequences for the XylR and RbsR binding sites vary in different subgroups of genomes. The functional coupling of gene clusters and the conservation of regulatory sites allow for detection of non-orthologous gene displacement of ribulose kinase in Enterococcus faecium and Clostridium acetobutylicum. Moreover, candidate catabolite responsive elements found upstream of most pentose-utilising genes suggest CcpA-mediated catabolite repression.

AraC Transcription Factor↗

PubMed2001

Computational analysis of the transcriptional regulation of pentose utilization systems in the gamma subdivision of Proteobacteria.

The comparative approach to the recognition of transcription regulatory sites is based on the assumption that as long as a regulator is conserved in several genomes, one can expect that sets of co-regulated genes (regulons) and regulatory sites for the regulator in these genomes are conserved as well. We used this approach to analyze the ribose (RbsR), arabinose (AraC), and xylose (XylR) regulons of gamma Proteobacteria for which (almost) completely sequenced genomes were available. Candidate binding sites for RbsR and AraC were detected. The improved XylR site consensus was proposed. Potential new members of the xylose regulons were found in the Escherichia coli, Salmonella typhi, and Klebsiella pneumoniae genomes. The function of these new xylose-regulated operons is likely to be the utilization of oligosaccharides containing xylose. Finally, candidate cAMP receptor-protein sites were identified in the regulatory regions of the majority of RbsR-, AraC-, and XylR-regulated operons.

AraC Transcription Factor↗

PubMed2001

Comparative analysis of FUR regulons in gamma-proteobacteria.

Iron is an essential element for the survival and pathogenesis of bacteria. The strict control of iron homeostasis is mediated by the FUR repressor, which is highly conserved among various bacterial species. Here we apply the comparative genomics approach to analyze candidate Fur-binding sites in the genomes of Escherichia coli (K12 and O157:H7), Salmonella typhi, Yersinia pestis and Vibrio cholerae. We describe a number of new loci encoding siderophore biosynthesis and transport proteins. A new regulator of iron-acquisition systems was found in S.typhi. We predict FUR regulation for several virulence systems. We also predict FUR regulation for the chemotaxis system of V.cholerae that is probably involved in the process of pathogenesis.

Amino Acid Sequence↗

PubMed2001

Computational analysis of candidate intron regulatory elements for tissue-specific alternative pre-mRNA splicing.

Alternative pre-mRNA splicing is a major cellular process by which functionally diverse proteins can be generated from the primary transcript of a single gene, often in tissue-specific patterns. The current study investigates the hypothesis that splicing of tissue-specific alternative exons is regulated in part by control sequences in adjacent introns and that such elements may be recognized via computational analysis of exons sharing a highly specific expression pattern. We have identified 25 brain-specific alternative cassette exons, compiled a dataset of genomic sequences encompassing these exons and their adjacent introns and used word contrast algorithms to analyze key features of these nucleotide sequences. By comparison to a control group of constitutive exons, brain-specific exons were often found to possess the following: divergent 5' splice sites; highly pyrimidine-rich upstream introns; a paucity of GGG motifs in the downstream intron; a highly statistically significant over-representation of the hexanucleotide UGCAUG in the proximal downstream intron. UGCAUG was also found at a high frequency downstream of a smaller group of muscle-specific exons. Intriguingly, UGCAUG has been identified previously in a few intron splicing enhancers. Our results indicate that this element plays a much wider role than previously appreciated in the regulated tissue-specific splicing of many alternative exons.

Algorithms↗

PubMed2001

Conservation of the binding site for the arginine repressor in all bacterial lineages.

BACKGROUND: The arginine repressor ArgR/AhrC is a transcription factor universally conserved in bacterial genomes. Its recognition signal (the ARG box), a weak palindrome, is also conserved between genomes, despite a very low degree of similarity between individual sites within a genome. Thus, the arginine repressor is different from two other universal transcription factors - HrcA, whose recognition signal is very strongly conserved both within and between genomes, and LexA/DinR, whose signal is strongly conserved within, but not between, genomes. The arginine regulon is well studied in Escherichia coli and to some extent in Bacillus subtilis and some other genomes. Here, we apply the comparative genomic approach to the prediction of the ArgR-binding sites in all completely sequenced bacterial genomes. RESULTS: Orthologs of ArgR/AhrC were identified in the complete genomes of E. coli, Haemophilus influenzae, Vibrio choleras, B. subtilis, Mycobacterium tuberculosis, Thermotoga maritima, Chlamydia pneumoniae and Deinococcus radiodurans. Candidate arginine repressor binding sites were identified upstream of arginine transport and metabolism genes. CONCLUSIONS: We found that the ArgR/AhrC recognition signal is conserved in all genomes that contain genes encoding orthologous transcription factors of this family. All genomes studied except M. tuberculosis contain ABC transport cassettes (related to the Art system of E. coli) belonging to the candidate arginine regulons.

Bacteria↗

PubMed2001

Pro-Frame: similarity-based gene recognition in eukaryotic DNA sequences with errors.

Performance of existing algorithms for similarity-based gene recognition in eukaryotes drops when the genomic DNA has been sequenced with errors. A modification of the spliced alignment algorithm allows for gene recognition in sequences with errors, in particular frameshifts. It tolerates up to 5% of sequencing errors without considerable drop of prediction reliability when a sufficiently close homologous protein is available (normalized evolutionary distance similarity score 50% or higher).

Algorithms↗

PubMed2001

Gene recognition in eukaryotic DNA by comparison of genomic sequences.

MOTIVATION: Sequencing of complete eukaryotic genomes and large syntenic fragments of genomes makes it possible to apply genomic comparison for gene recognition. RESULTS: This paper describes a spliced alignment algorithm that aligns candidate exon chains of two homologous genomic sequence fragments from different species. The algorithm is implemented in Pro-Gen software. Unlike other algorithms, Pro-Gen does not assume conservation of the exon-intron structure. Amino acid sequences obtained by the formal translation of candidate exons are aligned instead of nucleotide sequences, which allows for distant comparisons. The algorithm was tested on a sample of human-mammal (mouse), human-vertebrate (Xenopus ) and human-invertebrate (Drosophila ) gene pairs. Surprisingly, the best results, 97-98% correlation between the actual and predicted genes, were obtained for more distant comparisons, whereas the correlation on the human-mouse sample was only 93%. The latter value increases to 95% if conservation of the exon-intron structure is assumed. This is caused by a large amount of sequence conservation in non-coding regions of the human and mouse genes probably due to regulatory elements. AVAILABILITY: Pro-Gen v. 3.0 is available to academic researchers free of charge at http://www.anchorgen.com/pro_gen/pro_gen.html.

Algorithms↗

PubMed2001

Comparative approach to analysis of regulation in complete genomes: multidrug resistance systems in gamma-proteobacteria.

Comparative approach is a powerful tool for analysis of gene regulation in bacterial genomes. Here we apply it to analysis of regulation of the multidrug resistance transport (MDRT) systems in enterobacteria Escherichia coli, Salmonella typhi, Klebsiella pneumoniae and Yersinia pestis. Comparison of nucleotide sequences upstream of MDRT genes was performed in order to predict new regulatory sites (operators) and identify candidate regulons. Since the regulatory sites diverge slower than the non-coding regions in general, they are visible as strongly conserved islands. This analysis resulted in description of a regulatory network for known and hypothetical MDRT systems and porins. New candidate members of the MarA regulon were detected. Putative binding sites for EmrR and AcrR were suggested. A new hypothetical MarX regulon was described that includes some multidrug transporters and porins.

Bacterial Proteins↗

PubMed2001

Regulation of aromatic amino acid biosynthesis in gamma-proteobacteria.

Computational comparative techniques were applied to analysis of the aromatic amino acid regulons in gamma-proteobacteria. This resulted in characterization of the TrpR and TyrR regulons in the genomes of Yersinia pestis, Haemophilus influenzae, Vibrio cholerae and other bacteria and identification of new members of the PhhR regulon in the genome of Pseudomonas aeruginosa. Candidate attenuators were constructed for all studied genomes, including the trpBA operon of the very distantly related bacterium Chlamidia trachomatis. The pheA attenuator of Y. pestis is an integration site for the insertion element IS-200. It was shown that the triplication of the DAHP-synthase genes occurred prior to the divergence of families Enterobacteriaceae, Vibrionaceae and Alteromonadaceae. The candidate allosteric control site of the DAHP-syntheases was identified. This site is deteriorated in AroH of Buchnera sp. APS. The known DAHP-synthase of Bordetella pertussis is likely to be feedback-inhibited by phenylalanine, and the DAHP-synthase of Corynebacterium glutamicum could be inhibited by tyrosine. Overall, the most extensive regulation was observed in Escherichia coli, whereas the regulation in other genomes seems to be less developed. At the extreme, the tryptophan production in the aphid endosymbiont Buchnera sp. APS is free from transcriptional, attenuation, and allosteric control.

Allosteric Regulation↗

PubMed2000

Prediction of transcription regulatory sites in Archaea by a comparative genomic approach.

Intragenomic and intergenomic comparisons of upstream nucleotide sequences of archaeal genes were performed with the goal of predicting transcription regulatory sites (operators) and identifying likely regulons. Learning sets for the detection of regulatory sites were constructed using the available experimental data on archaeal transcription regulation or by analogy with known bacterial regulons, and further analysis was performed using iterative profile searches. The information content of the candidate signals detected by this method is insufficient for reliable predictions to be made. Therefore, this approach has to be complemented by examination of evolutionary conservation in different archaeal genomes. This combined strategy resulted in the prediction of a conserved heat shock regulon in all euryarchaea, a nitrogen fixation regulon in the methanogens Methanococcus jannaschii and Methanobacterium thermoautotrophicum and an aromatic amino acid regulon in M.thermoautotrophicum. Unexpectedly, the heat shock regulatory site was detected not only for genes that encode known chaperone proteins but also for archaeal histone genes. This suggests a possible function for archaeal histones in stress-related changes in DNA condensation. In addition, comparative analysis of the genomes of three Pyrococcus species resulted in the prediction of their purine metabolism and transport regulon. The results demonstrate the feasibility of prediction of at least some transcription regulatory sites by comparing poorly characterized prokaryotic genomes, particularly when several closely related genome sequences are available.

Archaea↗

PubMed2000

ASDB: database of alternatively spliced genes.

Version 2.1 of ASDB (Alternative Splicing Data Base) contains 1922 protein and 2486 DNA sequences. The protein entries from SWISS-PROT are joined into clusters corresponding to alternatively spliced variants of one gene. The DNA division consists of complete genes with alternative splicing mentioned or annotated in GenBank. The search engine allows one to search over SWISS-PROT and GenBank fields and then follow the links to all variants. The database can be assessed at the URL http://cbcg.nersc.gov/asdb

Alternative Splicing↗