Preferential codon usage in prokaryotic genes: the optimal codon-anticodon interaction energy and the selective codon usage in efficiently expressed genes

By considering the nucleotide sequence of several highly expressed coding regions in bacteriophage MS2 and mRNAs from Escherichia coli, it is possible to deduce some rules which govern the selection of the most appropriate synonymous codons NNU or NNC read by tRNAs having GNN, QNN or INN as anticodon. The rules fit with the general hypothesis that an efficient in-phase translation is facilitated by proper choice of degenerate codewords promoting a codon-anticodon interaction with intermediate strength (optimal energy) over those with very strong or very weak interaction energy. Moreover, codons corresponding to minor tRNAs are clearly avoided in these efficiently expressed genes. These correlations are clearcut in the normal reading frame but not in the corresponding frameshift sequences +1 and +2. We hypothesize that both the optimization of codon-anticodon interaction energy and the adaptation of the population to codon frequency or vice versa in highly expressed mRNAs of E. coli are part of a strategy that optimizes the efficiency of translation. Conversely, codon usage in weakly expressed genes such as repressor genes follows exactly the opposite rules. It may be concluded that, in addition to the need for coding an amino acid sequence, the energetic consideration for codon-anticodon pairing, as well as the adaptation of codons to the tRNA population, may have been important evolutionary constraints on the selection of the optimal nucleotide sequence.     

Organization of plastid-encoded ATPase genes and flanking regions including homologues of infB and tsf in the thermophilic red alga Galdieria sulphuraria

We have cloned and sequenced the plastid ATPase operons (atp1 and atp2) and flanking regions from the unicellular red alga Galdieria sulphuraria (Cyanidium caldarium). Six genes (5 atpI, H, G, F, D and A 3) are linked in atp1 encoding ATPase subunits a, c, b, b, and , respectively. The atpF gene does not contain an intron and overlaps atpD by 1 bp. As in the genome of chloroplasts from land plants, the cluster is located downstream of rps2, but between this gene and atp1 we found the gene for the prokaryotic translation elongation factor TS. Downstream of atpA, we detected two open reading frames, one encoding a putative transport protein. The genes atpB and atpE, encoding ATPase subunits and , respectively, are linked in atp2, seperated by a 2 bp spacer. Upstream of atpB, an uninterrupted orf167 was detected which is homologous to an intron-containing open reading frame in land plant chloroplasts. This orf167 is preceded on the opposite DNA strand by a homologue to initiation factor 2 in prokaryotes. The arrangement of atp1 and atp2 is the same as observed in the multicellular red alga Antithamnion sp. indicatiing a conserved genome arrangement in the red algal plastid genome. Differences compared to green chloroplast genomes suggest a large phylogenetic distance between red algae and green plants, while similarities in arrangement and sequence to chromophytic ATPase operons support a red algal origin of chlorophyll a/c-containing plastids or alternatively point to a common prokaryotic endosymbiont.     

Nucleotide sequence of the gene coding for the elongation factor Tu from the extremely thermophilic eubacterium Thermotoga maritima

Abstract The gene coding for the elongation factor Tu (EF-Tu) of Thermatoga maritima was cloned and sequenced. The predicted amino acid sequence was compared with those of other eubacteria, an archaebacterium and two eukaryotes as well. The similarity values and the distance matrix tree show that Thermotoga is more closely related to the eubacteria than to the representatives of the other urkingdoms. Thermotoga maritima represents the deepest branching within the tree of EF-Tu sequences from all eubacteria studied so far.     

Everything in moderation: Archaea as ‘non-extremophiles’

Well characterized and cultivated archaea are prokaryotic specialists that thrive in habitats of elevated temperature, low pH, high salinity, or strict anoxia. Recently, however, new groups of abundant, uncultivated archaea have been found to be widespread in more pedestrian biotopes, including marine plankton, terrestrial soils, lakes, marine and freshwater sediments, and in association with metazoa. Research efforts are presently focused on characterizing the physiology, biochemistry and genetics of these abundant and cosmopolitan but poorly understood archaea.     

Crystal structure of Ca2+ -free S100A2 at 1.6-A resolution

S100A2 is an EF hand-containing Ca^2 +-binding protein of the family of S100 proteins. The protein is localized exclusively in the nucleus and is involved in cell cycle regulation. It attracted most interest by its function as a tumor suppressor via p53 interaction. We determined the crystal structure of homodimeric S100A2 in the Ca^2 +-free state at 1.6-Å resolution. The structure revealed structural differences between subunits A and B, especially in the conformation of a loop that connects the N- and C-terminal EF hands and represents a part of the target-binding site in S100 proteins. Analysis of the hydrogen bonding network and molecular dynamics calculations indicate that one of the two observed conformations is more stable. The structure revealed Na⁺ bound to each N-terminal EF hand of both subunits coordinated by oxygen atoms of the backbone carbonyl and water molecules. Comparison with the structures of Ca^2 +-free S100A3 and S100A6 suggests that Na⁺ might occupy the S100-specific EF hand in the Ca^2 +-free state.     

