Biomedical subjects

Svetlana A Shabalina

Publications and source records attributed to Svetlana A Shabalina.

16 recordsLinked to original sources

PubMed2006

Three major haplotypes of the beta2 adrenergic receptor define psychological profile, blood pressure, and the risk for development of a common musculoskeletal pain disorder.

Adrenergic receptor beta(2) (ADRB2) is a primary target for epinephrine. It plays a critical role in mediating physiological and psychological responses to environmental stressors. Thus, functional genetic variants of ADRB2 will be associated with a complex array of psychological and physiological phenotypes. These genetic variants should also interact with environmental factors such as physical or emotional stress to produce a phenotype vulnerable to pathological states. In this study, we determined whether common genetic variants of ADRB2 contribute to the development of a common chronic pain condition that is associated with increased levels of psychological distress and low blood pressure, factors which are strongly influenced by the adrenergic system. We genotyped 202 female subjects and examined the relationships between three major ADRB2 haplotypes and psychological factors, resting blood pressure, and the risk of developing a chronic musculoskeletal pain condition-Temporomandibular Joint Disorder (TMD). We propose that the first haplotype codes for lower levels of ADRB2 expression, the second haplotype codes for higher ADRB2 expression, and the third haplotype codes for higher receptor expression and rapid agonist-induced internalization. Individuals who carried one haplotype coding for high and one coding for low ADRB2 expression displayed the highest positive psychological traits, had higher levels of resting arterial pressure, and were about 10 times less likely to develop TMD. Thus, our data suggest that either positive or negative imbalances in ADRB2 function increase the vulnerability to chronic pain conditions such as TMD through different etiological pathways that imply the need for tailored treatment options.

Adolescent↗

PubMed2006

A periodic pattern of mRNA secondary structure created by the genetic code.

Single-stranded mRNA molecules form secondary structures through complementary self-interactions. Several hypotheses have been proposed on the relationship between the nucleotide sequence, encoded amino acid sequence and mRNA secondary structure. We performed the first transcriptome-wide in silico analysis of the human and mouse mRNA foldings and found a pronounced periodic pattern of nucleotide involvement in mRNA secondary structure. We show that this pattern is created by the structure of the genetic code, and the dinucleotide relative abundances are important for the maintenance of mRNA secondary structure. Although synonymous codon usage contributes to this pattern, it is intrinsic to the structure of the genetic code and manifests itself even in the absence of synonymous codon usage bias at the 4-fold degenerate sites. While all codon sites are important for the maintenance of mRNA secondary structure, degeneracy of the code allows regulation of stability and periodicity of mRNA secondary structure. We demonstrate that the third degenerate codon sites contribute most strongly to mRNA stability. These results convincingly support the hypothesis that redundancies in the genetic code allow transcripts to satisfy requirements for both protein structure and RNA structure. Our data show that selection may be operating on synonymous codons to maintain a more stable and ordered mRNA secondary structure, which is likely to be important for transcript stability and translation. We also demonstrate that functional domains of the mRNA [5'-untranslated region (5'-UTR), CDS and 3'-UTR] preferentially fold onto themselves, while the start codon and stop codon regions are characterized by relaxed secondary structures, which may facilitate initiation and termination of translation.

3' Untranslated Regions↗

PubMed2006

Analysis of internal loops within the RNA secondary structure in almost quadratic time.

MOTIVATION: Evaluating all possible internal loops is one of the key steps in predicting the optimal secondary structure of an RNA molecule. The best algorithm available runs in time O(L(3)), L is the length of the RNA. RESULTS: We propose a new algorithm for evaluating internal loops, its run-time is O(M(*)log(2)L), M < L(2) is a number of possible nucleotide pairings. We created a software tool Afold which predicts the optimal secondary structure of RNA molecules of lengths up to 28 000 nt, using a computer with 2 Gb RAM. We also propose algorithms constructing sets of conditionally optimal multi-branch loop free (MLF) structures, e.g. the set that for every possible pairing (x, y) contains an optimal MLF structure in which nucleotides x and y form a pair. All the algorithms have run-time O(M(*)log(2)L).

Algorithms↗

PubMed2006

A putative RNA-interference-based immune system in prokaryotes: computational analysis of the predicted enzymatic machinery, functional analogies with eukaryotic RNAi, and hypothetical mechanisms of action.

