Optimal Bayesian Design for Patient Selection in a Clinical Study

Summary .  Bayesian experimental design for a clinical trial involves specifying a utility function that models the purpose of the trial, in this case the selection of patients for a diagnostic test. The best sample of patients is selected by maximizing expected utility. This optimization task poses difficulties due to a high-dimensional discrete design space and, also, to an expected utility formula of high complexity. A simulation-based optimal design method is feasible in this case. In addition, two deterministic algorithms that perform a systematic search over the design space are developed to address the computational issues.     

Analysis of recognition of transfer-messenger RNA by the ribosomal decoding center

Bacterial ribosomes stalled on defective mRNAs are rescued by tmRNA that functions as both tRNA and mRNA. The first ribosomal elongation cycle on tmRNA where tmRNA functions as tRNA is highly unusual: occupation of the ribosomal A site by tmRNA occurs without codon:anticodon pairing. Our analysis shows that in this case the role of a codon:anticodon duplex should be accomplished by a single unpaired triplet. In order that tmRNA could participate in the ribosomal elongation cycle, a triplet preceding the mRNA portion of tmRNA (the -1triplet) should be in the A-form and this form should be recognized by the ribosomal decoding center. A rule is derived that determines what triplets cannot be used as the -1triplet. The rule was tested with the -1triplets of all known 414 tmRNA species. All 23 observed -1triplets follow the formulated rule. The rule is also supported by the available data on base substitutions within the -1triplet.     

Relationships Among Stop Codon Usage Bias, Its Context, Isochores, and Gene Expression Level in Various Eukaryotes

It is well known that stop codons play a critical role in the process of protein synthesis. However, little effort has been made to investigate whether stop codon usage exhibits biases, such as widely seen for synonymous codon usage. Here we systematically investigate stop codon usage bias in various eukaryotes as well as its relationships with its context, GC3 content, gene expression level, and secondary structure. The results show that there is a strong bias for stop codon usage in different eukaryotes, i.e., UAA is overrepresented in the lower eukaryotes, UGA is overrepresented in the higher eukaryotes, and UAG is least used in all eukaryotes. Different conserved patterns for each stop codon in different eukaryotic classes are found based on information content and logo analysis. GC3 contents increase with increasing complexity of organisms. Secondary structure prediction revealed that UAA is generally associated with loop structures, whereas UGA is more uniformly present in loop and stem structures, i.e., UGA is less biased toward having a particular structure. The stop codon usage bias, however, shows no significant relationship with GC3 content and gene expression level in individual eukaryotes. The results indicate that genomic complexity and GC3 content might contribute to stop codon usage bias in different eukaryotes. Our results indicate that stop codons, like synonymous codons, exhibit biases in usage. Additional work will be needed to understand the causes of these biases and their relationship to the mechanism of protein termination. [Reviewing Editor: Dr. Manyuan Long]     

Intragenic Hill-Robertson interference influences selection intensity on synonymous mutations in Drosophila

Natural selection influences synonymous mutations and synonymouscodon usage in many eukaryotes to improve the efficiency oftranslation in highly expressed genes. Recent studies of genecomposition in eukaryotes have shown that codon usage also variesindependently of expression levels, both among genes and atthe intragenic level. Here, we investigate rates of evolution(Ks) and intensity of selection (_s) on synonymous mutationsin two groups of genes that differ greatly in the length oftheir exons, but with equivalent levels of gene expression andrates of crossing-over in Drosophila melanogaster. We estimate_s using patterns of divergence and polymorphism in 50 Drosophilagenes (100 kb of coding sequence) to take into account possiblevariation in mutation trends across the genome, among genesor among codons. We show that genes with long exons exhibithigher Ks and reduced _s compared to genes with short exons.We also show that Ks and _s vary significantly across long exons,with higher Ks and reduced _s in the central region comparedto flanking regions of the same exons, hence indicating thatthe difference between genes with short and long exons can bemostly attributed to the central region of these long exons.Although amino acid composition can also play a significantrole when estimating Ks and _s, our analyses show that the differencesin Ks and _s between genes with short and long exons and acrosslong exons cannot be explained by differences in protein composition.All these results are consistent with the Interference Selection(IS) model that proposes that the Hill-Robertson (HR) effectcaused by many weakly selected mutations has detectable evolutionaryconsequences at the intragenic level in genomes with recombination.Under the IS model, exon size and exon-intron structure influencethe effectiveness of selection, with long exons showing reducedeffectiveness of selection when compared to small exons andthe central region of long exons showing reduced intensity ofselection compared to flanking coding regions. Finally, ourresults further stress the need to consider selection on synonymousmutations and its variation—among and across genes andexons—in studies of protein evolution.     

