Kei Yura

Amino acid residue doublet propensity in the protein–RNA interface and its application to RNA interface prediction

Protein–RNA interactions play essential roles in a number of regulatory mechanisms for gene expression such as RNA splicing, transport, translation and post-transcriptional control. As the number of available protein–RNA complex 3D structures has increased, it is now possible to statistically examine protein–RNA interactions based on 3D structures. We performed computational analyses of 86 representative protein–RNA complexes retrieved from the Protein Data Bank. Interface residue propensity, a measure of the relative importance of different amino acid residues in the RNA interface, was calculated for each amino acid residue type (residue singlet interface propensity). In addition to the residue singlet propensity, we introduce a new residue-based propensity, which gives a measure of residue pairing preferences in the RNA interface of a protein (residue doublet interface propensity). The residue doublet interface propensity contains much more information than the sum of two singlet propensities alone. The prediction of the RNA interface using the two types of propensities plus a position-specific multiple sequence profile can achieve a specificity of about 80%. The prediction method was then applied to the 3D structure of two mRNA export factors, TAP (Mex67) and UAP56 (Sub2). The prediction enables us to point out candidate RNA interfaces, part of which are consistent with previous experimental studies and may contribute to elucidation of atomic mechanisms of mRNA export.

Novel types of two-domain multi-copper oxidases: possible missing links in the evolution.

An analysis of the genome sequence database revealed novel types of two‐domain multi‐copper oxidases. The two‐domain proteins have the conspicuous combination of blue‐copper and inter‐domain trinuclear copper binding residues, which is common in ceruloplasmin and ascorbate oxidase but not in nitrite reductase, and therefore are considered to retain the characteristics of the plausible ancestral form of ceruloplasmin and ascorbate oxidase. A possible evolutionary relationship of these proteins is proposed.

Human origin recognition complex binds preferentially to G-quadruplex-preferable RNA and single-stranded DNA

Background: ORC binds to replication origins, but human ORC does not exhibit apparent sequence-specificity. Results: G-quadruplex (G4)-preferable RNA or single-stranded DNA competes for DNA binding of ORC. Conclusion: Human ORC binds preferentially to RNA and single-stranded DNA that form G4, and the certain domain in ORC1 is involved in this binding. Significance: This ability may correlate with the G4-formable motif in human replication origins. Origin recognition complex (ORC), consisting of six subunits ORC1–6, is known to bind to replication origins and function in the initiation of DNA replication in eukaryotic cells. In contrast to the fact that Saccharomyces cerevisiae ORC recognizes the replication origin in a sequence-specific manner, metazoan ORC has not exhibited strict sequence-specificity for DNA binding. Here we report that human ORC binds preferentially to G-quadruplex (G4)-preferable G-rich RNA or single-stranded DNA (ssDNA). We mapped the G-rich RNA-binding domain in the ORC1 subunit, in a region adjacent to its ATPase domain. This domain itself has an ability to preferentially recognize G4-preferable sequences of ssDNA. Furthermore, we found, by structure modeling, that the G-rich RNA-binding domain is similar to the N-terminal portion of AdoMet_MTase domain of mammalian DNA methyltransferase 1. Therefore, in contrast with the binding to double-stranded DNA, human ORC has an apparent sequence preference with respect to its RNA/ssDNA binding. Interestingly, this specificity coincides with the common signature present in most of the human replication origins. We expect that our findings provide new insights into the regulations of function and chromatin binding of metazoan ORCs.

Essential role of the Arg112 residue of cytochrome P450cam for electron transfer from reduced putidaredoxin

Cytochrome P450cam (CYP101) of Pseudomonas putida PpGl in which Arg112 is substituted by Cys was isolated by in vitro random mutagenesis of the camC gene DNA coding for P450cam. The absorption spectra of the purified mutant enzyme were similar to those of the wild type enzyme, but its substrate‐dependent NADH oxidation activity in the presence of putidaredoxin (Pd) and putidaredoxin reductase (PdR) was extremely low. The rate constant of electron transfer from reduced Pd to the heme of the mutant P450cam, measured on an anaerobic stopped flow apparatus, was 1/400 of that of the wild type enzyme and the dissociation constant of the mutant P450cam for oxidized Pd was several fold higher than that of the wild type enzyme. A considerable decrease in mid‐point potential of the mutant enzyme was also noted. We conclude that Arg112, which is located on the surface of the P450cam molecule and hydrogen‐bonded to one of the heme propionate chains, plays an essential role in the electron transfer from Pd.

