Modelling the Efficiency of Codon–tRNA Interactions Based on Codon Usage Bias

The tRNA adaptation index (tAI) is a widely used measure of the efficiency by which a coding sequence is recognized by the intra-cellular tRNA pool. This index includes among others weights that represent wobble interactions between codons and tRNA molecules. Currently, these weights are based only on the gene expression in Saccharomyces cerevisiae. However, the efficiencies of the different codon–tRNA interactions are expected to vary among different organisms. In this study, we suggest a new approach for adjusting the tAI weights to any target model organism without the need for gene expression measurements. Our method is based on optimizing the correlation between the tAI and a measure of codon usage bias. Here, we show that in non-fungal the new tAI weights predict protein abundance significantly better than the traditional tAI weights. The unique tRNA–codon adaptation weights computed for 100 different organisms exhibit a significant correlation with evolutionary distance. The reported results demonstrate the usefulness of the new measure in future genomic studies.     

Arrangement of the Template 3′ of the A-Site Codon on the Human 80S Ribosome

The arrangement of the template sequence 3′ of the A-site codon on the 80S ribosome was studied using mRNA analogs containing Phe codon UUU at the 5′ end and a photoreactive perfluoroarylazido group linked to C5 of U or N7 of G. The analogs were positioned on the ribosome with the use of tRNA^Phe, which directed the UUU codon to the P site, bringing a modified nucleotide to position +9 or +12 relative to the first nucleotide of the P-site codon. Upon mild UV irradiation of ribosome complexes, the analogs of both types crosslinked to the 18S rRNA and proteins of the 40S subunit. Comparisons were made with the crosslinking patterns of complexes in which an mRNA analog contained a modified nucleotide in position +7 (the crosslinking to 18S rRNA in such complexes has been studied previously). The efficiency of crosslinking to ribosomal components depended on the nature of the modified nucleotide of an mRNA analog and its position on the ribosome. The extent of crosslinking to the 18S rRNA drastically decreased as the modified nucleotide was transferred from position +7 to position +12. The 18S rRNA nucleotides involved in crosslinking were identified. A modified nucleotide in position +9 crosslinked to the invariant dinucleotide A1824/A1825 and variable A1823 in the 3′ minidomain of the 18S rRNA and to S15. The same ribosomal components have earlier been shown to crosslink to modified nucleotides in positions +4 to +7. In addition, all mRNA analogs crosslinked to invariant C1698 in the 3′ minidomain and to conserved region 605–620, which closes helix 18 in the 5′ domain.     

Template Location on the Human Ribosome: Environment of the mRNA Nucleotide Adjacent to the A-Site Codon on the 3′-Side

The 18S rRNA nucleotides close to the template nucleotide adjacent to the 80S ribosomal A-site codon on the 3′-end (i.e., the nucleotide in position +7 relative to the first nucleotide of the P-site codon) were identified using the affinity crosslinking approach. For this purpose, the photoreactive mRNA analogues with a perfluorophenylazide group attached through various linkers to the uridine C5, 3′-terminal phosphate or guanosine N7 were used. The position of the mRNA analogues on the ribosome was preset using tRNA^Phe, which recognized the phenylalanine codon directed to the P-site. An analysis of the rRNAs isolated from the irradiated complexes of 80S ribosomes showed that all the analogues are almost equally crosslinked to the 18S rRNA nucleotides we attributed to the A-site codon environment: namely, to nucleotides A1823, A1824, and A1825 of the 3′-minidomain and to the 620–630 fragment of the 18S rRNA 5′-domain. In addition, we identified a new component of the mRNA binding site of human ribosomes, nucleotide C1698, belonging to the 18S rRNA 3′-minidomain, using analogues bearing a perfluorophenylazide group on uridine and guanine residues.__________Translated from Bioorganicheskaya Khimiya, Vol. 31, No. 3, 2005, pp. 295–302.Original Russian Text Copyright © 2005 by Demeshkina, Styazhkina, Bulygin, Repkova, Ven’yaminova, Karpova.     

Codon–anticodon interaction and the genetic code evolution

