Search

Translation termination is a finishing step of protein biosynthesis. The significant role in this process belongs not only to protein factors of translation termination but also to the nearest nucleotide environment of stop codons. There are numerous descriptions of stop codons readthrough, which is due to specific nucleotide sequences behind them. However, represented data are segmental and don't explain the mechanism of the nucleotide context influence on translation termination. It is well known that stop codon UAA usage is preferential for A/T-rich genes, and UAG, UGA-for G/C-rich genes, which is related to an expression level of these genes. We investigated the connection between a frequency of nucleotides occurrence in 3' area of stop codons in the human genome and their influence on translation termination efficiency. We found that 3' context motif, which is cognate to the sequence of a stop codon, stimulates translation termination. At the same time, the nucleotide composition of 3' sequence that differs from stop codon, decreases translation termination efficiency.

The genetic code is considered to use five nucleic bases (adenine, guanine, cytosine, thymine and uracil), which form two pairs for encoding information in DNA and two pairs for encoding information in RNA. Nevertheless, in recent years several artificial base pairs have been developed in attempts to expand the genetic code. Employment of these additional base pairs increases the information capacity and variety of DNA sequences, and provides a platform for the site-specific, enzymatic incorporation of extra functional components into DNA and RNA. As a result, of the development of such expanded systems, many artificial base pairs have been synthesized and tested under various conditions. Following many stages of enhancement, unnatural base pairs have been modified to eliminate their weak points, qualifying them for specific research needs. Moreover, the first attempts to create a semi-synthetic organism containing DNA with unnatural base pairs seem to have been successful. This further extends the possible applications of these kinds of pairs. Herein, we describe the most significant qualities of unnatural base pairs and their actual applications.

Gammaproteobacteria get energy for their growth from different carbon sources using either glycolysis or alternative metabolic pathways induced in stress conditions. These metabolic switches are coordinated by complex interplay of regulatory proteins sensing concentrations of available metabolites by mechanisms yet to be understood. Here, we use two transcriptional regulators, ExuR and UxuR, controlling d-galacturonate (d-gal) and d-glucuronate metabolism in Escherichia coli, as the targets for computational search of low-molecular compounds capable to bind their ligand-binding domains. Using a flexible molecular docking, we modeled the interactions of these proteins with substrates and intermediates of glycolysis, Ashwell and Entner-Doudoroff pathways. For UxuR, the two preferred sites of ligand binding were found: one is located within the C-terminal domain, while another occupies the interdomain space. For ExuR, the only one preferred site was detected in the interdomain area. Availability of this area to different ligands suggests that, similar to the Lac repressor, the DNA-binding properties of UxuR and ExuR may be changed by repositioning of their domains. Experimental assays confirmed the ability of ligands with highest affinities to bind the regulatory proteins and affect their interaction with DNA. d-gal that is carried into the cell by the ExuT transporter appeared to be the best ligand for repressor of the exuT transcription, ExuR. For UxuR, the highest affinity was found for d-fructuronate transported by GntP, which biosynthesis is repressed by UxuR. Providing a feedback loop to balance the concentrations of different nutrients, such ligand-mediated modulation can also coordinate switching between different metabolic pathways in bacteria.

Fitness landscapes depict how genotypes manifest at the phenotypic level and form the basis of our understanding of many areas of biology, yet their properties remain elusive. Previous studies have analysed specific genes, often using their function as a proxy for fitness, experimentally assessing the effect on function of single mutations and their combinations in a specific sequence or in different sequences. However, systematic high-throughput studies of the local fitness landscape of an entire protein have not yet been reported. Here we visualize an extensive region of the local fitness landscape of the green fluorescent protein from Aequorea victoria (avGFP) by measuring the native function (fluorescence) of tens of thousands of derivative genotypes of avGFP. We show that the fitness landscape of avGFP is narrow, with 3/4 of the derivatives with a single mutation showing reduced fluorescence and half of the derivatives with four mutations being completely non-fluorescent. The narrowness is enhanced by epistasis, which was detected in up to 30% of genotypes with multiple mutations and mostly occurred through the cumulative effect of slightly deleterious mutations causing a threshold-like decrease in protein stability and a concomitant loss of fluorescence. A model of orthologous sequence divergence spanning hundreds of millions of years predicted the extent of epistasis in our data, indicating congruence between the fitness landscape properties at the local and global scales. The characterization of the local fitness landscape of avGFP has important implications for several fields including molecular evolution, population genetics and protein design., Competing Interests: The authors declare no interests.