Large-scale identification and characterization of alternative splicing variants of human gene transcripts using 56,419 completely sequenced and manually annotated full-length cDNAs.

We report the first genome-wide identification and characterization of alternative splicing in human gene transcripts based on analysis of the full-length cDNAs. Applying both manual and computational analyses for 56 419 completely sequenced and precisely annotated full-length cDNAs selected for the H-Invitational human transcriptome annotation meetings, we identified 6877 alternative splicing genes with 18 297 different alternative splicing variants. A total of 37 670 exons were involved in these alternative splicing events. The encoded protein sequences were affected in 6005 of the 6877 genes. Notably, alternative splicing affected protein motifs in 3015 genes, subcellular localizations in 2982 genes and transmembrane domains in 1348 genes. We also identified interesting patterns of alternative splicing, in which two distinct genes seemed to be bridged, nested or having overlapping protein coding sequences (CDSs) of different reading frames (multiple CDS). In these cases, completely unrelated proteins are encoded by a single locus. Genome-wide annotations of alternative splicing, relying on full-length cDNAs, should lay firm groundwork for exploring in detail the diversification of protein function, which is mediated by the fast expanding universe of alternative splicing variants.

Tuning glycosidase inhibition through aglycone interactions: pharmacological chaperones for Fabry disease and GM1 gangliosidosis.

Competitive inhibitors of either α-galactosidase (α-Gal) or β-galactosidase (β-Gal) with high affinity and selectivity have been accessed by exploiting aglycone interactions with conformationally locked sp(2)-iminosugars. Selected compounds were profiled as potent pharmacological chaperones for mutant lysosomal α- and β-Gal associated with Fabry disease and GM(1) gangliosidosis.

Module-intron correlation and intron sliding in family F/10 xylanase genes

Xylanases are classified into two families, numbered F/10 and G/11 according to the similarity of amino acid sequences of their catalytic domain (Henrissat, B., Bairoch, A., 1993. New families in the classification of glycosyl hydrolases based on amino acid sequence similarities. Biochem. J. 293, 781-788). Three-dimensional structure of the catalytic domain of the family F/10 xylanase was reported (White, A., Withers, S.G., Gilkes, N.R., Rose, D.R., 1994. Crystal structure of the catalytic domain of the beta-1,4-glycanase Cex from Cellulomonas fimi. Biochemistry 33, 12546-12552). The domain was decomposed into 22 modules by centripetal profiles (Go, M., Nosaka, M., 1987. Protein architecture and the origin of introns. Cold Spring Harbor Symp. Quant. Biol. 52, 915-924; Noguti, T., Sakakibara, H., Go, M., 1993. Localization of hydrogen-bonds within modules in barnase. Proteins 16, 357-363). A module is a contiguous polypeptide segment of amino acid residues having a compact conformation within a globular domain. Collected 31 intron sites of the family F/10 xylanase genes from fungus were found to be correlated to module boundaries with considerable statistical force (p values <0.001). The relationship between the intron locations and protein structures provides supporting evidence for the ancient origin of introns, because such a relationship cannot be expected by random insertion of introns into eukaryotic genes, but it rather suggests pre-existence of introns in the ancestral genes of prokaryotes and eukaryotes. A phylogenetic tree of the fungal and bacterial xylanase sequences made two clusters; one includes both the bacterial and fungal genes, but the other consists of only fungal genes. The mixed cluster of bacterial genes without introns and the fungal genes with introns further supports the ancient origin of introns. Comparison of the conserved base sequences of introns indicates that sliding of a splice site occurred in Aspergillus kawachii gene by one base from the ancestral position. Substrate-binding sites of xylanase are localized on eight modules, and introns are found at both termini of six out of these functional modules. This result suggests that introns might play a functional role in shuffling the exons encoding the substrate-binding modules.

An investigation of the nature and function of module 10 in a family F/10 xylanase FXYN of Streptomyces olivaceoviridis E-86 by module shuffling with the Cex of Cellulomonas fimi and by site-directed mutagenesis