The evolution of the genetic code, with 20 amino acids encoded from the beginning, is analyzed from the viewpoint of codon–anticodon interaction. Imposing a minimum principle for the interaction, in the framework of the so called crystal basis model of the genetic code, we determine the structure of the anticodons in the ancient, archetypal and early genetic codes, that are all reconciled in a unique frame. Most of our results agree with the generally accepted scheme.     

Highly Conserved Region 1816–1831 of the 18S Ribosomal RNA Is Close to the First Nucleotide of the A-Site Codon in Elongation and Termination of Translation

Two mRNA analogs, pUUCUAAA (with stop codon UAA) and pUUCUCAA (with Ser codon UCA) containing a perfluoroarylazido group at U4, were used to study the position relative to the 18S rRNA for the first nucleotide of the codon located in the A site of the human 80S ribosome. To place UAA or UCA in the A site, UCC-recognizing tRNA^Phe was bound in the P site. With each analog, crosslinking was detected for highly conserved fragment 1816–1831, which contains invariant dinucleotide A1823/A1824 and is in helix 44 at the 3" end of the 18S rRNA. Since 18S rRNA modification did not depend on whether the U4 photoreactive group was in the sense or stop codon, it was assumed that polypeptide chain release factor 1 directly recognizes the trinucleotide of a stop codon located in the A site.     

The Template Region 5′ of the E-Site Codon Is Close to S26 in the Human 80S Ribosome

The environment of the template sequence 5 of the E-site codon on the 80S ribosome was studied with nonaribonucleotide or dodecaribonucleotide derivatives containing Phe codon UUU at the 3 end and a perfluoroarylazido group at the first or third nucleotide. A photoreactive group was linked to C5 of U or N7 of G. The analogs were positioned on the ribosome with the use of tRNA^Phe, which is cognate to the UUU codon and directs it to the P site, bringing a modified nucleotide in position –4 to –9 relative to the first nucleotide of the P-site codon. Upon irradiation of ribosome complexes with tRNA^Phe and the mRNA analogs with mild UV light, the analogs crosslinked predominantly to the 40S subunit, modifying the proteins. The major target of modification was S26 in all cases. In addition, S3 was modified to a low extent when the reactive nucleotide was in position –4 and S14 was in position –6. In the absence of tRNA, all mRNA analogs modified S3.     

Evidence of Gene Conversion in the Evolutionary Process of the Codon 41/42 (-CTTT) Mutation Causing β-Thalassemia in Southern China

The 4-bp deletion (-CTTT) at codon 41/42 (CD41/42) of the human beta-globin gene represents one of the most common beta-thalassemia mutations in East Asia and Southeast Asia, which is historically afflicted with endemic malaria, thus hypothetically evolving under natural selection by malaria infection. To understand the evolutionary process of generating the beta(CD41/42) allele and its maintenance, including the effect of natural selection on the pattern of linkage disequilibrium (LD), we sequenced a 15.933-kb region spanning 20.693 kb of the beta-globin cluster surrounding the 4-bp deletion using a sample from a Chinese population consisting of 24 normal individuals and 16 heterozygotes for the deletion. Forty-nine polymorphic sites were found. Analysis of the data, using a variety of methods including formal population genetics analysis and visual approaches, suggests that the spread of the CD41/42 (-CTTT) deletion is most likely mediated by interallelic gene conversion, although independent deletions in different lineages are also possible. The neutrality test resulted in a significant positive Tajima's D for the beta-globin locus, which is consistent with the existence of balancing selection. This suggests that the 4-bp deletion that occurred at this locus may be an event that is subject to natural selection, due to malaria, which leads to the heterozygote advantage, spread widely with further help by conversion and migration. The evolutionary process of this mutant through gene conversion that could conceivably take place between the 4-bp deletion and the normal sequence in the respective region is discussed in detail.     

Synthesis and Properties of an Initiation Codon Analog Consisting of 2′-O-Methyl Nucleotides

Abstract

An initiation codon analog consisting of 2′-O-methyl nucleotides (AmUmG) was synthesized and examined for its binding efficiency to E. coli ribosomes with fMetRNA^fMet and also its stability in binding assay systems for comparison with those of r(AUG) and d(ATG). AmUmG was found completely resistant to nucleases under the conditions used.     

Correlation Between Shine–Dalgarno Sequence Conservation and Codon Usage of Bacterial Genes

