Avoidance of Long Mononucleotide Repeats in Codon Pair Usage

Protein is an essential component for life, and its synthesis is mediated by codons in any organisms on earth. While some codons encode the same amino acid, their usage is often highly biased. There are many factors that can cause the bias, but a potential effect of mononucleotide repeats, which are known to be highly mutable, on codon usage and codon pair preference is largely unknown. In this study we performed a genomic survey on the relationship between mononucleotide repeats and codon pair bias in 53 bacteria, 68 archaea, and 13 eukaryotes. By distinguishing the codon pair bias from the codon usage bias, four general patterns were revealed: strong avoidance of five or six mononucleotide repeats in codon pairs; lower observed/expected (o/e) ratio for codon pairs with C or G repeats (C/G pairs) than that with A or T repeats (A/T pairs); a negative correlation between genomic GC contents and the o/e ratios, particularly for C/G pairs; and avoidance of C/G pairs in highly conserved genes. These results support natural selection against long mononucleotide repeats, which could induce frameshift mutations in coding sequences. The fact that these patterns are found in all kingdoms of life suggests that this is a general phenomenon in living organisms. Thus, long mononucleotide repeats may play an important role in base composition and genetic stability of a gene and gene functions.

A novel computer simulation method for simulating the multiscale transduction dynamics of signal proteins

Signal proteins are able to adapt their response to a change in the environment, governing in this way a broad variety of important cellular processes in living systems. While conventional molecular-dynamics (MD) techniques can be used to explore the early signaling pathway of these protein systems at atomistic resolution, the high computational costs limit their usefulness for the elucidation of the multiscale transduction dynamics of most signaling processes, occurring on experimental timescales. To cope with the problem, we present in this paper a novel multiscale-modeling method, based on a combination of the kinetic Monte-Carlo- and MD-technique, and demonstrate its suitability for investigating the signaling behavior of the photoswitch light-oxygen-voltage-2-Jα domain from Avena Sativa (AsLOV2-Jα) and an AsLOV2-Jα-regulated photoactivable Rac1-GTPase (PA-Rac1), recently employed to control the motility of cancer cells through light stimulus. More specifically, we show that their signaling pathways begin with a residual re-arrangement and subsequent H-bond formation of amino acids near to the flavin-mononucleotide chromophore, causing a coupling between β-strands and subsequent detachment of a peripheral α-helix from the AsLOV2-domain. In the case of the PA-Rac1 system we find that this latter process induces the release of the AsLOV2-inhibitor from the switchII-activation site of the GTPase, enabling signal activation through effector-protein binding. These applications demonstrate that our approach reliably reproduces the signaling pathways of complex signal proteins, ranging from nanoseconds up to seconds at affordable computational costs.

Clinical, structural and functional implications of mutations and polymorphisms in human NADPH P450 oxidoreductase

Cytochrome P450 proteins are involved in metabolism of drugs and xenobiotics. In the endoplasmic reticulum a single nicotinamide adenine dinucleotide phosphate (NADPH) P450 oxidoreductase (POR) supplies electrons to all microsomal P450s for catalytic activity. POR is a flavoprotein that contains both flavin mononucleotide and flavin adenine dinucleotide as cofactors and uses NADPH as the source of electrons. We have recently reported a number of POR mutations in the patients with disordered steroidogenesis. In the first report we had described missense mutations (A287P, R457H, V492E, C569Y, and V608F) identified in four patients with defects in steroid production. Each POR variant was produced as recombinant N-27 form of the enzyme in bacteria and as full-length form in yeast. Membranes from bacteria or yeast expressing normal or variant POR were purified and their activities were characterized in cytochrome c and CYP17A1 assays. Later we have published a larger study that described a whole range of POR mutations and characterized the mutants/polymorphisms A115V, T142A, M263V, Y459H, A503V, G539R, L565P, R616X, V631I, and F646del from the sequencing of patient DNA. We also studied POR variants Y181D, P228L, R316W, G413S, and G504R that were available in public databases or published literature. Three-dimensional structure of rat POR is known and we have used this structure to deduce the structure-function correlation of POR mutations in human. The missense mutations found in patients with disordered steroidogenesis are generally in the co-factor binding and functionally important domains of POR and the apparent polymorphisms are found in regions with lesser structural importance. A variation in POR can alter the activity of all microsomal P450s, and therefore, can affect the metabolism of drugs and xenobiotics even when the P450s involved are otherwise normal. It is important to study the genetic and biochemical basis of POR variants in human population to gain information about possible differences in P450 mediated reactions among the individuals carrying a variant or polymorphic form of POR that could impact their metabolism.