Although the amino acid homology in the catalytic domain of FXYN xylanase from Streptomyces olivaceoviridis E‐86 and Cex xylanase from Cellulomonas fimi is only 50%, an active chimeric enzyme was obtained by replacing module 10 in FXYN with module 10 from Cex. In the family F/10 xylanases, module 10 is an important region as it includes an acid/base catalyst and a substrate binding residue. In FXYN, module 10 consists of 15 amino acid residues, while in Cex it consists of 14 amino acid residues. The K m and k cat values of the chimeric xylanase FCF‐C10 for PNP‐xylobioside (PNP‐X2) were 10‐fold less than those for FXYN. CD spectral data indicated that the structure of the chimeric enzyme was similar to that of FXYN. Based on the comparison of the amino acid sequences of FXYN and Cex in module 10, we constructed four mutants of FXYN. When D133 or S135 of FXYN was deleted, the kinetic properties were not changed from those of FXYN. By deletion of both D133 and S135, the K m value for PNP‐X2 decreased from the 2.0 mM of FXYN to 0.6 mM and the k cat value decreased from the 20 s−1 of FXYN to 8.7 s−1. Insertion of Q140 into the doubly deleted mutant further reduced the K m value to 0.3 mM and the k cat value to 3.8 s−1. These values are close to those for the chimeric enzyme FCF‐C10. These results indicate that module 10 itself is able to accommodate changes in the sequence position of amino acids which are critical for enzyme function. Since changes of the spatial position of these amino acids would be expected to result in enzyme inactivation, module 10 must have some flexibility in its tertiary structure. The structure of module 10 itself also affects the substrate specificity of the enzyme.

Sequence analysis of the gliding protein Gli349 in Mycoplasma mobile

The motile mechanism of Mycoplasma mobile remains unknown but is believed to differ from any previously identified mechanism in bacteria. Gli349 of M. mobile is known to be responsible for both adhesion to glass surfaces and mobility. We therefore carried out sequence analyses of Gli349 and its homolog MYPU2110 from M. pulmonis to decipher their structures. We found that the motif “YxxxxxGF” appears 11 times in Gli349 and 16 times in MYPU2110. Further analysis of the sequences revealed that Gli349 contains 18 repeats of about 100 amino acid residues each, and MYPU2110 contains 22. No sequence homologous to any of the repeats was found in the NCBI RefSeq non-redundant sequence database, and no compatible fold structure was found among known protein structures, suggesting that the repeat found in Gli349 and MYPU2110 is novel and takes a new fold structure. Proteolysis of Gli349 using chymotrypsin revealed that cleavage positions were often located between the repeats, implying that regions connecting repeats are unstructured, flexible and exposed to the solvent. Assuming that each repeat folds into a structural domain, we constructed a model of Gli349 that fits well the shape and size of images obtained with electron microscopy.

Correlation between amino acid residues converted by RNA editing and functional residues in protein three-dimensional structures in plant organelles

BackgroundIn plant organelles, specific messenger RNAs (mRNAs) are subjected to conversion editing, a process that often converts the first or second nucleotide of a codon and hence the encoded amino acid. No systematic patterns in converted sites were found on mRNAs, and the converted sites rarely encoded residues located at the active sites of proteins. The role and origin of RNA editing in plant organelles remain to be elucidated.ResultsHere we study the relationship between amino acid residues encoded by edited codons and the structural characteristics of these residues within proteins, e.g., in protein-protein interfaces, elements of secondary structure, or protein structural cores. We find that the residues encoded by edited codons are significantly biased toward involvement in helices and protein structural cores. RNA editing can convert codons for hydrophilic to hydrophobic amino acids. Hence, only the edited form of an mRNA can be translated into a polypeptide with helix-preferring and core-forming residues at the appropriate positions, which is often required for a protein to form a functional three-dimensional (3D) structure.ConclusionWe have performed a novel analysis of the location of residues affected by RNA editing in proteins in plant organelles. This study documents that RNA editing sites are often found in positions important for 3D structure formation. Without RNA editing, protein folding will not occur properly, thus affecting gene expression. We suggest that RNA editing may have conferring evolutionary advantage by acting as a mechanism to reduce susceptibility to DNA damage by allowing the increase in GC content in DNA while maintaining RNA codons essential to encode residues required for protein folding and activity.