In this study, we analyzed the correlation between codon usage bias and Shine–Dalgarno (SD) sequence conservation, using complete genome sequences of nine prokaryotes. For codon usage bias, we adopted the codon adaptation index (CAI), which is based on the codon usage preference of genes encoding ribosomal proteins, elongation factors, heat shock proteins, outer membrane proteins, and RNA polymerase subunit proteins. To compute SD sequence conservation, we used SD motif sequences predicted by Tompa and systematically aligned them with 5′UTR sequences. We found that there exists a clear correlation between the CAI values and SD sequence conservation in the genomes of Escherichia coli, Bacillus subtilis, Haemophilus influenzae, Archaeoglobus fulgidus, Methanobacterium thermoautotrophicum, and Methanococcus jannaschii, and no relationship is found in M. genitalium, M. pneumoniae, and Synechocystis. That is, genes with higher CAI values tend to have more conserved SD sequences than do genes with lower CAI values in these organisms. Some organisms, such as M. thermoautotrophicum, do not clearly show the correlation. The biological significance of these results is discussed in the context of the translation initiation process and translation efficiency. Received: 22 June 2000 / Accepted: 18 October 2000     

New Insights into the Fructosyltransferase Activity of Schwanniomyces occidentalis β-Fructofuranosidase,Emerging from Nonconventional Codon Usage and Directed Mutation