BACKGROUND: All archaeal and many bacterial genomes contain Clustered Regularly Interspaced Short Palindrome Repeats (CRISPR) and variable arrays of the CRISPR-associated (cas) genes that have been previously implicated in a novel form of DNA repair on the basis of comparative analysis of their protein product sequences. However, the proximity of CRISPR and cas genes strongly suggests that they have related functions which is hard to reconcile with the repair hypothesis. RESULTS: The protein sequences of the numerous cas gene products were classified into approximately 25 distinct protein families; several new functional and structural predictions are described. Comparative-genomic analysis of CRISPR and cas genes leads to the hypothesis that the CRISPR-Cas system (CASS) is a mechanism of defense against invading phages and plasmids that functions analogously to the eukaryotic RNA interference (RNAi) systems. Specific functional analogies are drawn between several components of CASS and proteins involved in eukaryotic RNAi, including the double-stranded RNA-specific helicase-nuclease (dicer), the endonuclease cleaving target mRNAs (slicer), and the RNA-dependent RNA polymerase. However, none of the CASS components is orthologous to its apparent eukaryotic functional counterpart. It is proposed that unique inserts of CRISPR, some of which are homologous to fragments of bacteriophage and plasmid genes, function as prokaryotic siRNAs (psiRNA), by base-pairing with the target mRNAs and promoting their degradation or translation shutdown. Specific hypothetical schemes are developed for the functioning of the predicted prokaryotic siRNA system and for the formation of new CRISPR units with unique inserts encoding psiRNA conferring immunity to the respective newly encountered phages or plasmids. The unique inserts in CRISPR show virtually no similarity even between closely related bacterial strains which suggests their rapid turnover, on evolutionary scale. Corollaries of this finding are that, even among closely related prokaryotes, the most commonly encountered phages and plasmids are different and/or that the dominant phages and plasmids turn over rapidly. CONCLUSION: We proposed previously that Cas proteins comprise a novel DNA repair system. The association of the cas genes with CRISPR and, especially, the presence, in CRISPR units, of unique inserts homologous to phage and plasmid genes make us abandon this hypothesis. It appears most likely that CASS is a prokaryotic system of defense against phages and plasmids that functions via the RNAi mechanism. The functioning of this system seems to involve integration of fragments of foreign genes into archaeal and bacterial chromosomes yielding heritable immunity to the respective agents. However, it appears that this inheritance is extremely unstable on the evolutionary scale such that the repertoires of unique psiRNAs are completely replaced even in closely related prokaryotes, presumably, in response to rapidly changing repertoires of dominant phages and plasmids.

Journal Article↗

PubMed2006

Computational models with thermodynamic and composition features improve siRNA design.

BACKGROUND: Small interfering RNAs (siRNAs) have become an important tool in cell and molecular biology. Reliable design of siRNA molecules is essential for the needs of large functional genomics projects. RESULTS: To improve the design of efficient siRNA molecules, we performed a comparative, thermodynamic and correlation analysis on a heterogeneous set of 653 siRNAs collected from the literature. We used this training set to select siRNA features and optimize computational models. We identified 18 parameters that correlate significantly with silencing efficiency. Some of these parameters characterize only the siRNA sequence, while others involve the whole mRNA. Most importantly, we derived an siRNA position-dependent consensus, and optimized the free-energy difference of the 5' and 3' terminal dinucleotides of the siRNA antisense strand. The position-dependent consensus is based on correlation and t-test analyses of the training set, and accounts for both significantly preferred and avoided nucleotides in all sequence positions. On the training set, the two parameters' correlation with silencing efficiency was 0.5 and 0.36, respectively. Among other features, a dinucleotide content index and the frequency of potential targets for siRNA in the mRNA added predictive power to our model (R = 0.55). We showed that our model is effective for predicting the efficiency of siRNAs at different concentrations. We optimized a neural network model on our training set using three parameters characterizing the siRNA sequence, and predicted efficiencies for the test siRNA dataset recently published by Novartis. On this validation set, the correlation coefficient between predicted and observed efficiency was 0.75. Using the same model, we performed a transcriptome-wide analysis of optimal siRNA targets for 22,600 human mRNAs. CONCLUSION: We demonstrated that the properties of the siRNAs themselves are essential for efficient RNA interference. The 5' ends of antisense strands of efficient siRNAs are U-rich and possess a content similarity to the pyrimidine-rich oligonucleotides interacting with the polypurine RNA tracks that are recognized by RNase H. The advantage of our method over similar methods is the small number of parameters. As a result, our method requires a much smaller training set to produce consistent results. Other mRNA features, though expensive to compute, can slightly improve our model.

