Selenoprotein synthesis: an expansion of the genetic code

Asymptotic structural properties of quasi-random saturated structures of RNA

Background

RNA folding depends on the distribution of kinetic traps in the landscape of all secondary structures. Kinetic traps in the Nussinov energy model are precisely those secondary structures that are saturated, meaning that no base pair can be added without introducing either a pseudoknot or base triple. In previous work, we investigated asymptotic combinatorics of both random saturated structures and of quasi-random saturated structures, where the latter are constructed by a natural stochastic process.

Results

We prove that for quasi-random saturated structures with the uniform distribution, the asymptotic expected number of external loops is O(logn) and the asymptotic expected maximum stem length is O(logn), while under the Zipf distribution, the asymptotic expected number of external loops is O(log2n) and the asymptotic expected maximum stem length is O(logn/log logn).

Conclusions

Quasi-random saturated structures are generated by a stochastic greedy method, which is simple to implement. Structural features of random saturated structures appear to resemble those of quasi-random saturated structures, and the latter appear to constitute a class for which both the generation of sampled structures as well as a combinatorial investigation of structural features may be simpler to undertake.

Wobble decoding by the Escherichia coli selenocysteine insertion machinery

Selenoprotein expression in Escherichia coli redefines specific single UGA codons from translational termination to selenocysteine (Sec) insertion. This process requires the presence of a Sec Insertion Sequence (SECIS) in the mRNA, which forms a secondary structure that binds a unique Sec-specific elongation factor that catalyzes Sec insertion at the predefined UGA instead of release factor 2-mediated termination. During overproduction of recombinant selenoproteins, this process nonetheless typically results in expression of UGA-truncated products together with the production of recombinant selenoproteins. Here, we found that premature termination can be fully avoided through a SECIS-dependent Sec-mediated suppression of UGG, thereby yielding either tryptophan or Sec insertion without detectable premature truncation. The yield of recombinant selenoprotein produced with this method approached that obtained with a classical UGA codon for Sec insertion. Sec-mediated suppression of UGG thus provides a novel method for selenoprotein production, as here demonstrated with rat thioredoxin reductase. The results also reveal that the E. coli selenoprotein synthesis machinery has the inherent capability to promote wobble decoding.

NADP-specific electron-bifurcating [FeFe]-hydrogenase in a functional complex with formate dehydrogenase in Clostridium autoethanogenum grown on CO