Schwanniomyces occidentalis β-fructofuranosidase (Ffase) releases β-fructose from the nonreducing ends of β-fructans and synthesizes 6-kestose and 1-kestose, both considered prebiotic fructooligosaccharides. Analyzing the amino acid sequence of this protein revealed that it includes a serine instead of a leucine at position 196, caused by a nonuniversal decoding of the unique mRNA leucine codon CUG. Substitution of leucine for Ser196 dramatically lowers the apparent catalytic efficiency (k_cat/K_m) of the enzyme (approximately 1,000-fold), but surprisingly, its transferase activity is enhanced by almost 3-fold, as is the enzymes'' specificity for 6-kestose synthesis. The influence of 6 Ffase residues on enzyme activity was analyzed on both the Leu196/Ser196 backgrounds (Trp47, Asn49, Asn52, Ser111, Lys181, and Pro232). Only N52S and P232V mutations improved the transferase activity of the wild-type enzyme (about 1.6-fold). Modeling the transfructosylation products into the active site, in combination with an analysis of the kinetics and transfructosylation reactions, defined a new region responsible for the transferase specificity of the enzyme.β-Fructofuranosidases (EC 3.2.1.26) are enzymes of biotechnological interest that catalyze the release of β-fructose from the nonreducing termini of various β-d-fructofuranoside substrates. In general, they exhibit a high degree of sequence homology, and based on their amino acid sequences, they fall into family 32 of the glycosyl-hydrolases (GH), along with invertases, inulinases, and fructosyltransferases (http://www.cazy.org). The GH32 family has been studied intensely, and some three-dimensional structures are now available, such as that of inulinase from Aspergillus awamorii (26), fructan-exohydrolase from Cichorium intybus (CiFEH) (34, 36), or invertase from Thermotoga maritima (2, 3) and Arabidopsis thaliana (35). These proteins contain a five-blade β-propeller N-terminal catalytic module and a C-terminal β-sandwich domain (19). Multiple-sequence alignment of GH32 proteins, which are included in the GH-J clan together with the GH68 proteins of the inulosucrase family, reveals the presence of three conserved motifs, each containing a key acidic residue (in boldface) implicated in substrate binding and hydrolysis: Asn-Asp-Pro-Asn-Gly (NDPNG), Arg-Asp-Pro (RDP), and Glu-Cys (EC) (28). These conserved residues are implicated in a double-displacement reaction in which a covalent glycosyl-enzyme intermediate is formed. Thus, the catalytic mechanism proposed for the Saccharomyces cerevisiae invertase implies that Asp23 (NDPNG) acts as a nucleophile and Glu204 (EC) acts as the acid/base catalyst (29), whereas Asp309 (RDP) of Acetobacter diazotropicus levansucrase influences the efficiency of sucrose hydrolysis (7) and Arg188 and Asp189 of the latter motif define the substrate binding and specificity of exoinulinase from A. awamorii toward fructopyranosyl residues (26).As well as hydrolyzing sucrose, β-fructofuranosidases may also catalyze the synthesis of short-chain fructooligosaccharides (FOS), in which one to three fructosyl moieties are linked to the sucrose skeleton by different glycosidic bonds, depending on the source of the enzyme (12, 21, 31). FOS act as prebiotics, and they exert a beneficial effect on human health, participating in the prevention of cardiovascular diseases, colon cancer, and osteoporosis (16). Currently, FOS are mainly produced by Aspergillus fructosyltransferase in industry (10, 31), providing a mixture of FOS with an inulin-type structure that contains β-(2→1)-linked fructose oligomers (¹F-FOS: 1-kestose or nystose). Curiously, when the link between two fructose units (⁶F-FOS: 6-kestose) or between fructose and the glucosyl moiety (⁶G-FOS: neokestose) involves a β-(2→6) link, the prebiotic properties of the FOS may be enhanced beyond that of commercial FOS (23).The yeast Schwanniomyces occidentalis (also called Debaryomyces occidentalis) produces a number of extracellular enzymes that make it of interest in biotechnology. Several of its amylolytic enzymes have been characterized, including amylases and glucoamylase (1, 9), as well as an invertase (17). In addition, we also characterized an extracellular β-fructofuranosidase (Ffase) from this yeast that hydrolyzes sucrose, 1-kestose, and nystose (5). This enzyme exhibited a transfructosylating activity that efficiently produces the trisaccharides 6-kestose and 1-kestose in the ratio 3:1, generating the highest 6-kestose yield yet reported, as far as we know. The Ffase three-dimensional structure has recently been solved (6) and represented as a homodimer, each modular subunit arranged like other GH32 enzymes. The Asp50 (NDPNG) and Glu230 (EC) located at the center of the propeller are the catalytic residues implicated in substrate binding and hydrolysis, whereas Arg178 and Asp179 form the RDP motif (6).The genetic codes of some yeasts incorporate certain variations. For example, while CUG was believed to be a universal codon for leucine, in the cytoplasm of certain species of the genus Candida (15) it encodes a serine, as in Pichia farinosa (33). The reassignment of this codon is mediated by a novel serine-tRNA that acquired a leucine 5′-CAG-3′ anticodon (25).Here, we show that deviation from the standard use of the CUG leucine codon to encode serine was correlated with the transferase capacity and specificity of the Ffase enzyme. Indeed, the S196L substitution enhanced the transferase activity of the enzyme 3-fold. Several site-directed mutants were generated and characterized to study their transferase capacities. These results are considered on the basis of the enzymes'' three-dimensional structure, which enables a novel putative binding site of sucrose that serves as a water substitute donor in the hydrolytic reaction yielding the tranglycosylation product 6-kestose to be identified.     

The Highly Conserved Codon following the Slippery Sequence Supports ?1 Frameshift Efficiency at the HIV-1 Frameshift Site

HIV-1 utilises −1 programmed ribosomal frameshifting to translate structural and enzymatic domains in a defined proportion required for replication. A slippery sequence, U UUU UUA, and a stem-loop are well-defined RNA features modulating −1 frameshifting in HIV-1. The GGG glycine codon immediately following the slippery sequence (the ‘intercodon’) contributes structurally to the start of the stem-loop but has no defined role in current models of the frameshift mechanism, as slippage is inferred to occur before the intercodon has reached the ribosomal decoding site. This GGG codon is highly conserved in natural isolates of HIV. When the natural intercodon was replaced with a stop codon two different decoding molecules—eRF1 protein or a cognate suppressor tRNA—were able to access and decode the intercodon prior to −1 frameshifting. This implies significant slippage occurs when the intercodon is in the (perhaps distorted) ribosomal A site. We accommodate the influence of the intercodon in a model of frame maintenance versus frameshifting in HIV-1.     

Optimization of the AT-content of Codons Immediately Downstream of the Initiation Codon and Evaluation of Culture Conditions for High-level Expression of Recombinant Human G-CSF in <Emphasis Type="Italic">Escherichia coli</Emphasis>

Enhanced therapeutic importance of recombinant human granulocyte colony stimulating factor (rhG-CSF) has encouraged us to develop a processing method for its high-level expression in E. coli. In this study, we established a high-yielding clone by incorporation of silent mutations at N-terminal region of human G-CSF gene. We studied and optimized various parameters of culture conditions connected with the expression of rhG-CSF. The maximum expression was obtained in a defined medium supplemented with 1% glucose. The gene in pET-3a vector in E. coli BL21 (DE3) PLysS host strain was induced with 2 mM isopropyl β-d-1-thiogalacto pyronoside. The cell growth and productivity was enhanced about 1.6- and 1.5-folds, respectively when inducing the culture at OD₆₀₀ value of 6 than 2. The protein expression was significantly increased by addition of rifampicin at concentration of 200 μg/ml. The AT content of 51.8% with suitable codon sequences at N-terminal region and the concentration of rifampicin were identified as the key factors with a significant impact on protein expression. The specific productivity of 104 mg/OD/l (68.7% of total cellular protein) of rhG-CSF was obtained toward the end of the study, which is almost 1.5 times higher yield than reported so far in the literature.     

A Stable Upstream Stem-loop Structure Enhances Selection of the First 5′-ORF-AUG as a Main Start Codon for Translation Initiation of Human ACAT1 mRNA

Human ACAT1 cDNA K1 was first cloned and functionally expressed in 1993. There are two adjacent in-frame AUG codons, AUG1397-1399 and AUG1415-1417, at 5′-terminus of the open reading frame (ORF,nt 1397-3049) of human ACAT1 mRNA corresponding to cDNA K1. In current work, these two adjacent in-frame AUGs at 5′-terminus of the predicted ORF (5′-ORF-AUGs) as start codons for translation initiation of human ACAT1 mRNA were characterized in detail. Codon mutations indicated that both of these two adjacent 5′-ORF-AUGs can be selected as start codons but the first 5′-ORF-AUG1397-1399 is a main start codon consistent with that of the predicted ORF of human ACAT1 mRNA. Further deletion and mutation analyses demonstrated that a stable upstream stem-loop structure enhanced the selection of the first 5′-ORF-AUG1397-1399 as a main start codon, in addition to upstream nucleotide A in the -3 position, which is a key site of Kozak sequence. In addition, result of ACAT1 enzymatic activity assay showed no obvious difference between these two ACAT1 proteins respectively initiated from the two adjacent 5′-ORF-AUGs. This work showed that astable upstream stem-loop structure could modulate the start codon selection during translation initiation of mRNAs that contain adjacent multi-5′-ORF-AUGs.     

A Generalized Michaelis–Menten Equation in Protein Synthesis: Effects of Mis-Charged Cognate tRNA and Mis-Reading of Codon

The sequence of amino acid monomers in the primary structure of a protein is decided by the corresponding sequence of codons (triplets of nucleic acid monomers) on the template messenger RNA (mRNA). The polymerization of a protein, by incorporation of the successive amino acid monomers, is carried out by a molecular machine called ribosome. We develop a stochastic kinetic model that captures the possibilities of mis-reading of mRNA codon and prior mis-charging of a tRNA. By a combination of analytical and numerical methods, we obtain the distribution of the times taken for incorporation of the successive amino acids in the growing protein in this mathematical model. The corresponding exact analytical expression for the average rate of elongation of a nascent protein is a ‘biologically motivated’ generalization of the Michaelis–Menten formula for the average rate of enzymatic reactions. This generalized Michaelis–Menten-like formula (and the exact analytical expressions for a few other quantities) that we report here display the interplay of four different branched pathways corresponding to selection of four different types of tRNA.     

Effect of cortisone on the developmental pattern of the neutral and the acid β-galactosidases of the small intestine of the rat

1. The developmental pattern and effect of cortisone on acid beta-galactosidase and neutral beta-galactosidase were studied in postnatal rats by a recently proposed method for their independent determination. 2. After birth the acid beta-galactosidase activity increases in the ileum, whereas it decreases slightly in the jejunum. On day 16 after birth the activity in the ileum decreases and in 20-day-old rats activity in both parts of the intestine decreases to adult values. In suckling animals the activity in the ileum exceeds the jejunal activity severalfold and in adult animals the activity in the jejunum is slightly higher than that in the ileum. 3. Neutral beta-galactosidase activity is high after birth and decreases in both jejunum and ileum after day 20 after birth. In 12-20-day-old rats activity in both parts is essentially the same, but in adult animals jejunal activity exceeds ileal activity four-to five-fold. 4. Cortisone (0.5, 2.0 or 5.0mg/100g body wt. daily for 4 days) does not influence the activity of either enzyme in 60-day-old rats. Acid beta-galactosidase activity is decreased after cortisone treatment in 8-, 12-, 16-and 18-day-old rats, with sensitivity to cortisone increasing with the approach of weaning. No effect of cortisone on acid beta-galactosidase is seen in 8-day-old rats. Neutral beta-galactosidase activity is increased in the ileum of 8-, 12-, 16- and 18-day old rats, but only in the jejunum of 8-and 12-day-old rats.     

Association of the subunits of the calcium-independent receptor of α-latrotoxin

CIRL-1 also called latrophilin 1 or CL belongs to the family of adhesion G protein-coupled receptors (GPCRs). As all members of adhesion GPSR family CIRL-1 consists of two heterologous subunits, extracellular hydrophilic p120 and heptahelical membrane protein p85. Both CIRL-1 subunits are encoded by one gene but as a result of intracellular proteolysis of precursor, mature receptor has two-subunit structure. It was also shown that a minor portion of the CIRL-1 receptor complexes dissociates, producing the soluble receptor ectodomain, and this dissociation is due to the second cleavage at the site between the site of primary proteolysis and the first transmembrane domain. Recently model of independent localization p120 and p85 on the cell surface was proposed. In this article we evaluated the amount of p120-p85 complex still presented on the cellular membrane and confirmed that on cell surface major amount of mature CIRL-1 presented as a p120-p85 subunit complex.     

Phylogeography of the harlequin beetle-riding pseudoscorpion and the rise of the Isthmus of Panamá

Molecular and geological evidence indicates that the emergence of the Isthmus of Panamá influenced the historical biogeography of the Neotropics in a complex, staggered manner dating back at least 9 Myr bp. To assess the influence of Isthmus formation on the biogeography of the harlequin beetle-riding pseudoscorpion, Cordylochernes scorpioides, we analysed mitochondrial COI sequence data from 71 individuals from 13 locations in Panamá and northern South America. Parsimony and likelihood-based phylogenies identified deep divergence between South American and Panamanian clades. In contrast to low haplotype diversity in South America, the Panamanian Cordylochernes clade is comprised of three highly divergent lineages: one clade consisting predominantly of individuals from central Panamá (PAN A), and two sister clades (PAN B1 and PAN B2) of western Panamanian pseudoscorpions. Breeding experiments demonstrated a strictly maternal mode of inheritance, indicating that our analyses were not confounded by nuclear-mitochondrial pseudogenes. Haplotype diversity is striking in western Atlantic Panamá, where all three Panamanian clades can occur in a single host tree. This sympatry points to the existence of a cryptic species hybrid zone in western Panamá, a conclusion supported by interclade crosses and coalescence-based migration rates. Molecular clock estimates yield a divergence time of approximately 3 Myr between the central and western Panamanian clades. Taken together, these results are consistent with a recent model in which a transitory proto-Isthmus enabled an early wave of colonization out of South America at the close of the Miocene, followed by sea level rise, inundation of the terrestrial corridor and then a second wave of colonization that occurred when the Isthmus was completed approximately 3 Myr bp.     

Synonymous Codon Usage Bias and Overexpression of a Synthetic Gene Encoding Interferon α2b in Yeast

To achieve higher level expression of Interferon α2b(IFN-α2b)in methylotrophic yeast(Pichia pastoris),a cDNA fragment coding for the mature IFN-α2b was designed and synthesized based on the synonymous codon bias of P.pastoris and optimized G C content.The synthetic IFN-α2b was inserted into the secreted expression vector pPICZαA,and then integrated into P.pastoris GS115 genome by electroporation.Multi-copy integrants in the Mut recombinant P.pastoris strain were screened by high concentrations of Zeocin.120 hours culturing allowed expression of the IFN-α2b transformant up to 810 mg/L as detected by SDS-PAGE and quantitative methods.In addition,Western blot analysis showed that the recombinant proteins had immunogenicity.The significant antiviral activity of the recombinant IFN-α2b protein was verified by WISH/VSV system,which was 3.3×105 IU/mL.