Artificial Intelligence↗

PubMed2005

Persistence time of loss-of-function mutations at nonessential loci affecting eye color in Drosophila melanogaster.

Persistence time of a mutant allele, the expected number of generations before its elimination from the population, can be estimated as the ratio of the number of segregating mutations per individual over the mutation rate per generation. We screened two natural populations of Drosophila melanogaster for mutations causing clear-cut eye phenotypes and detected 25 mutant alleles, falling into 19 complementation groups, in 1164 haploid genomes, which implies 0.021 eye mutations/genome. The de novo haploid mutation rate for the same set of loci was estimated as 2 x 10(-4) in a 10-generation mutation-accumulation experiment. Thus, the average persistence time of all mutations causing clear-cut eye phenotypes is approximately 100 generations (95% confidence interval: 61-219). This estimate shows that the strength of selection against phenotypically drastic alleles of nonessential loci is close to that against recessive lethals. In both cases, deleterious alleles are apparently eliminated by selection against heterozygous individuals, which show no visible phenotypic differences from wild type.

Animals↗

PubMed2004

Genetic basis for individual variations in pain perception and the development of a chronic pain condition.

Pain sensitivity varies substantially among humans. A significant part of the human population develops chronic pain conditions that are characterized by heightened pain sensitivity. We identified three genetic variants (haplotypes) of the gene encoding catecholamine-O-methyltransferase (COMT) that we designated as low pain sensitivity (LPS), average pain sensitivity (APS) and high pain sensitivity (HPS). We show that these haplotypes encompass 96% of the human population, and five combinations of these haplotypes are strongly associated (P=0.0004) with variation in the sensitivity to experimental pain. The presence of even a single LPS haplotype diminishes, by as much as 2.3 times, the risk of developing myogenous temporomandibular joint disorder (TMD), a common musculoskeletal pain condition. The LPS haplotype produces much higher levels of COMT enzymatic activity when compared with the APS or HPS haplotypes. Inhibition of COMT in the rat results in a profound increase in pain sensitivity. Thus, COMT activity substantially influences pain sensitivity, and the three major haplotypes determine COMT activity in humans that inversely correlates with pain sensitivity and the risk of developing TMD.

Adolescent↗

PubMed2004

Identification of regions in multiple sequence alignments thermodynamically suitable for targeting by consensus oligonucleotides: application to HIV genome.

BACKGROUND: Computer programs for the generation of multiple sequence alignments such as "Clustal W" allow detection of regions that are most conserved among many sequence variants. However, even for regions that are equally conserved, their potential utility as hybridization targets varies. Mismatches in sequence variants are more disruptive in some duplexes than in others. Additionally, the propensity for self-interactions amongst oligonucleotides targeting conserved regions differs and the structure of target regions themselves can also influence hybridization efficiency. There is a need to develop software that will employ thermodynamic selection criteria for finding optimal hybridization targets in related sequences. RESULTS: A new scheme and new software for optimal detection of oligonucleotide hybridization targets common to families of aligned sequences is suggested and applied to aligned sequence variants of the complete HIV-1 genome. The scheme employs sequential filtering procedures with experimentally determined thermodynamic cut off points: 1) creation of a consensus sequence of RNA or DNA from aligned sequence variants with specification of the lengths of fragments to be used as oligonucleotide targets in the analyses; 2) selection of DNA oligonucleotides that have pairing potential, greater than a defined threshold, with all variants of aligned RNA sequences; 3) elimination of DNA oligonucleotides that have self-pairing potentials for intra- and inter-molecular interactions greater than defined thresholds. This scheme has been applied to the HIV-1 genome with experimentally determined thermodynamic cut off points. Theoretically optimal RNA target regions for consensus oligonucleotides were found. They can be further used for improvement of oligo-probe based HIV detection techniques. CONCLUSIONS: A selection scheme with thermodynamic thresholds and software is presented in this study. The package can be used for any purpose where there is a need to design optimal consensus oligonucleotides capable of interacting efficiently with hybridization targets common to families of aligned RNA or DNA sequences. Our thermodynamic approach can be helpful in designing consensus oligonucleotides with consistently high affinity to target variants in evolutionary related genes or genomes.