Flavin-based electron bifurcation is a recently discovered mechanism of coupling endergonic to exergonic redox reactions in the cytoplasm of anaerobic bacteria and archaea. Among the five electron-bifurcating enzyme complexes characterized to date, one is a heteromeric ferredoxin- and NAD-dependent [FeFe]-hydrogenase. We report here a novel electron-bifurcating [FeFe]-hydrogenase that is NADP rather than NAD specific and forms a complex with a formate dehydrogenase. The complex was found in high concentrations (6% of the cytoplasmic proteins) in the acetogenic Clostridium autoethanogenum autotrophically grown on CO, which was fermented to acetate, ethanol, and 2,3-butanediol. The purified complex was composed of seven different subunits. As predicted from the sequence of the encoding clustered genes (fdhA/hytA-E) and from chemical analyses, the 78.8-kDa subunit (FdhA) is a selenocysteine- and tungsten-containing formate dehydrogenase, the 65.5-kDa subunit (HytB) is an iron-sulfur flavin mononucleotide protein harboring the NADP binding site, the 51.4-kDa subunit (HytA) is the [FeFe]-hydrogenase proper, and the 18.1-kDa (HytC), 28.6-kDa (HytD), 19.9-kDa (HytE1), and 20.1-kDa (HytE2) subunits are iron-sulfur proteins. The complex catalyzed both the reversible coupled reduction of ferredoxin and NADP(+) with H2 or formate and the reversible formation of H2 and CO2 from formate. We propose the complex to have two functions in vivo, namely, to normally catalyze CO2 reduction to formate with NADPH and reduced ferredoxin in the Wood-Ljungdahl pathway and to catalyze H2 formation from NADPH and reduced ferredoxin when these redox mediators get too reduced during unbalanced growth of C. autoethanogenum on CO (E0' = -520 mV).

Characterization of cytosolic glutathione peroxidase and phospholipid-hydroperoxide glutathione peroxidase genes in rainbow trout (Oncorhynchus mykiss) and their modulation by in vitro selenium exposure

Selenium (Se) is an oligonutrient with both essential biological functions and recognized harmful effects. As the selenocysteine (SeCys) amino acid, selenium is integrated in several Se-containing proteins (selenoproteins), many of which are fundamental for cell homeostasis. Nevertheless, selenium may exert toxic effects at levels marginally above those required, mainly through the generation of reactive oxygen species (ROS). The selenium chemical speciation can strongly affect the bioavailability of this metal and its impact on metabolism, dictating the levels that can be beneficial or detrimental towards an organism. Glutathione peroxidase (GPxs) is the largest and the most studied selenoprotein family. Cytosolic glutathione peroxidase (cGPx, GPx1) and phospholipid hydroperoxide glutathione peroxidase (PHGPx, GPx4) are widely distributed throughout tissues, and play a pivotal role in regulating the oxidative status in the cell. In this study we have cloned GPx1 and GPx4 genes in rainbow trout (Oncorhynchus mykiss). The constitutive mRNA expression of these GPx genes was examined in 18 trout tissues and their responsiveness to Se availability was analysed using a rainbow trout liver cell line (RTL). An inorganic (sodium selenite, Na2SeO3) and organic (selenocysteine, Cys-Se-Se-Cys) selenocompound have been used as Se sources. GPx1 activity was also tested to verify the impact of transcript changes on the enzymatic function of these molecules. To understand if the results obtained from the transcript expression analysis were due to Se bioavailability or generation of ROS, the cytoxicity of the two selenocompounds was tested by measuring the impact of Se on cell membrane integrity. Lastly, Se availability was quantified by mass spectrophotometry to determine the amount of Se in the cell culture media, the Se background due to the foetal calf serum supplement and the contribution from the two selenocompounds used in the treatments. Three isoforms of genes for both GPx1 (GPx1a, 1b1 and 1b2) and GPx4 (GPx4a1, a2 and b) have been identified. The discovery of a third gene encoding for GPx1 and GPx4 hints that salmonids may have the biggest selenoproteome amongst all vertebrates. Transcripts of GPx4 genes were more highly expressed in most tissues examined in vivo (except blood, head kidney and spleen), whereas those of the GPx1 genes were more responsive to selenium exposure in vitro, especially to the organic form. Interestingly, GPx1a was the most sensitive to selenium availability in non stressful conditions, whereas GPx1b1 and GPx1b2 were highly induced by exposure to selenium levels that had some toxic effects on the cells. Although the different concentrations tested of the two selenocompounds modulate GPx1 transcript expression to various degrees, no significant change of GPx1 enzymatic activity was detectable. Our results lead us to conclude that trout GPx1 transcripts expression level may represent a sensitive biomarker for selenium intake, helping to evaluate if selenium concentration and chemical speciation impact on cell homeostasis.

Structure of dipeptides having N-terminal selenocysteine residues: a DFT study in gas and aqueous phase

Over the last few decades, dipeptides as well as their analogues have served as important model systems for the computational studies concerning the structure of protein and energetics of protein folding. Here, we present a density functional structural study on a set of seven dipeptides having N-terminal selenocysteine residues (the component in the C-terminus is varied with seven different combinations viz. Ala, Phe, Glu, Thr, Asn, Arg and Sec) in gas and simulated aqueous phase using a polarizable continuum model (PCM). The molecular geometries of the dipeptides are fully optimized at B3LYP/6-311++G(d,p) level and subsequent frequency calculations confirm them as true minima. The effects of solvation and identity of the varying C-terminal residue on the energetics, structural features of the peptide planes, values of the ψ and ф dihedrals, geometry around the α-carbon atoms and theoretically predicted vibrational spectra of the dipeptides are investigated. Two types of intramolecular H-bonds, namely N…H-N and O…H-C, are found to play important roles in influencing the planarity of the peptide planes and geometry around the α-carbon atoms of the dipeptides. The identity of the varying C-terminal residue influences the values of ф, planarity of the peptide planes and geometry around the C₇ α-carbon atoms while the solvation effects are evident on the values of bond lengths and bond angles of the amide planes.

Selenium. Role of the essential metalloid in health

Selenium is an essential micronutrient in mammals, but is also recognized as toxic in excess. It is a non-metal with properties that are intermediate between the chalcogen elements sulfur and tellurium. Selenium exerts its biological functions through selenoproteins. Selenoproteins contain selenium in the form of the 21st amino acid, selenocysteine (Sec), which is an analog of cysteine with the sulfur-containing side chain replaced by a Se-containing side chain. Sec is encoded by the codon UGA, which is one of three termination codons for mRNA translation in non-selenoprotein genes. Recognition of the UGA codon as a Sec insertion site instead of stop requires a Sec insertion sequence (SECIS) element in selenoprotein mRNAs and a unique selenocysteyl-tRNA, both of which are recognized by specialized protein factors. Unlike the 20 standard amino acids, Sec is biosynthesized from serine on its tRNA. Twenty-five selenoproteins are encoded in the human genome. Most of the selenoprotein genes were discovered by bioinformatics approaches, searching for SECIS elements downstream of in-frame UGA codons. Sec has been described as having stronger nucleophilic and electrophilic properties than cysteine, and Sec is present in the catalytic site of all selenoenzymes. Most selenoproteins, whose functions are known, are involved in redox systems and signaling pathways. However, several selenoproteins are not well characterized in terms of their function. The selenium field has grown dramatically in the last few decades, and research on selenium biology is providing extensive new information regarding its importance for human health.

Selenium-dependent biogenesis of formate dehydrogenase in Campylobacter jejuni is controlled by the fdhTU accessory genes

The food-borne bacterial pathogen Campylobacter jejuni efficiently utilizes organic acids such as lactate and formate for energy production. Formate is rapidly metabolized via the activity of the multisubunit formate dehydrogenase (FDH) enzyme, of which the FdhA subunit is predicted to contain a selenocysteine (SeC) amino acid. In this study we investigated the function of the cj1500 and cj1501 genes of C. jejuni, demonstrate that they are involved in selenium-controlled production of FDH, and propose the names fdhT and fdhU, respectively. Insertional inactivation of fdhT or fdhU in C. jejuni resulted in the absence of FdhA and FdhB protein expression, reduced fdhABC RNA levels, the absence of FDH enzyme activity, and the lack of formate utilization, as assessed by (1)H nuclear magnetic resonance. The fdhABC genes are transcribed from a single promoter located two genes upstream of fdhA, and the decrease in fdhABC RNA levels in the fdhU mutant is mediated at the posttranscriptional level. FDH activity and the ability to utilize formate were restored by genetic complementation with fdhU and by supplementation of the growth media with selenium dioxide. Disruption of SeC synthesis by inactivation of the selA and selB genes also resulted in the absence of FDH activity, which could not be restored by selenium supplementation. Comparative genomic analysis suggests a link between the presence of selA and fdhTU orthologs and the predicted presence of SeC in FdhA. The fdhTU genes encode accessory proteins required for FDH expression and activity in C. jejuni, possibly by contributing to acquisition or utilization of selenium.

Identification and cloning of a selenophosphate synthetase (SPS) from tiger shrimp, Penaeus monodon, and its transcription in relation to molt stages and following pathogen infection

Complementary (c)DNA encoding selenophosphate synthetase (SPS) messenger (m)RNA of the tiger shrimp Penaeus monodon, designated PmSPS, was obtained from the hepatopancreas by a reverse-transcription polymerase chain reaction (RT-PCR) and rapid amplification of cDNA ends (RACE). The 1582-bp cDNA contained an open reading frame (ORF) of 1248 bp, a 103-bp 5'-untranslated region (UTR), and a 231-bp 3'-UTR, which contained a conserved selenocysteine insertion sequence (SECIS) element, a conventional polyadenylation signal, and a poly A tail. The molecular mass of the deduced amino acid (aa) sequence (416 aa) was 45.5 kDa with an estimated pI of 4.85. It contained a putative selenocysteine residue which was encoded by the unusual stop codon, (275)TGA(277), which formed at the active site with residues Sec(58) and Lys(61). A comparison of amino acid sequences showed that PmSPS was more closely related to invertebrate SPS1, such as those of Heliothis virescens and Drosophila melanogaster a and b. PmSPS cDNA was synthesized in all tested tissues, especially in the hepatopancreas. PmSPS in the hepatopancreas of shrimp significantly increased after an injection with either Photobacterium damsela or white spot syndrome virus (WSSV) in order to protect cells against damage from oxidation, and enhance the recycling of selenocysteine or selenium metabolism, indicating that PmSPS is involved in the disease-resistance response. The PmSPS expression by hemocytes significantly increased in stage C, and then gradually decreased until stage A, suggesting that the cloned PmSPS in hemocytes might play a role in viability by renewing hemocytes and antioxidative stress response for new exoskeleton synthesis during the molt cycle of shrimp.

The respiratory molybdo-selenoprotein formate dehydrogenases of Escherichia coli have hydrogen: benzyl viologen oxidoreductase activity

Background

Escherichia coli synthesizes three membrane-bound molybdenum- and selenocysteine-containing formate dehydrogenases, as well as up to four membrane-bound [NiFe]-hydrogenases. Two of the formate dehydrogenases (Fdh-N and Fdh-O) and two of the hydrogenases (Hyd-1 and Hyd-2) have their respective catalytic subunits located in the periplasm and these enzymes have been shown previously to oxidize formate and hydrogen, respectively, and thus function in energy metabolism. Mutants unable to synthesize the [NiFe]-hydrogenases retain a H₂: benzyl viologen oxidoreductase activity. The aim of this study was to identify the enzyme or enzymes responsible for this activity.

Results

Here we report the identification of a new H₂: benzyl viologen oxidoreductase enzyme activity in E. coli that is independent of the [NiFe]-hydrogenases. This enzyme activity was originally identified after non-denaturing polyacrylamide gel electrophoresis and visualization of hydrogen-oxidizing activity by specific staining. Analysis of a crude extract derived from a variety of E. coli mutants unable to synthesize any [NiFe]-hydrogenase-associated enzyme activity revealed that the mutants retained this specific hydrogen-oxidizing activity. Enrichment of this enzyme activity from solubilised membrane fractions of the hydrogenase-negative mutant FTD147 by ion-exchange, hydrophobic interaction and size-exclusion chromatographies followed by mass spectrometric analysis identified the enzymes Fdh-N and Fdh-O. Analysis of defined mutants devoid of selenocysteine biosynthetic capacity or carrying deletions in the genes encoding the catalytic subunits of Fdh-N and Fdh-O demonstrated that both enzymes catalyze hydrogen activation. Fdh-N and Fdh-O can also transfer the electrons derived from oxidation of hydrogen to other redox dyes.

Conclusions

The related respiratory molybdo-selenoproteins Fdh-N and Fdh-O of Escherichia coli have hydrogen-oxidizing activity. These findings demonstrate that the energy-conserving selenium- and molybdenum-dependent formate dehydrogenases Fdh-N and Fdh-O exhibit a degree of promiscuity with respect to the electron donor they use and identify a new class of dihydrogen-oxidizing enzyme.

Disorders of selenium metabolism and selenoprotein function

PURPOSE OF REVIEW

Inborn errors of metabolism are increasingly recognized as underlying causes in pediatric diseases. Selenium and selenoproteins have only recently been identified as causes of inherited defects. Respective case reports have broadened our understanding of selenoprotein function and their developmental importance. This review presents the characterized defects and tries to attract attention to the spectrum of potential phenotypes.

RECENT FINDINGS

The characterization of patients with inherited mutations in selenoprotein N has corroborated the physiological importance of selenium for muscle function. Individuals with inherited defects in selenocysteine insertion sequence (SECIS)-binding protein 2 display a syndrome of selenoprotein-related defects including abnormal thyroid hormone metabolism, delayed bone maturation, and other more individual phenotypes. The recent identification of mutations in selenocysteine synthase causing progressive cerebello-cerebral atrophy underlines the central role of selenoproteins in brain development and protection from neurodegeneration.

SUMMARY

The spectrum of diseases related to inborn defects of selenium utilization, transport, and metabolism is expanding. However, only few examples are already known, resulting from defects in one selenoprotein gene and two genes involved in selenoprotein biosynthesis, respectively. Complex syndromes with impaired muscle function, stunted growth, neurosensory and/or immune defects may point to the involvement of impaired selenium metabolism and selenoprotein function, necessitating specific diagnostic procedures.

Redox biology: computational approaches to the investigation of functional cysteine residues

Cysteine (Cys) residues serve many functions, such as catalysis, stabilization of protein structure through disulfides, metal binding, and regulation of protein function. Cys residues are also subject to numerous post-translational modifications. In recent years, various computational tools aiming at classifying and predicting different functional categories of Cys have been developed, particularly for structural and catalytic Cys. On the other hand, given complexity of the subject, bioinformatics approaches have been less successful for the investigation of regulatory Cys sites. In this review, we introduce different functional categories of Cys residues. For each category, an overview of state-of-the-art bioinformatics methods and tools is provided, along with examples of successful applications and potential limitations associated with each approach. Finally, we discuss Cys-based redox switches, which modify the view of distinct functional categories of Cys in proteins.

Selenoprotein-P is down-regulated in prostate cancer, which results in lack of protection against oxidative damage

BACKGROUND

Oxidative stress plays a role in prostate cancer (PrCa) initiation and development. Selenoprotein-P (SepP; a protein involved in antioxidant defence) mRNA levels are down-regulated in PrCa. The main goal of our study was to assess whether SepP protects prostate cells from reactive oxygen species (ROS) in prostate carcinogenesis.

METHODS

Modification of SepP levels and ROS conditions in C3(1)/Tag-derived cell lines representing prostate epithelial neoplasia (PIN) lesions (Pr-111, with high SepP expression); and invasive tumors (Pr-14, with very low SepP expression).

RESULTS

Both Pr-111 and Pr-14 cells express ApoER2 (SepP receptor), which suggests that they may uptake SepP. Pr-14 cells had much higher ROS levels than Pr-111 cells and were highly sensitive to H(2)O(2)-mediated cytotoxicity. When SepP mRNA levels were knocked down with siRNAs in Pr-111 cells, a significant increase in ROS and cell growth inhibition upon H(2)O(2) exposure was found. Subsequent administration of purified SepP in the culture medium of these cells was able to rescue the original phenotype. Similarly, administration of SepP to Pr-14 cells was able to reduce ROS concentrations. Administration of flutamide decreased SepP mRNA levels whereas dihydrotestosterone or synthetic androgens induced SepP expression, indicating the importance of androgens for SepP expression. Immunohistochemical analysis using a PrCa tissue microarray further revealed that SepP protein was reduced in 60.8% prostate tumors compared to benign prostates.

CONCLUSIONS

Levels of SepP in prostate cells determine basal ROS levels and sensitivity to H(2)O(2)-induced cytotoxicity. Deregulation of SepP during prostate carcinogenesis may increase free radicals, thus promoting tumor development and de-differentiation.

Conformational preferences and pK(a) value of selenocysteine residue

The conformational preferences of the L-selenocysteine (Sec) dipeptides with selenol and selenolate groups (Ac-Sec-NHMe and Ac-Sec(-) -NHMe, respectively) and the apparent (i.e., macroscopic) pK(a) value of the Sec residue have been studied using the dispersion-corrected density functionals M06-2X and B2PLYP-D with the implicit solvation method in the gas phase and in water. In the gas phase, the backbone-to-backbone and/or side chain-to-backbone hydrogen bonds are found to contribute in stabilizing the most preferred conformations for the Sec and Sec(-) residues, as seen for the Cys and Cys(-) residues. However, the polyproline II-like conformations prevail over the conformations with the backbone-to-backbone hydrogen bonds in water because of the weakened hydrogen bonds by the favorable direct interactions between the backbone CO and HN groups and water molecules. The Sec and Sec(-) residues are found to adopt more various conformations than the Cys and Cys(-) residues in water, although the most preferred conformations of the neutral and/or anionic forms of the two residues are similar each other in the gas phase and in water. Using the statistically weighted free energies of the Sec and Sec(-) dipeptides in the gas phase and their solvation free energies, the pK(a) value of the Sec residue is estimated to be 5.47 at 25°C, which is in good agreement with the experimental value of 5.43 ± 0.02. It is found that the lower pK(a) value of the selenol side chain for the Sec residue by ∼3 units than the thiol side chain for the Cys residue is ascribed to the higher gas-phase acidity of the Sec residue.

Thermodynamics of RNA structures by Wang-Landau sampling

Motivation: Thermodynamics-based dynamic programming RNA secondary structure algorithms have been of immense importance in molecular biology, where applications range from the detection of novel selenoproteins using expressed sequence tag (EST) data, to the determination of microRNA genes and their targets. Dynamic programming algorithms have been developed to compute the minimum free energy secondary structure and partition function of a given RNA sequence, the minimum free-energy and partition function for the hybridization of two RNA molecules, etc. However, the applicability of dynamic programming methods depends on disallowing certain types of interactions (pseudoknots, zig-zags, etc.), as their inclusion renders structure prediction an nondeterministic polynomial time (NP)-complete problem. Nevertheless, such interactions have been observed in X-ray structures.

Results: A non-Boltzmannian Monte Carlo algorithm was designed by Wang and Landau to estimate the density of states for complex systems, such as the Ising model, that exhibit a phase transition. In this article, we apply the Wang-Landau (WL) method to compute the density of states for secondary structures of a given RNA sequence, and for hybridizations of two RNA sequences. Our method is shown to be much faster than existent software, such as RNAsubopt. From density of states, we compute the partition function over all secondary structures and over all pseudoknot-free hybridizations. The advantage of the WL method is that by adding a function to evaluate the free energy of arbitary pseudoknotted structures and of arbitrary hybridizations, we can estimate thermodynamic parameters for situations known to be NP-complete. This extension to pseudoknots will be made in the sequel to this article; in contrast, the current article describes the WL algorithm applied to pseudoknot-free secondary structures and hybridizations.

Availability: The WL RNA hybridization web server is under construction at http://bioinformatics.bc.edu/clotelab/.

Contact:clote@bc.edu

Regulation of redox signaling by selenoproteins

The unique chemistry of oxygen has been both a resource and threat for life on Earth for at least the last 2.4 billion years. Reduction of oxygen to water allows extraction of more metabolic energy from organic fuels than is possible through anaerobic glycolysis. On the other hand, partially reduced oxygen can react indiscriminately with biomolecules to cause genetic damage, disease, and even death. Organisms in all three superkingdoms of life have developed elaborate mechanisms to protect against such oxidative damage and to exploit reactive oxygen species as sensors and signals in myriad processes. The sulfur amino acids, cysteine and methionine, are the main targets of reactive oxygen species in proteins. Oxidative modifications to cysteine and methionine can have profound effects on a protein's activity, structure, stability, and subcellular localization. Non-reversible oxidative modifications (oxidative damage) may contribute to molecular, cellular, and organismal aging and serve as signals for repair, removal, or programmed cell death. Reversible oxidation events can function as transient signals of physiological status, extracellular environment, nutrient availability, metabolic state, cell cycle phase, immune function, or sensory stimuli. Because of its chemical similarity to sulfur and stronger nucleophilicity and acidity, selenium is an extremely efficient catalyst of reactions between sulfur and oxygen. Most of the biological activity of selenium is due to selenoproteins containing selenocysteine, the 21st genetically encoded protein amino acid. The most abundant selenoproteins in mammals are the glutathione peroxidases (five to six genes) that reduce hydrogen peroxide and lipid hydroperoxides at the expense of glutathione and serve to limit the strength and duration of reactive oxygen signals. Thioredoxin reductases (three genes) use nicotinamide adenine dinucleotide phosphate to reduce oxidized thioredoxin and its homologs, which regulate a plethora of redox signaling events. Methionine sulfoxide reductase B1 reduces methionine sulfoxide back to methionine using thioredoxin as a reductant. Several selenoproteins in the endoplasmic reticulum are involved in the regulation of protein disulfide formation and unfolded protein response signaling, although their precise biological activities have not been determined. The most widely distributed selenoprotein family in Nature is represented by the highly conserved thioredoxin-like selenoprotein W and its homologs that have not yet been assigned specific biological functions. Recent evidence suggests selenoprotein W and the six other small thioredoxin-like mammalian selenoproteins may serve to transduce hydrogen peroxide signals into regulatory disulfide bonds in specific target proteins.

Hierarchical regulation of selenoprotein expression and sex-specific effects of selenium

The expression of selenoproteins is controlled on each one of the textbook steps of protein biosynthesis, i.e., during gene transcription, RNA processing, translation and posttranslational events as well as via control of the stability of the involved intermediates and final products. Selenoproteins are unique in their dependence on the trace element Se which they incorporate as the 21st proteinogenic amino acid, selenocysteine. Higher mammals have developed unique pathways to enable a fine-tuned expression of all their different selenoproteins according to developmental stage, actual needs, and current availability of the trace element. Tightly controlled and dynamic expression patterns of selenoproteins are present in different tissues. Interestingly, these patterns display some differences in male and female individuals, and can be grossly modified during disease, e.g. in cancer, inflammation or neurodegeneration. Likewise, important health issues related to the selenium status show unexpected sexual dimorphisms. Some detailed molecular insights have recently been gained on how the hierarchical Se distribution among the different tissues is achieved, how the selenoprotein biosynthesis machinery discriminates among the individual selenoprotein transcripts and how impaired selenoprotein biosynthesis machinery becomes phenotypically evident in humans. This review tries to summarize these fascinating findings and highlights some interesting and surprising sex-specific differences.

Probing the role of the proximal heme ligand in cytochrome P450cam by recombinant incorporation of selenocysteine

The unique monooxygenase activity of cytochrome P450cam has been attributed to coordination of a cysteine thiolate to the heme cofactor. To investigate this interaction, we replaced cysteine with the more electron-donating selenocysteine. Good yields of the selenoenzyme were obtained by bacterial expression of an engineered gene containing the requisite UGA codon for selenocysteine and a simplified yet functional selenocysteine insertion sequence (SECIS). The sulfur-to-selenium substitution subtly modulates the structural, electronic, and catalytic properties of the enzyme. Catalytic activity decreases only 2-fold, whereas substrate oxidation becomes partially uncoupled from electron transfer, implying a more complex role for the axial ligand than generally assumed.

How an obscure archaeal gene inspired the discovery of selenocysteine biosynthesis in humans

Selenocysteine (Sec) is the 21st genetically encoded amino acid found in organisms from all three domains of life. Sec biosynthesis is unique in that it always proceeds from an aminoacyl-tRNA precursor. Even though Sec biosynthesis in bacteria was established almost two decades ago, only recently the pathway was elucidated in archaea and eukaryotes. While other aspects of Sec biology have been reviewed previously (Allmang and Krol, Biochimie 2006;88:1561-1571, Hatfield et al., Prog Nucleic Acid Res Mol Biol 2006;81:97-142, Squires and Berry, IUBMB Life 2008;60:232-235), here we review the biochemistry and evolution of Sec biosynthesis and coding and show how the knowledge of an archaeal cysteine biosynthesis pathway helped to uncover the route to Sec formation in archaea and eukaryotes.

Aluminum accumulation induced testicular oxidative stress and altered selenium metabolism in mice

Present work was carried out to investigate how testicular selenium (Se) metabolisms respond to oxidative stress induced by aluminum (Al). Mice were intraperitoneally exposed to 0, 7, or 35mg Al/kg/d for 14 days (CNL, LAL and HAL groups). Al administration significantly increased Al, reactive oxygen radical and malondialdehyde (MDA) levels, as well as decreased glutathione peroxidase (GPx) and glutathione reductase (GR) activities in serum and testes. The serum concentrations of Se were remarkably lower at LAL and HAL groups compared to the controls, whereas the testicular Se levels significantly reduced only in the HAL group. In addition, RT-PCR analysis revealed an increased testicular selenoprotein P (SelP) expression by Al treatment. Western blot analysis showed increased levels of SelP protein expression in the LAL group, but the expression levels were significantly reduced in HAL group. It was suggested that altered metabolism of Se, further stimulated testicular SelP transcription that may compensate for the loss of SelP protein resulted from Al-induced oxidative damage.

Structural modifications into diphenyl diselenide molecule do not cause toxicity in mice

The aim of the present study was to evaluate toxicological parameters of following compounds: 1a (4,4'-dichloro-diphenyl diselenide [(ClPhSe)(2)]), 1b (3,3'-ditrifluoromethyl-diphenyl diselenide [(F(3)CPhSe)(2)]) and 1c (4,4'-dimethoxyl-diphenyl diselenide [(CH(3)OPhSe)(2)]). Calculated lethal dose (LD(50)) values for mice exposed, by oral route, to a single application of compounds 1a, 1b or 1c were estimated to be >381, 278 and >372mg/kg, respectively. Compounds 1a and 1b significantly reduced body weight gain as well as food and water intake in mice. δ-Aminolevulinate dehydratase (δ-ALA-D) and catalase activities were inhibited in mice which received the highest dose of compounds 1a or 1b. Exposure to compounds 1a, 1b and 1c did not modify lipid peroxidation, vitamin C levels, cerebral Na(+)/K(+)-ATPase activity and the biochemical parameters evaluated. The important point for medicinal chemistry is that the structural modifications are not introducing toxicity for the compounds in mice.

Selenoproteins

Selenium is an essential micronutrient for man and animals. The role of selenium has been attributed largely to its presence in selenoproteins as the 21st amino acid, selenocysteine (Sec, U). Sec is encoded by TGA in DNA. A unique mechanism is used to decode the UGA codon in mRNA to co-translationally incorporate Sec into the growing polypeptide because there is no free pool of Sec. In the human genome, 25 genes for selenoproteins have been identified. Selenoproteins such as glutathione peroxidases, thioredoxin reductases, and iodothyronine deiodinases are involved in redox reactions, and Sec is an active-site residue essential for catalytic activity. Selenoproteins have biological functions in oxidoreductions, redox signaling, antioxidant defense, thyroid hormone metabolism, and immune responses. They thus possess a strong correlation with human diseases such as cancer, Keshan disease, virus infections, male infertility, and abnormalities in immune responses and thyroid hormone function.

Evolutionary patterns in the sequence and structure of transfer RNA: a window into early translation and the genetic code

Transfer RNA (tRNA) molecules play vital roles during protein synthesis. Their acceptor arms are aminoacylated with specific amino acid residues while their anticodons delimit codon specificity. The history of these two functions has been generally linked in evolutionary studies of the genetic code. However, these functions could have been differentially recruited as evolutionary signatures were left embedded in tRNA molecules. Here we built phylogenies derived from the sequence and structure of tRNA, we forced taxa into monophyletic groups using constraint analyses, tested competing evolutionary hypotheses, and generated timelines of amino acid charging and codon discovery. Charging of Sec, Tyr, Ser and Leu appeared ancient, while specificities related to Asn, Met, and Arg were derived. The timelines also uncovered an early role of the second and then first codon bases, identified codons for Ala and Pro as the most ancient, and revealed important evolutionary take-overs related to the loss of the long variable arm in tRNA. The lack of correlation between ancestries of amino acid charging and encoding indicated that the separate discoveries of these functions reflected independent histories of recruitment. These histories were probably curbed by co-options and important take-overs during early diversification of the living world.

Introduction to special issue on RNA

In this introduction to the special issue on RNA, we provide a brief overview of some of the novel and exciting biological discoveries concerning diverse roles played by RNA, and subsequently we give a rapid summary of some algorithmic aspects of RNA structure and alignment. Each of the contributions to this special issue is briefly described.

Trends in selenium utilization in marine microbial world revealed through the analysis of the global ocean sampling (GOS) project

Selenium is an important trace element that occurs in proteins in the form of selenocysteine (Sec) and in tRNAs in the form of selenouridine. Recent large-scale metagenomics projects provide an opportunity for understanding global trends in trace element utilization. Herein, we characterized the selenoproteome of the microbial marine community derived from the Global Ocean Sampling (GOS) expedition. More than 3,600 selenoprotein gene sequences belonging to 58 protein families were detected, including sequences representing 7 newly identified selenoprotein families, such as homologs of ferredoxin–thioredoxin reductase and serine protease. In addition, a new eukaryotic selenoprotein family, thiol reductase GILT, was identified. Most GOS selenoprotein families originated from Cys-containing thiol oxidoreductases. In both Pacific and Atlantic microbial communities, SelW-like and SelD were the most widespread selenoproteins. Geographic location had little influence on Sec utilization as measured by selenoprotein variety and the number of selenoprotein genes detected; however, both higher temperature and marine (as opposed to freshwater and other aquatic) environment were associated with increased use of this amino acid. Selenoproteins were also detected with preference for either environment. We identified novel fusion forms of several selenoproteins that highlight redox activities of these proteins. Almost half of Cys-containing SelDs were fused with NADH dehydrogenase, whereas such SelD forms were rare in terrestrial organisms. The selenouridine utilization trait was also analyzed and showed an independent evolutionary relationship with Sec utilization. Overall, our study provides insights into global trends in microbial selenium utilization in marine environments.

Selenium (Se) is an essential micronutrient due to its requirement for biosynthesis and function of the 21^st amino acid, selenocysteine (Sec). Sec is found in the active sites of selenoproteins, most of which exhibit redox function, in all three domains of life. In recent years, genome sequencing projects provided a large volume of nucleotide and protein sequence information. Identification of complete sets of selenoproteins (selenoproteomes) of individual organisms and environmental samples is important for better understanding of Se utilization, biological functions of this element, and changes in Se use during evolution. Here, we describe a comprehensive analysis of the selenoproteome of the microbial marine community derived from the Global Ocean Sampling (GOS) expedition. More than 3,600 selenoprotein gene sequences belonging to 58 protein families were detected and analyzed. Our study generated the largest selenoproteome reported to date and provided important insights into microbial Se utilization and its evolutionary trends in marine environments.

Towards understanding selenocysteine incorporation into bacterial proteins

In bacteria, UGA stop codons can be recoded to direct the incorporation of selenocysteine into proteins on the ribosome. Recoding requires a selenocysteine incorporation sequence (SECIS) downstream of the UGA codon, a specialized translation factor SelB, and the non-canonical Sec-tRNASec, which is formed from Ser-tRNASec by selenocysteine synthase, SelA, using selenophosphate as selenium donor. Here we describe a rapid-kinetics approach to study the mechanism of selenocysteine insertion into proteins on the ribosome. Labeling of SelB, Sec-tRNASec and other components of the translational machinery allows direct observation of the formation or dissociation of complexes by monitoring changes in the fluorescence of single dyes or fluorescence resonance energy transfer between two fluorophores. Furthermore, the structure of SelA was studied by electron cryomicroscopy (cryo-EM). We report that intact SelA from the thermophilic bacterium Moorella thermoacetica (mthSelA) can be vitrified for cryo-EM using a controlled-environment vitrification system. Two-dimensional image analysis of vitrified mthSelA images shows that SelA can adopt the wide range of orientations required for high-resolution structure determination by cryo-EM. The results indicate that mthSelA forms a homodecamer that has a ring-like structure with five bilobed wings, similar to the structure of the E. coli complex determined previously.

Molecular Biology

Genetic information is passed with high accuracy from the parental organism to the offspring and its expression governs the biochemical and physiological tasks of the cell. Although different types of cells exist and are shaped by development to fill different physiological niches, all cells have fundamental similarities and share common principles of organization and biochemical activities. This chapter gives an overview of general principles of the storage and flow of genetic information. It aims to summarize and describe in a broadly approachable way, from the point of view of molecular biology, some general terms, mechanisms and processes used as a base for the molecular computing in the subsequent chapters.

Transfer RNA and the origins of diversified life

The evolution of the transfer RNA (tRNA) molecule is controversial but embeds the history of protein biosynthesis, the genetic code, and the origins of diversified life. A new phylogenetic method based on RNA structure that we developed provides new lines of evidence to support the genome tag hypothesis and confirms that the 'top half' of tRNA is more ancient than the 'bottom half'. Timelines of amino acid charging function generated from constraint analyses showed that selenocysteine, tyrosine, serine, and leucine specificities were ancient, while those related to asparagine, methionine, and arginine were more recent. The timelines also uncovered an early role of the second and then first codon bases, identified codons for alanine and proline as the most ancient, and revealed important evolutionary take-overs related to the loss of the long variable arm of tRNA. Furthermore, organismal timelines showed Archaea was the oldest superkingdom, followed by viruses, and superkingdoms Eukarya and Bacteria in that order supporting conclusions from recent phylogenomic studies of protein architecture. Strikingly, results showed that the origin of viruses was not only ancient but was linked to Archaea, supporting the notion that the archaeal lineage is the most ancient on earth and its origin predated diversification of tRNA function and specificity.

Structure and catalytic mechanism of eukaryotic selenocysteine synthase

In eukaryotes and Archaea, selenocysteine synthase (SecS) converts O-phospho-L-seryl-tRNA [Ser]Sec into selenocysteyl-tRNA [Ser]Sec using selenophosphate as the selenium donor compound. The molecular mechanisms underlying SecS activity are presently unknown. We have delineated a 450-residue core of mouse SecS, which retained full selenocysteyl-tRNA [Ser]Sec synthesis activity, and determined its crystal structure at 1.65 A resolution. SecS exhibits three domains that place it in the fold type I family of pyridoxal phosphate (PLP)-dependent enzymes. Two SecS monomers interact intimately and together build up two identical active sites around PLP in a Schiff-base linkage with lysine 284. Two SecS dimers further associate to form a homotetramer. The N terminus, which mediates tetramer formation, and a large insertion that remodels the active site set SecS aside from other members of the family. The active site insertion contributes to PLP binding and positions a glutamate next to the PLP, where it could repel substrates with a free alpha-carboxyl group, suggesting why SecS does not act on free O-phospho-l-serine. Upon soaking crystals in phosphate buffer, a previously disordered loop within the active site insertion contracted to form a phosphate binding site. Residues that are strictly conserved in SecS orthologs but variant in related enzymes coordinate the phosphate and upon mutation corrupt SecS activity. Modeling suggested that the phosphate loop accommodates the gamma-phosphate moiety of O-phospho-l-seryl-tRNA [Ser]Sec and, after phosphate elimination, binds selenophosphate to initiate attack on the proposed aminoacrylyl-tRNA [Ser]Sec intermediate. Based on these results and on the activity profiles of mechanism-based inhibitors, we offer a detailed reaction mechanism for the enzyme.

Minimalist active-site redesign: teaching old enzymes new tricks

Although nature evolves its catalysts over millions of years, enzyme engineers try to do it a bit faster. Enzyme active sites provide highly optimized microenvironments for the catalysis of biologically useful chemical transformations. Consequently, changes at these centers can have large effects on enzyme activity. The prediction and control of these effects provides a promising way to access new functions. The development of methods and strategies to explore the untapped catalytic potential of natural enzyme scaffolds has been pushed by the increasing demand for industrial biocatalysts. This Review describes the use of minimal modifications at enzyme active sites to expand their catalytic repertoires, including targeted mutagenesis and the addition of new reactive functionalities. Often, a novel activity can be obtained with only a single point mutation. The many successful examples of active-site engineering through minimal mutations give useful insights into enzyme evolution and open new avenues in biocatalyst research.

A recoding element that stimulates decoding of UGA codons by Sec tRNA[Ser]Sec

Selenocysteine insertion during decoding of eukaryotic selenoprotein mRNA requires several trans-acting factors and a cis-acting selenocysteine insertion sequence (SECIS) usually located in the 3' UTR. A second cis-acting selenocysteine codon redefinition element (SRE) has recently been described that resides near the UGA-Sec codon of selenoprotein N (SEPN1). Similar phylogenetically conserved elements can be predicted in a subset of eukaryotic selenoprotein mRNAs. Previous experimental analysis of the SEPN1 SRE revealed it to have a stimulatory effect on readthrough of the UGA-Sec codon, which was not dependent upon the presence of a SECIS element in the 3' UTR; although, as expected, readthrough efficiency was further elevated by inclusion of a SECIS. In order to examine the nature of the redefinition event stimulated by the SEPN1 SRE, we have modified an experimentally tractable in vitro translation system that recapitulates efficient selenocysteine insertion. The results presented here illustrate that the SRE element has a stimulatory effect on decoding of the UGA-Sec codon by both the methylated and unmethylated isoforms of Sec tRNA([Ser]Sec), and confirm that efficient selenocysteine insertion is dependent on the presence of a 3'-UTR SECIS. The variation in recoding elements predicted near UGA-Sec codons implies that these elements may play a differential role in determining the amount of selenoprotein produced by acting as controllers of UGA decoding efficiency.

Natural history and experimental evolution of the genetic code

The standard genetic code is a set of rules that relates the 20 canonical amino acids in proteins to groups of three bases in the mRNA. It evolved from a more primitive form and the attempts to reconstruct its natural history are based on its present-day features. Genetic code engineering as a new research field was developed independently in a few laboratories during the last 15 years. The main intention is to re-program protein synthesis by expanding the coding capacities of the genetic code via re-assignment of specific codons to un-natural amino acids. This article focuses on the question as to which extent hypothetical scenarios that led to codon re-assignments during the evolution of the genetic code are relevant for its further evolution in the laboratory. Current attempts to engineer the genetic code are reviewed with reference to theoretical works on its natural history. Integration of the theoretical considerations into experimental concepts will bring us closer to designer cells with target-engineered genetic codes that should open not only tremendous possibilities for the biotechnology of the twenty-first century but will also provide a basis for the design of novel life forms.

Computing the Partition Function and Sampling for Saturated Secondary Structures of RNA, with Respect to the Turner Energy Model

An RNA secondary structure is saturated if no base pairs can be added without violating the definition of secondary structure. Here we describe a new algorithm, RNAsat, which for a given RNA sequence a, an integral temperature 0 <or= T <or= 100 in degrees Celsius, and for all integers k, computes the Boltzmann partition function Z(k)(T)(a) = SigmaSepsilonSAT(k)(a) exp(-E(S)/RT), where the sum is over all saturated secondary structures of a which have exactly k base pairs, R is the universal gas constant and E(S) denotes the free energy with respect to the Turner nearest neighbor energy model. By dynamic programming, we compute Z(k)(T)simultaneously for all values of k in time O(n(5)) and space O(n(3)).Additionally, RNAsat computes the partition function Q(k)(T)(a) = SigmaSepsilonS(k)(a) exp(-E(S)/RT), where the sum is over all secondary structures of a which have k base pairs; the latter computation is performed simultaneously for all values of k in O(n(4)) time and O(n(3)) space. Lastly, using the partition function Z(k)(T) [resp. Q(k)(T)] with stochastic backtracking, RNAsat rigorously samples the collection of saturated secondary structures [resp. secondary structures] having k base pairs; for Q(k)(T) this provides a parametrized form of Sfold sampling (Ding and Lawrence, 2003). Using RNAsat, (i) we compute the ensemble free energy for saturated secondary structures having k base pairs, (ii) show cooperativity of the Turner model, (iii) demonstrate a temperature-dependent phase transition, (iv) illustrate the predictive advantage of RNAsat for precursor microRNA cel-mir-72 of C. elegans and for the pseudoknot PKB 00152 of Pseudobase (van Batenburg et al., 2001), (v) illustrate the RNA shapes (Giegerich et al., 2004) of sampled secondary structures [resp. saturated structures] having exactly k base pairs. A web server for RNAsat is under construction at bioinformatics.bc.edu/clotelab/RNAsat/.

Neuronal and ependymal expression of selenoprotein P in the human brain

Selenoprotein P (SePP) is central to selenium (Se) metabolism in the mammalian organism. Human SePP contains 10 Se atoms that are covalent constituents of the polypeptide chain incorporated as the rare amino acid selenocysteine (Sec). Since hepatocytes secrete SePP into plasma, SePP is commonly regarded as a Se transport protein, although SePP mRNA is expressed in many organs. Gene targeting of SePP in mice leads to neurological dysfunction resulting from Se deficiency and associated reduction of selenoenzyme activities in the brain. However, more recent data revealed that isolated hepatic SePP deficiency does not alter brain Se levels, suggesting a role for SePP locally expressed in the brain. Some of the best characterized and most abundant selenoenzymes, glutathione peroxidases, thioredoxin reductases, and methionine sulfoxide reductase B, play major roles in the cellular defense against reactive oxygen species. Therefore, it was hypothesized that reduced brain Se bioavailability may be involved in the pathogenesis of neurodegenerative disease and normal ageing. We present evidence that human CSF contains SePP and that the human brain expresses SePP mRNA. Moreover, SePP-like immunoreactivity localizes to neurons and ependymal cells and thus appears strategically situated for maintenance and control of Se-dependent anti-oxidative defense systems.

Combinatorics of saturated secondary structures of RNA

Following Zuker (1986), a saturated secondary structure for a given RNA sequence is a secondary structure such that no base pair can be added without violating the definition of secondary structure, e.g., without introducing a pseudoknot. In the Nussinov-Jacobson energy model (Nussinov and Jacobson, 1980), where the energy of a secondary structure is -1 times the number of base pairs, saturated secondary structures are local minima in the energy landscape, hence form kinetic traps during the folding process. Here we present recurrence relations and closed form asymptotic limits for combinatorial problems related to the number of saturated secondary structures. In addition, Python source code to compute the number of saturated secondary structures having k base pairs can be found at the web servers link of bioinformatics.bc.edu/clotelab/.

Identification of Leishmania selenoproteins and SECIS element

Selenoproteins result from the incorporation of selenocysteine (Sec-U) at an UGA-stop codon positioned within a gene's open reading frame and directed by selenocysteine insertion sequence (SECIS) elements. Although the selenocysteine incorporation pathway has been identified in a wide range of organisms it has not yet been reported in the Kinetoplastida Leishmania and Trypanosoma. Here we present evidence consistent with the presence of a selenocysteine biosynthetic pathway in Kinetoplastida. These include the existence of SECIS-containing coding sequences in Leishmania major and Leishmania infantum, the incorporation of (75)Se into Leishmania proteins, the occurrence of selenocysteine-tRNA (tRNA (sec) (uca)) in both Leishmania and Trypanosoma and in addition the finding of all genes necessary for selenocysteine synthesis such as SELB, SELD, PSTK and SECp43. As in other eukaryotes, the Kinetoplastids have no identifiable SELA homologue. To our knowledge this is the first report on the identification of selenocysteine insertion machinery in Kinetoplastida, more specifically in Leishmania, at the sequence level.