A -1 frameshift in the HIV-1 env gene is enhanced by arginine deficiency via a hungry codon mechanism

Ribosomal frameshifting is used by various organisms to maximize protein coding potential of genomic sequences. It is commonly exploited by RNA viruses to overcome the constraint of their limited genome size. Frameshifting requires specific RNA structural features, such as a suitable heptanucleotide “slippery” sequence and an RNA pseudoknot. Previous genomic analysis of HIV-1 indicated the potential for several hidden genes encoded through frameshifting; one of these, overlapping the envelope gene, has an RNA pseudoknot just downstream from a slippery sequence, AAAAAGA that features an adenine quadruplet prior to a potential hungry arginine codon (AGA). This env-frameshift (env-fs) gene has been shown to encode a truncated glutathione peroxidase homologue, with both antioxidant and anti-apoptotic activities in transfected cells. Using a dual reporter cell-based frameshift assay, we demonstrate that the env-fs frameshift sequence is active in vitro. Furthermore, in arginine deficient media, env-fs frameshifting increased over 100% (p < 0.005), consistent with the hypothesized hungry codon mechanism. As a response to arginine deficiency, increased expression of the antioxidant viral GPx gene (env-fs) by upregulation of frameshifting could be protective to HIV-infected cells, as a countermeasure to the increased oxidative stress induced by arginine deficiency (because NO is a known scavenger of hydroxyl radical).     

The evolution of repeated mating under sexual conflict

In insects, repeated mating by females may have direct effects on female fecundity, fertility, and longevity. In addition, a female's remating rate affects her fitness through mortality costs of male harassment and ecological risks of mating such as predation. We analyse a model where these female fitness factors are put into their life-history context, and traded against each other, while accounting for limitations because of mate availability. We solve analytically for the condition when female multiple mating will evolve. We show that the probability that a female mates with a courting male decreases with increases in population density. The extent of conflict between the sexes thus automatically becomes larger at higher densities. However, because at higher densities females meet males at a higher rate, the resulting ESS female remating rate is independent of population density. The female remating probability is in conflict with male adaptations that increase male mating rate by persuading or forcing females to mate, and also in conflict with male adaptations for protecting the own sperm from being removed by future female mates. We show that the relative importance of these conflicts depends on population density.     

High level expression of a recombinant acid phytase gene in Pichia pastoris

AIMS: To achieve high phytase yield with improved enzymatic activity in Pichia pastoris. METHODS AND RESULTS: The 1347-bp phytase gene of Aspergillus niger SK-57 was synthesized using a successive polymerase chain reaction and was altered by deleting intronic sequences, optimizing codon usage and replacing its original signal sequence with a synthetic signal peptide (designated MF4I) that is a codon-modified Saccharomyces cerevisiae mating factor alpha-prepro-leader sequence. The gene constructs containing wild type or modified phytase gene coding sequences under the control of the highly-inducible alcohol oxidase gene promoter with the MF4I- or wild type alpha-signal sequence were used to transform Pichia pastoris. The P. pastoris strain that expressed the modified phytase gene (phyA-sh) with MF4I sequence produced 6.1 g purified phytase per litre of culture fluid, with the phytase activity of 865 U ml(-1). The expressed phytase varied in size (64, 67, 87, 110 and 120 kDa), but could be deglycosylated to produce a homogeneous 64 kDa protein. The recombinant phytase had two pH optima (pH 2.5 and pH 5.5) and an optimum temperature of 60 degrees C. CONCLUSIONS: The P. pastoris strain with the genetically engineered phytase gene produced 6.1 g l(-1) of phytase or 865 U ml(-1) phytase activity, a 14.5-fold increase compared with the P. pastoris strain with the wild type phytase gene. SIGNIFICANCE AND IMPACT OF THE STUDY: The P. pastoris strain expressing the modified phytase gene with the MF4I signal peptide showed great potential as a commercial phytase production system.     

Effect of energy source on the efficiency of translational termination during cell-free protein synthesis

We studied how the fidelity of translation termination is affected by the method of ATP regeneration during cell-free protein synthesis. During the in vivo expression of hEPO, whose termination is directed by the UGA codon, we found that substantial proportions of the translational products showed a larger molecular weight than expected. Similar results were obtained in a cell-free synthesis reaction using phosphoenol pyruvate (PEP) or 3-phosphoglycerate (3PG) for ATP regeneration. However, when the energy source was switched to creatine phosphate (CP), the readthrough of the UGA codon was completely repressed and only the target protein of the correct size was expressed in a high yield. To the best of our knowledge, this is the first report describing the relationship between the regeneration of nucleotide triphosphates and protein readthrough, and we also believe that the discovery would pave the way to the selective and efficient expression of target proteins in cell-free protein synthesis systems.