Base Sequence↗

PubMed2004

The mammalian transcriptome and the function of non-coding DNA sequences.

For decades, researchers have focused most of their attention on protein-coding genes and proteins. With the completion of the human and mouse genomes and the accumulation of data on the mammalian transcriptome, the focus now shifts to non-coding DNA sequences, RNA-coding genes and their transcripts. Many non-coding transcribed sequences are proving to have important regulatory roles, but the functions of the majority remain mysterious.

Animals↗

PubMed2004

Comparative analysis of orthologous eukaryotic mRNAs: potential hidden functional signals.

Sequencing of multiple, nearly complete eukaryotic genomes creates opportunities for detecting previously unnoticed, subtle functional signals in non-coding regions. A genome-wide comparative analysis of orthologous sets of mammalian and yeast mRNAs revealed distinct patterns of evolutionary conservation at the boundaries of the untranslated regions (UTRs) and the coding region (CDS). Elevated sequence conservation was detected in approximately 30 nt regions around the start codon. There seems to be a complementary relationship between sequence conservation in the approximately 30 nt regions of the 5'-UTR immediately upstream of the start codon and that in the synonymous positions of the 5'-terminal 30 nt of the CDS: in mammalian mRNAs, the 5'-UTR shows a greater conservation than the CDS, whereas the opposite trend holds for yeast mRNAs. Unexpectedly, a approximately 30 nt region downstream of the stop codon shows a substantially lower level of sequence conservation than the downstream portions of the 3'-UTRs. However, the sequence in this poorly conserved 30 nt portion of the 3'-UTR is non-random in that it has a higher GC content than the rest of the UTR. It is hypothesized that the elevated sequence conservation in the region immediately upstream of the start codon is related to the requirement for initiation factor binding during pre-initiation ribosomal scanning. In contrast, the poorly conserved region downstream of the stop codon could be involved in the post- termination scanning and dissociation of the ribosomes from the mRNA, which requires only the mRNA-ribosome interaction. Additionally, it was found that the choice of the stop codon in mammals, but not in yeasts, and the context in the immediate vicinity of the stop codons in both mammals and yeasts are subject to strong selection. Thus, genome-wide analysis of orthologous gene sets allows detection of previously unrecognized patterns of sequence conservation, which are likely to reflect hidden functional signals, such as ribosomal filters that could regulate translation by modulating the interaction between the mRNA and ribosomes.

3' Untranslated Regions↗

PubMed2003

Patterns in interspecies similarity correlate with nucleotide composition in mammalian 3'UTRs.

Post-transcriptional regulation and the formation of mRNA 3' ends are crucial for gene expression in eukaryotes. Interspecies conservation of many sequences within 3'UTRs reveals selective constraint due to similar function. To study the pattern of conservation within 3'UTRs, we compiled and aligned 50 sets of complete orthologous 3'UTRs from four orders of mammals. We observed a mosaic pattern of conservation, with alternating regions of high (phylogenetic footprints) and low similarity. Conservation in 3'UTRs correlates with their base composition and also with the synonymous substitution rate in corresponding coding regions. The non-uniform distribution of conservation is more pronounced for 3'UTRs with a moderate or low level of overall conservation, where invariant nucleotides are more numerous, and their runs of lengths 4-7 occur more frequently than if conservation were random. Many runs of invariant nucleotides are AU-rich or pyrimidine-rich. Some of these runs coincide with known functional cis- elements of eukaryotic mRNAs, such as the U-rich upstream element, polyadenylation signal and DICE regulatory signal. More divergent regions of multiple alignments of 3'UTRs are often more G- and/or C-rich. Our results provide evidence on the importance of moderately conserved regions in 3'UTRs and suggest that regulatory functions of 3'UTRs might utilize gene-specific information in these regions.

3' Untranslated Regions↗

PubMed2003

A significant fraction of conserved noncoding DNA in human and mouse consists of predicted matrix attachment regions.