The nature of factors governing the tempo and mode of protein evolution is a fundamental issue in evolutionary biology. Specifically, whether or not interactions between different sites, or epistasis, are important in directing the course of evolution became one of the central questions. Several recent reports have scrutinized patterns of long-term protein evolution claiming them to be compatible only with an epistatic fitness landscape. However, these claims have not yet been substantiated with a formal model of protein evolution. Here, we formulate a simple covarion-like model of protein evolution focusing on the rate at which the fitness impact of amino acids at a site changes with time. We then apply the model to the data on convergent and divergent protein evolution to test whether or not the incorporation of epistatic interactions is necessary to explain the data. We find that convergent evolution cannot be explained without the incorporation of epistasis and the rate at which an amino acid state switches from being acceptable at a site to being deleterious is faster than the rate of amino acid substitution. Specifically, for proteins that have persisted in modern prokaryotic organisms since the last universal common ancestor for one amino acid substitution approximately ten amino acid states switch from being accessible to being deleterious, or vice versa. Thus, molecular evolution can only be perceived in the context of rapid turnover of which amino acids are available for evolution., (© The Author 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

The main forces directing long-term molecular evolution remain obscure. A sizable fraction of amino-acid substitutions seem to be fixed by positive selection, but it is unclear to what degree long-term protein evolution is constrained by epistasis, that is, instances when substitutions that are accepted in one genotype are deleterious in another. Here we obtain a quantitative estimate of the prevalence of epistasis in long-term protein evolution by relating data on amino-acid usage in 14 organelle proteins and 2 nuclear-encoded proteins to their rates of short-term evolution. We studied multiple alignments of at least 1,000 orthologues for each of these 16 proteins from species from a diverse phylogenetic background and found that an average site contained approximately eight different amino acids. Thus, without epistasis an average site should accept two-fifths of all possible amino acids, and the average rate of amino-acid substitutions should therefore be about three-fifths lower than the rate of neutral evolution. However, we found that the measured rate of amino-acid substitution in recent evolution is 20 times lower than the rate of neutral evolution and an order of magnitude lower than that expected in the absence of epistasis. These data indicate that epistasis is pervasive throughout protein evolution: about 90 per cent of all amino-acid substitutions have a neutral or beneficial impact only in the genetic backgrounds in which they occur, and must therefore be deleterious in a different background of other species. Our findings show that most amino-acid substitutions have different fitness effects in different species and that epistasis provides the primary conceptual framework to describe the tempo and mode of long-term protein evolution.

Background: The evolution and genomic stop codon frequencies have not been rigorously studied with the exception of coding of non-canonical amino acids. Here we study the rate of evolution and frequency distribution of stop codons in bacterial genomes., Results: We show that in bacteria stop codons evolve slower than synonymous sites, suggesting the action of weak negative selection. However, the frequency of stop codons relative to genomic nucleotide content indicated that this selection regime is not straightforward. The frequency of TAA and TGA stop codons is GC-content dependent, with TAA decreasing and TGA increasing with GC-content, while TAG frequency is independent of GC-content. Applying a formal, analytical model to these data we found that the relationship between stop codon frequencies and nucleotide content cannot be explained by mutational biases or selection on nucleotide content. However, with weak nucleotide content-dependent selection on TAG, -0.5 < Nes < 1.5, the model fits all of the data and recapitulates the relationship between TAG and nucleotide content. For biologically plausible rates of mutations we show that, in bacteria, TAG stop codon is universally associated with lower fitness, with TAA being the optimal for G-content < 16% while for G-content > 16% TGA has a higher fitness than TAG., Conclusions: Our data indicate that TAG codon is universally suboptimal in the bacterial lineage, such that TAA is likely to be the preferred stop codon for low GC content while the TGA is the preferred stop codon for high GC content. The optimization of stop codon usage may therefore be useful in genome engineering or gene expression optimization applications.

Abstract It is well known, that local B→A transformation in DNA is involved in several biological processes. In vitro B↔A transition is sequence-specific. The physical basis of this specificity is not known yet. Here we analyze the effect of intramolecular interactions on the structural behavior of the GG/CC and AA/TT steps. These steps exemplify sequence specific bias to the B- or A-form structure. Optimization of potential energy of the molecular systems composed of an octanucle-otide, neutralized by Na(+) and solvated with TIP3P water molecules in rectangular box with periodic boundary conditions gives the statistically representative sets of low energy structures for GG/CC and AA/TT steps in the middle of the diverse flanking sequences. Permissible 3D variations of GG/CC and AA/TT, and correlation of the relative motion of base pairs in these steps were analyzed. AA/TT step permits high variability for low energy conformers in the B-form DNA and small variability for low energy conformers in the A-form DNA. In contrast GG/CC step permits high variability for low energy conformers in the A-form DNA and small variability for low energy conformers in the B-form DNA. The relative motion of base pairs in GG/CC step is high correlated, while in AA/TT step this correlation is notably less. Atom-atom interactions inside-the-step always favors the B-form and their component - stacking interactions (atomatom interactions between nucleic bases) is crucial for the duplex stabilization. Formation of the A-form for both steps is a result of interactions with the flanking sequences and water-cation environment in the box. The average energy difference between conformations presenting B-form and A-form for the GG/CC step is high, while for the AA/TT step it is rather low. Thus, intramolecular interactions in GG/CC and AA/TT steps affect the possible conformational diversity ("conformational entropy") of the A- and B- type structures of DNA step. This determines the known bias of the A-form DNA depending on the enrichment of sequences with GG/CC. If structural tuning during the process of protein-DNA complex formation lead to the local B→A transformation of DNA, it is largely directed by high conformational diversity of GG/CC step in the A-form. In such a case the presence in the target site of both kinds of examined steps ensures the reversible character of ligand binding.

It is well known, that local B--> A transformation in DNA is involved in several biological processes. In vitro B<--> A transition is sequence-specific. The physical basis of this specificity is not known yet. Here we analyze the effect of intramolecular interactions on the structural behavior of the GG/CC and AA/TT steps. These steps exemplify sequence specific bias to the B- or A-form structure. Optimization of potential energy of the molecular systems composed of an octanucleotide, neutralized by Na(+) and solvated with TIP3P water molecules in rectangular box with periodic boundary conditions gives the statistically representative sets of low energy structures for GG/CC and AA/TT steps in the middle of the diverse flanking sequences. Permissible 3D variations of GG/CC and AA/TT, and correlation of the relative motion of base pairs in these steps were analyzed. AA/TT step permits high variability for low energy conformers in the B-form DNA and small variability for low energy conformers in the A-form DNA. In contrast GG/CC step permits high variability for low energy conformers in the A-form DNA and small variability for low energy conformers in the B-form DNA. The relative motion of base pairs in GG/CC step is high correlated, while in AA/TT step this correlation is notably less. Atom-atom interactions inside-the-step always favors the B-form and their component - stacking interactions (atom-atom interactions between nucleic bases) is crucial for the duplex stabilization. Formation of the A-form for both steps is a result of interactions with the flanking sequences and water-cation environment in the box. The average energy difference between conformations presenting B-form and A-form for the GG/CC step is high, while for the AA/TT step it is rather low. Thus, intramolecular interactions in GG/CC and AA/TT steps affect the possible conformational diversity ("conformational entropy") of the A- and B- type structures of DNA step. This determines the known bias of the A-form DNA depending on the enrichment of sequences with GG/CC. If structural tuning during the process of protein-DNA complex formation lead to the local B--> A transformation of DNA, it is largely directed by high conformational diversity of GG/CC step in the A-form. In such a case the presence in the target site of both kinds of examined steps ensures the reversible character of ligand binding.

An aim of the study was to assess clinical-economic effectiveness of therapy in the stages of examination of an epileptic patient by a local neurologist and then by an epileptologist after therapy optimization. Three hundreds and ten patients with different forms of epilepsy and types of epileptic seizures were treated. Costs-of-illness and medical services were assessed using common parameters on the basis of current prices in the health-care system. Special attention was drawn to expenses related to purchasing of antiepileptic drugs (traditional and new ones). The expenses for treatment were compared to treatment effectiveness. It is concluded that mean costs of treatment of epilepsy, including new brands of antiepileptic drugs with high prices, have substantial medical-social benefits that compensate expenses for their purchasing at the expense of other direct and indirect costs.

It has been found that 1500 tetrapeptides out of 160000 possible combinations occurring in proteins exhibit preference for particular conformational states. Most conformationally stable tetrapeptides obtained in the analysis of a sampling containing 706 proteins are in the alpha-helical form. The features of the amino acid composition of conformationally stable oligopeptides have been studied.

Regions with a quasiperiodical location of exon--intron sites have been found in the loci of genes of I and VII type collagens (with a total length of exons more than 15% of the entire size of the locus). The periods observed are similar to periods typical for the nucleosomal level of the organization of chromatin. It was shown that the sites consisting of successively arranged exons and introns form groups involving two to five such regions of the same length. The groups encoding the fibrillar regions of the gene product contain more than 50% of exons. The regions are on the average 165 nt long, which is close to the minimal nucleosomal repeat length observed in some regions of the eukaryotic genome. In the nonfibrillar region of the gene of VII type collagen, groups of several exon-intron pairs with an average length of 227 nt were identified. The change in the length of exon-introns sites on going from the nonfibrillar to the fibrillar moiety occurs in a jump, which is clearly seen on a periodogram of the locus.

Recognition of coding regions within eukaryotic genomes is one of oldest but yet not solved problems of bioinformatics. New high-accuracy methods of splicing sites recognition are needed to solve this problem. A question of current interest is to identify specific features of nucleotide sequences nearby splicing sites and recognize sites in sequence context. We performed a statistical analysis of human genes fragment database and revealed some characteristics of nucleotide sequences in splicing sites neighborhood. Frequencies of all nucleotides and dinucleotides in splicing sites environment were computed and nucleotides and dinucleotides with extremely high\low occurrences were identified. Statistical information obtained in this work can be used in further development of the methods of splicing sites annotation and exon-intron structure recognition.

We identified a mutation in the CRYGD gene (P23S) of the gamma-crystallin gene cluster that is associated with a polymorphic congenital cataract that occurs with frequency of approximately 0.3% in a human population. To gain insight into the molecular mechanism of the pathogenesis of gamma-crystallin isoforms, we undertook an evolutionary analysis of the available mammalian and newly obtained primate sequences of the gamma-crystallin genes. The cataract-associated serine at site 23 corresponds to the ancestral state, since it was found in CRYGD of a lower primate and all the surveyed nonprimate mammals. Crystallin proteins include two structurally similar domains, and substitutions in mammalian CRYGD protein at site 23 of the first domain were always associated with substitutions in the structurally reciprocal sites 109 and 136 of the second domain. These data suggest that the cataractogenic effect of serine at site 23 in the N-terminal domain of CRYGD may be compensated indirectly by amino acid changes in a distal domain. We also found that gene conversion was a factor in the evolution of the gamma-crystallin gene cluster throughout different mammalian clades. The high rate of gene conversion observed between the functional CRYGD gene and two primate gamma-crystallin pseudogenes (CRYGEP1 and CRYGFP1) coupled with a surprising finding of apparent negative selection in primate pseudogenes suggest a deleterious impact of recently derived pseudogenes involved in gene conversion in the gamma-crystallin gene cluster.

Having an extra copy of a gene is thought to provide some functional redundancy, which results in a higher rate of evolution in duplicated genes. In this article, we estimate the impact of gene duplication on the selection of tuf paralogs, and we find that in the absence of gene conversion, tuf paralogs have evolved significantly slower than when gene conversion has been a factor in their evolution. Thus, tuf gene copies evolve under a selective pressure that ensures their functional uniformity, and gene conversion reduces selection against amino acid substitutions that affect the function of the encoded protein, EF-Tu.

We describe a new method for polyproline II-type (PPII) secondary structure prediction based on tetrapeptide conformation properties using data obtained from all globular proteins in the Protein Data Bank (PDB). This is the first method for PPII prediction with a relatively high level of accuracy (approximately 60%). Our method uses only frequencies of different conformations among oligopeptides without any additional parameters. We also attempted to predict alpha-helices and beta-strands using the same approach. We find that the application of our method reveals interrelation between sequence and structure even for very short oligopeptides (tetrapeptides)., (Proteins 2005. 2005 Wiley-Liss, Inc.)

Alignment of protein sequences is a key step in most computational methods for prediction of protein function and homology-based modeling of three-dimensional (3D)-structure. We investigated correspondence between "gold standard" alignments of 3D protein structures and the sequence alignments produced by the Smith-Waterman algorithm, currently the most sensitive method for pair-wise alignment of sequences. The results of this analysis enabled development of a novel method to align a pair of protein sequences. The comparison of the Smith-Waterman and structure alignments focused on their inner structure and especially on the continuous ungapped alignment segments, "islands" between gaps. Approximately one third of the islands in the gold standard alignments have negative or low positive score, and their recognition is below the sensitivity limit of the Smith-Waterman algorithm. From the alignment accuracy perspective, the time spent by the algorithm while working in these unalignable regions is unnecessary. We considered features of the standard similarity scoring function responsible for this phenomenon and suggested an alternative hierarchical algorithm, which explicitly addresses high scoring regions. This algorithm is considerably faster than the Smith-Waterman algorithm, whereas resulting alignments are in average of the same quality with respect to the gold standard. This finding shows that the decrease of alignment accuracy is not necessarily a price for the computational efficiency., (Copyright 2003 Wiley-Liss, Inc.)

In order to estimate the influence of point mutations on the protein structure, 552 pairs of structurally aligned proteins with pairwise sequence identities greater than 60% were analyzed. We selected all isolated point mismatches from these alignments. The statistics of local conformational changes corresponding to these mismatches was studied. This study revealed two surprising aspects: (a) only an extremely minor fraction (< 3%) of all observed mutations leads to significant changes in local conformations; (b) observed substitutions, which affect local conformation are more frequently those of similar (high scoring) amino acid pairs rather than the substitutions thought of as "hazardous" (low scoring), although the overall average score of the changing substitutions is still lower than the corresponding average score for the substitutions that leave the structure unchanged.