Exploring structurally conserved solvent sites in protein families

Protein-bound water molecules are important components of protein structure, and therefore, protein function and energetics. Although structural conservation of solvent has been studied in a few protein families, a lack of suitable computational tools has hindered more comprehensive analyses. Herein we present a semiautomated computational approach for identifying solvent sites that are conserved among proteins sharing a common three-dimensional structure. This method is tested on six protein families: (1) monodomain cytochrome c, (2) fatty-acid binding protein, (3) lactate/malate dehydrogenase, (4) parvalbumin, (5) phospholipase A2, and (6) serine protease. For each family, the method successfully identified previously known conserved solvent sites. Moreover, the method discovered 22 novel conserved solvent sites, some of which have higher degrees of conservation than the previously known sites. All six families studied had solvent sites with more than 90% conservation and these sites were invariably located in regions of the protein with very high sequence conservation. These results suggest that highly conserved solvent sites, by virtue of their proximity to conserved residues, should be considered as one of the defining three-dimensional structural characteristics of protein families and folds.     

Centrin isoforms in mammals. Relation to calmodulin

In mammals, three calmodulin (CaM) genes code for 100% identical proteins. In these species, four centrin (Cetn) genes have been reported to exist. They are examined in this paper. While the gene for Cetn 1 contains no introns and appears to be derived from Cetn 2 by retroposition, a gene product for Cetn 1 is expressed. Cetn 2, 3, and 4 represent bona fide genes. The major difference between the members of the CaM and the Cetn subfamilies is the presence (usually) in Cetn of an approximately 23 amino acids long (but occasionally much longer) protruding amino acid end. In all members of these two subgroups, four EF hand motifs (in this paper taken as loops containing 12 amino acids) are separated by 24, 25 and 24 amino acids (each a helix–loop–helix) positioned between motifs 1and 2, 2 and 3, and 3 and 4, respectively. This rule applies not only to CaM and Cetn in mammals but also to these two subfamilies in simpler eukaryotes such as Saccharomyces cerevisiae and Giardia lamblia. The various mRNA products can be identified most readily by their characteristic 3′ UTRs. While CaM is an ancient molecule that is expressed in all cells and is ubiquitous within these cells and interacts therein with almost 100 different proteins, many of which display the IQ or related binding motifs, the distribution and function of Cetn (an equally ancient molecule) is restricted mostly to basal bodies (e.g. in rods of the retina), axonemes, flagella, cilia and centrosomes. Are these two subclasses of calcium carriers (each molecule possessing four EF hands which possibly interact with different association constants)—if they are both present within a cell—randomly chosen for their service to the specific proteins with which they interact?     

Exit Strategies for Charged tRNA from GluRS

For several class I aminoacyl-tRNA synthetases (aaRSs), the rate-determining step in aminoacylation is the dissociation of charged tRNA from the enzyme. In this study, the following factors affecting the release of the charged tRNA from aaRSs are computationally explored: the protonation states of amino acids and substrates present in the active site, and the presence and the absence of AMP and elongation factor Tu.Through molecular modeling, internal pK_a calculations, and molecular dynamics simulations, distinct, mechanistically relevant post-transfer states with charged tRNA bound to glutamyl-tRNA synthetase from Thermus thermophilus (Glu-tRNA^Glu) are considered. The behavior of these nonequilibrium states is characterized as a function of time using dynamical network analysis, local energetics, and changes in free energies to estimate transitions that occur during the release of the tRNA. The hundreds of nanoseconds of simulation time reveal system characteristics that are consistent with recent experimental studies.Energetic and network results support the previously proposed mechanism in which the transfer of amino acid to tRNA is accompanied by the protonation of AMP to H-AMP. Subsequent migration of proton to water reduces the stability of the complex and loosens the interface both in the presence and in the absence of AMP. The subsequent undocking of AMP or tRNA then proceeds along thermodynamically competitive pathways. Release of the tRNA acceptor stem is further accelerated by the deprotonation of the α-ammonium group on the charging amino acid. The proposed general base is Glu41, a residue binding the α-ammonium group that is conserved in both structure and sequence across nearly all class I aaRSs. This universal handle is predicted through pK_a calculations to be part of a proton relay system for destabilizing the bound charging amino acid following aminoacylation. Addition of elongation factor Tu to the aaRS·tRNA complex stimulates the dissociation of the tRNA core and the tRNA acceptor stem.