Noncoding DNA in the human-mouse orthologous intergenic regions contains "islands" of conserved sequences, the functions of which remain largely unknown. We hypothesized that some of these regions might be matrix-scaffold attachment regions, MARs (or S/MARs). MARs comprise one of the few classes of eukaryotic noncoding DNA with an experimentally characterized function, being involved in the attachment of chromatin to the nuclear matrix, chromatin remodeling and transcription regulation. To test our hypothesis, we analyzed the co-occurrence of predicted MARs with highly conserved noncoding DNA regions in human-mouse genomic alignments. We found that 11% of the conserved noncoding DNA consists of predicted MARs. Conversely, more than half of the predicted MARs co-occur with one or more independently identified conserved sequence blocks. An excess of conserved predicted MARs is seen in intergenic regions preceding 5' ends of genes, suggesting that these MARs are primarily involved in transcriptional control.

Animals↗

PubMed2002

Classification of common conserved sequences in mammalian intergenic regions.

Comparisons between orthologous intergenic regions of related genomes reveal numerous hits, i.e. pairs of relatively short highly similar sequences that evolved slowly, perhaps due to selective constraint. We analyzed and classified 2638 hits found within 100 pairs of complete, orthologous intergenic regions of human and murine genomes. We identified all common fragments of hits that align well with many other hits and constructed their classification. Our analysis revealed 20 abundant classes each containing 10 or more fragments. Fragments of the same class may perform the same function, e.g. bind a particular protein. Ten of the abundant classes apparently correspond to known functional consensuses, whereas others may represent novel conserved sites. Thus, large-scale comparative analysis of slowly evolving intergenic sequences can provide valuable insights into their function.

Animals↗

PubMed2002

Analysis of similarity within 142 pairs of orthologous intergenic regions of Caenorhabditis elegans and Caenorhabditis briggsae.

Patterns of similarity between genomes of related species reflect the distribution of selective constraint within DNA. We analyzed alignments of 142 orthologous intergenic regions of Caenorhabditis elegans and Caenorhabditis briggsae and found a mosaic pattern with regions of high similarity (phylogenetic footprints) interspersed with non-alignable sequences. Footprints cover approximately 20% of intergenic regions, often occur in clumps and are rare within 5' UTRs but common within 3' UTRs. The footprints have a higher ratio of transitions to transversions than expected at random and a higher GC content than the rest of the intergenic region. The number of footprints and the GC content of footprints within an intergenic region are higher when genes are oriented so that their 5' ends form the boundaries of the intergenic region. Overall, the patterns and characteristics identified here, along with other comparative and experimental studies, suggest that many footprints have a regulatory function, although other types of function are also possible. These conclusions may be quite general across eukaryotes, and the characteristics of conserved regulatory elements determined from genomic comparisons can be useful in prediction of regulation sites within individual DNA sequences.

Animals↗

PubMed2002

A hierarchical approach to aligning collinear regions of genomes.

MOTIVATION: As a first approximation, similarity between two long orthologous regions of genomes can be represented by a chain of local similarities. Within such a chain, pairs of successive similarities are collinear (non-conflicting), i.e. segments involved in the nth similarity precede in both sequences segments involved in the (n+1)th similarity. However, when all similarities between two long sequences are considered, usually there are many conflicts between them. Although some conflicts can be avoided by masking transposons or low-complexity sequences, selecting only those similarities that reflect orthology and, thus, belong to the evolutionarily true chain is not trivial. RESULTS: We propose a simple, hierarchical algorithm of finding the true chain of local similarities. Starting from similarities with low P-values, we resolve each pairwise conflict by deleting a similarity with a higher P-value. This greedy approach constructs a chain of similarities faster than when a chain optimal with respect to some global criterion is sought, and makes more sense biologically.

Algorithms↗

PubMed2002

OWEN: aligning long collinear regions of genomes.

OWEN is an interactive tool for aligning two long DNA sequences that represents similarity between them by a chain of collinear local similarities. OWEN employs several methods for constructing and editing local similarities and for resolving conflicts between them. Alignments of sequences of lengths over 10(6) can often be produced in minutes. OWEN requires memory below 20 L, where L is the sum of lengths of the compared sequences.

Algorithms↗