Correlation between amino acid residues converted by RNA editing and functional residues in protein three-dimensional structures in plant organelles
BMC Plant Biology volume 8, Article number: 79 (2008)
In 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.
Here 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.
We 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.
RNA editing is a process that inserts, deletes and converts nucleotides in RNA after transcription, distinct from RNA splicing and 3' processing [1, 2]. The insertion/deletion type of RNA editing was first discovered in protozoan kinetoplastid mRNAs , and the conversion type of RNA editing was first discovered in the mammalian mRNA encoding apolipoprotein B (apoB) [4, 5], followed by its discovery in the plant mitochondrial mRNA for coxII  and chloroplast mRNA for rpl2 . Since then, conversion editing has been mostly found in mRNAs transcribed from the mitochondrial and chloroplast genomes of land plants [8–12]. In the mitochondrion of Arabidopsis thaliana, at least 441 nucleotides in mRNAs are subjected to RNA editing . In the chloroplast of black pine, at least 26 sites are edited ; in the chloroplast of the hornwort Anthoceros formosae, 942 RNA editing sites have been identified . The majority of conversion editing events in plant organelles occur within protein-coding regions of mRNAs, and involve cytidine-to-uridine (C-U) and sometimes uridine-to-cytidine (U-C) conversion . Earlier analysis  of the location of RNA targets within transcripts did not detect any rules to explain why particular residues within a protein sequence were affected by codon changes while other residues were not altered.
RNA editing often increases the percent identity of the encoded amino acid sequence to the homologous sequences, implying an important role for RNA editing in the function of encoded proteins. In maize chloroplast rpl2, the AUG initiation codon is generated by conversion of ACG . In cytochrome c oxidase subunit II, encoded by the mitochondrial DNA of Zea mays, a codon for a copper ligand residue was converted from the codon for Thr to that for Met; Met is required at the site to bind a copper ion, which is in turn prerequisite for electron transfer, the biological function of cytochrome c oxidase . In wheat mitochondrion ORF240, equivalent to cytochrome c biosynthesis protein, RNA editing converts a codon for one of the heme-binding residues to encode an amino acid appropriate for the heme interaction . Unedited psbF mRNA of spinach chloroplast causes a photosystem II-deficient phenotype ; unedited petB mRNA of tobacco chloroplast causes a defect in heme attachment to cytochrome b6 ; and unedited acetyl-coA carboxylase carboxyl transferase β of pea is not functional . These examples are exceptional cases that demonstrate the functional importance of RNA editing; for the majority of RNA editing events in organelles, however, the functional importance has not been specifically elucidated.
Genome sequencing and structural genomics projects have produced massive amounts of data, including RNA editing sites, organelle genome sequences, and protein 3D structures. In this report, we combine these data and computationally investigate implications for the functional roles of RNA editing. We define protein functions through protein 3D structures, and find that residues converted by RNA editing have significant bias toward structurally important sites.
Conversion-type RNA editing in DNA sequence databases
In Genbank release 158, there are 365 genes in plant organelles that undergo conversion editing at 3,560 nucleotides within their protein-coding regions (Table 1). Among these, 1,219 RNA editing events are observed in the first nucleotide of a codon, 1,983 events in the second nucleotide, and 358 events in the third nucleotide. The number of edited nucleotides in a codon is not limited to one. Out of the 3,560 events, 129 pairs of RNA editing events are targeted to the same codon within a gene. The first and second nucleotides of a codon are edited in 86 cases, the second and third in 33 cases, and the first and third in 10 cases.
Classification of the 365 gene products by sequence identity results in 88 protein families (Additional files 1 and 2). There are 1,923 unique RNA editing events; most of these are C-U conversions (Table 2). In Table 2, there are 13 events involving other types of conversions; all of these events were observed in the mRNA encoding mitochondrial cytochrome b6 from Pfiesterra piscicida .
These data show that ~90% (= (1219 + 1983)/3560) of RNA editing events are observed on the first or the second nucleotide of a codon; this observation suggests that RNA editing events often change the identity of the encoded amino acid. Table 3 shows the patterns of conversion of amino acid residues in the 1,923 RNA editing events. The top five patterns are Ser-Leu conversion (333) followed by Pro-Leu (325), Ser-Phe (248), Pro-Ser (101) and Arg-Trp (83). These conversions mostly restore evolutionarily conserved amino acid residues found in the multiple sequence alignment of homologous proteins (Fig. 1A). The patterns of conversion shown in Tables 2 and 3 are similar to the ones reported in previous studies [16, 23, 24].
Protein 3D structures of edited transcripts
Out of the 88 protein families whose mRNAs undergo conversion editing, 52 families contain members for which 3D structure data have been deposited in Protein Data Bank ; hence, their 3D structures can be modeled (Table 1). Out of the 1,923 RNA editing sites described above, 755 sites encode residues present in these 3D structures.
Correlation between functional residues and RNA editing sites
We have assigned functional residues based on the 3D structures of the 52 proteins, and collated lists of both the functional residues and those residues encoded by edited codons. The correspondence between the edited sites and functional sites is summarized in Figure 2.
We compared the ligand-binding residues and the residues converted by RNA editing in RuBisCO, for which binding sites for a ligand analogue are known , and found that none of the residues at the ligand-binding sites are converted by RNA editing (Fig. 1A). Out of the 755 RNA editing sites, only seven sites were found in ligand-binding sites (Additional file 1). First, a Zn-ligand residue in acetyl-CoA carboxylase carboxyltransferase β is encoded by an edited codon. In acetyl-CoA carboxylase carboxyltransferase β encoded by the chloroplast accD gene of Adiantum capillus-veneris, a CGC codon encoding Arg73 is converted to UGC encoding Cys. The Cys is a part of a zinc finger motif, and the deletion of the motif in plastid accD abrogated enzymatic activity . Second, Cu-ligand residues in cytochrome c oxidases I and II are encoded by codons undergoing RNA editing. In the mitochondrial coxI gene of Picea abies, Larix sp., Megaceros sp. and Zamia sp., a CAU codon encoding His, is converted to UAU encoding Tyr; this Tyr residue bridges an electron transfer pathway from heme to Cu. Replacing Tyr to Phe relocated the Cu, suggesting a role for Tyr in positioning the Cu into the proper electron pathway [28, 29]. In many mitochondrial coxII genes, an ACG encoding Thr is edited to AUG encoding Met; the Met residue is a ligand for Cu. The Cu ion in coxII is also on an electron transfer pathway. When this residue is mutated to Thr, CoxII almost completely loses its ability to bind copper . Third, Fe-S ligand residue in NADH dehydrogenase subunit I is targeted by RNA editing. A CGC codon encoding Arg70 is converted to a UGC encoding Cys in ndhI mRNA of chloroplast A. capillus-veneris. A Fe-S cluster binding to the Cys residue is a part of an electron transfer pathway . Finally, Mg-ATP binding residue in F1 ATPase α is encoded by an edited codon. A CCU encoding Pro in the mitochondrial gene of Nicotiana tabacum is converted to a UCU encoding Ser. However, the residue encodes Pro in the sequence whose 3D structure was determined , hence the editing seems to have no effect on the binding of ATP.
Residues in secondary structures
Many residues whose codons are converted by RNA editing are found to be a part of secondary structures in proteins, especially helices (Additional file 1). Out of the 755 RNA editing sites, 456 (about 60%) sites are found to locate in an α- or 310-helix structure and 89 locate in a β-sheet structure. The skewed distribution of RNA editing sites on helix is statistically significant (p < 8.38 × 10-9). In the current data, there are 131 cases in which the residue converted by RNA editing is in a helix but is originally encoded as Pro in the DNA. In these cases, a non-edited amino acid sequence would have a kinked or truncated helix, which would cause a defect in the protein 3D structure (because Pro is a helix-breaker ) and presumably result in a less stable protein. Thus, editing repairs plausible defects of the proteins by restoring residues to ones that are appropriate for secondary structure formation.
Residues in subunit interfaces
All of the proteins to which we assigned 3D structures in the current study are multi-subunit complexes; stable subunit interactions are an important factor for these proteins to function. RuBisCO consists of an octamer of heterodimers (large and small subunits) with an overall barrel structure. Each large subunit has interfaces to two small subunits and two other large subunits. Residues in the interfaces are shown in Fig. 1A by black triangles. Out of 30 unique residues encoded by codons converted by RNA editing, 13 residues are in the interfaces, and most of them are converted from hydrophilic to hydrophobic residues (Fig. 1). Out of the 755 RNA editing sites in the 52 protein families, 325 sites are located in the interfaces of subunits (Additional file 1). Conversion of Ser to Leu, Ser to Phe and Arg to Trp are respectively found 64, 57 and 23 times in the subunit interfaces. These conversions all switch the physicochemical properties of residues from hydrophilic to hydrophobic. This type of conversion increases the hydrophobicity of interface residues, one of the important properties of protein-protein interfaces .
Protein structural core
A key step in constructing the 3D structure of a protein is the formation of the structural core. RuBisCO large subunit has 139 (~30%) residues in structural cores (Fig. 1A), and residues converted by RNA editing are biased toward core-forming residues. In Figs. 1B–D, residues converted by RNA editing are shown in white, and those white residues are clustered inside the core-forming residues in red. Out of the 30 unique RNA editing sites in RuBisCO, 18 residues (60%) are located in structural cores, and these sites are over-represented in protein structural cores, considering that only 30% of all residues are in structural cores. The residue of Anthoceros formosae RuBisCO corresponding to residue 339 of Chlamydomonas reinhardtii RuBisCO (3D known) is originally encoded as Pro, and RNA editing of the codon converts the encoded residue to Leu (Fig. 1A). The residue is a part of a structural core and is located on an α helix. If the residue were not converted by RNA editing, then Pro would disrupt the α helix and the structural core, altering the local volume of the protein even if the resulting peptide was able to fold (Fig. 1D).
Out of the 755 RNA editing sites in the 52 proteins, 648 sites of RNA editing resulted in residue conversion. The other 107 (= 755-648) RNA editing sites are either involved in a stop codon or RNA editing that did not convert amino acid residues. Of the 648 RNA editing sites where the residue is converted, 191 (~30%) are targeted to residues in a protein structural core. In the 52 proteins with 3D structures, there are 12,370 residues and 2,331 residues (about 19%) are included in protein structural cores. Thus, the probability of obtaining the observed distribution is less than 4.54 × 10-10. However, since residues in a protein structural core tend to be hydrophobic and RNA editing often converts a residue to a hydrophobic one (Table 3), it could be a natural consequence that the RNA editing sites tend to be located in a protein structural core. We therefore performed the significance test only on Leu and Phe residues of the 52 proteins on all the sequences, and found that the probability is 7.81 × 10-5, still significant enough to support the hypothesis that the RNA editing sites are selectively located in codons for residues located within protein structural cores.
We have gathered the known RNA editing sites within protein-coding regions from the nucleotide sequence databases, and compared the location of RNA editing sites to functionally and structurally important sites within these proteins. In previous studies, residues converted by RNA editing were thought to be distributed without any rules in amino acid sequences , and only a few cases of RNA editing had been found to target active sites of proteins . In this study, we demonstrate that residues encoded by edited codons are biased toward helices, protein-protein interfaces and protein structural cores. The skew toward helices and protein structural cores are particularly statistically significant. The distributions suggest that the RNA editing sites are located on helices and protein core for some biological reason.
The recent study by Mulligan et al.  demonstrated with a sophisticated statistical method that edited codons are frequently grouped within nucleotide sequences, and these groups are separated by long gaps that contain no editing sites. Our analyses here may give structural context to their non-random distribution of edited sites. The groups of edited codons along the amino acid sequence might correspond to a cluster of residues within the protein structural core, and the long gaps between the groups of edited codons may correspond to the region that do not form a protein structural core, i.e., surfaces.
Functional impact of RNA editing through protein 3D structures
Mutation in a helix is known to have significant impact on protein stability. In T4 lysozyme, introduction of Pro into a helix resulted in significant kink structure and decreased stability by approximately 2.5 kcal/mol . Considering that the overall thermodynamic stability of protein is achieved by a free energy difference of ~5–15 kcal/mol, such a decrease in free energy difference between folded and unfolded states of the peptide could easily have an impact on the stability of a protein. An artificial back-mutation of Leu to Pro in the RNA editing site in a helix of cytochrome b6 was carried out; the authors of that study found that the mutated protein could not form a protein complex, apparently because cytochrome b6 itself did not fold appropriately . The direct translation of the original DNA sequences, therefore, seems to have repercussion in biological function through protein stability and/or quaternary structure formation; RNA editing restores the residues to make the protein stable and functional.
A mutation in a protein structural core can also have significant impact on protein stability. A mutation of Leu to Ser in a core, one of the typical conversions repaired by RNA editing, is equivalent to the loss of one methyl group and introduction of one hydroxyl group in a side chain of the residue; hence formation of hydrophilic cavity in a protein structural core is expected. The effect of cavity formation in a protein has been experimentally measured; it decreases stability ~3.3 kcal/mol . Alteration of Val to Ser in the core of ribonuclease T1 decreased stability of the protein by 4.7 kcal/mol . Two to three cavities in a protein structural core, expected in the unedited polypeptide, may therefore have a significant effect on protein stability. Sakaki et al.  expressed both edited and unedited acetyl-coA carboxylase carboxyl transferases β of pea, which had Leu and Ser, respectively, at residue 267, and measured the function of the protein in a complex with the α subunit. In addition to detecting no activity in the unedited complex, the authors found that the solubility of the unedited complex was low compared to the edited complex. Freezing and thawing of the eluate affected the unedited complex and resulted in an insoluble complex, whereas the disturbance did not affect the edited complex. Based on our 3D structure analysis, residue 267 is in a protein structural core; hence, this experiment directly showed that RNA editing in a structural core governs the stability of a protein. In the mitochondrion of Z. mays, the unedited mRNA encoding ribosomal protein S12 is translated into a polypeptide that cannot be incorporated into ribosomes . The product from the unedited mRNA had hydrophilic residues in place of hydrophobic residues required for protein core formation, and therefore would not form a structure stable enough to participate in supramolecule formation. Islas-Osuna et al.  suggested that several of the RNA editing sites on cytochrome b from grapevine were located either in the protein structural core or the interfaces with other subunits. Their analysis of one specific protein is consistent with our whole data analyses.
The conversion of amino acid residues within helices and structural cores by RNA editing restores residues that contribute to the formation of stable 3D structure. The unedited products are generally unstable or do not fold, and could disturb the protein networks in organelles that the proteins involve. A mechanism to post-transcriptionally switch amino acids within a protein core could be used to regulate functionality of the protein without regulating transcription of the gene for the protein. In addition, a mixture of edited and unedited products in a cell could yield multiple proteins with different stabilities, all derived from a single gene. If the ratio of edited products in the cell can be modified by regulation of RNA editing enzymes encoded in the nuclear genome [42, 43], then the efficiency of the biological functions involving the edited products can be controlled post-transcriptionally.
Implications for the evolution of RNA editing
The similarities between conversion editing observed in mitochondrion and chloroplast suggests that both mechanisms originated from a common source . Combining the discussions in the previous studies summarized below, we propose an evolutionary scenario for the origins of RNA editing in plant organelles that differs from the previously proposed ones [44, 45].
Three sets of observations may bear on discussions of the evolution of the mechanism of RNA editing in organelles. First, it has long been noted that the target C of RNA editing is often followed by a pyrimidine nucleotide, usually U . We counted the number of each nucleotide before and after the edited C on DNA in our dataset, and found that at the preceding site, T far exceeds the other types of nucleotides (Table 4). This bias of T suggests that the RNA editing recovers TT dinucleotide sequences that are underrepresented on DNA. Second, genome sequences of plant mitochondria are known to be GC rich, and there is a positive correlation between the number of RNA editing sites in mitochondrial genomes and their G+C content . RNA editing was then suggested to be a mechanism to compensate for the genetic drift from T to C , namely to counteract for GC pressure. The cause of GC pressure may vary, but one of the suggested physical causes is to avoid pyrimidine dimmer formation . Third, RNA editing in organelles has been almost exclusively found in land plants. Yoshinaga et al. proposed that RNA editing in chloroplasts had been acquired to effect the land adaptation of plants . When plants started to migrate to the land (about 0.4 billion years ago), the land was yet to be protected against ultraviolet (UV) light by the ozone layer [48, 49]. Whereas aquatic plants are protected by water from the hazardous effects of UV light (i.e., the formation of pyrimidine dimers on DNA) , land plants required a novel means to protect themselves.
The studies cited above conjectured that introduction of RNA editing in plant organelle was a positively selected countermeasure against pyrimidine dimer formation. In that case, RNA editing should have been introduced at almost all TT sites. After the reduction of UV light by ozone layer formation, RNA editing sites at the protein cores remained, while RNA editing sites that were not essential to protein structure formation, namely ones located on the surface, could disappear in a neutral manner. If direct translation of the encoded residue has little effect on protein 3D structure stability, then RNA editing on the site could disappear. The correlation between residues converted by RNA editing and their positions in protein structures could have emerged as an outcome of the process.
A different scenario for emergence of RNA editing can be considered. Mitochondria and chloroplasts are the primary sites for oxygen burning and oxygen production, respectively. Molecular oxygen can causes alterations in the chemical structure of the bases in DNA. In order to mitigate the effect of oxygen molecules on organellar DNAs, it would have been effective to both reduce the number of DNAs in the organellar genomes by gene transfer to the nucleus, and to compensate for mutations within organellar DNAs by RNA editing. In animal mitochondria, gene transfer to the nuclear genome has been observed. In plant organelles, both gene transfer and alteration of T to C on DNAs have been observed. T to C conversion on DNA would reduce the number of AT base pairs, which are more susceptible to oxygen molecules than GC base pairs .
We found statistically significant correlations between residues encoded by edited codons and the residues responsible for secondary structure and protein structural core formation. The correlation suggests that RNA editing affects protein functions indirectly by regulating protein stability as well as sometimes being essential for protein enzymatic activity. By repairing the nucleotide sequence of mRNA to encode a stable protein, the RNA editing machinery may regulate expression of protein functions in plant organelles.
Collection of data for conversion RNA editing
The descriptions of RNA editing sites in Genbank/EMBL/DDBJ [52–54] are not standardized, and there are a number of efforts (including the one here) to launch a database for RNA editing sites [55, 56]. In order to identify RNA editing sites, we performed a full-text search to find a string of characters that matches both "RNA" and "editing" in "/note" of "misc_feature" line of Genbank database release 158. We wrote a computer program to extract entries with RNA editing in protein-coding regions, and translated both edited and unedited mRNAs into amino acid sequences. We have encountered a significant number of errors in the Genbank annotation during this process, and corrected these annotations based on either literature or communication with depositors. Most of the errors took the form of discrepancies between the nucleotide position number described in misc_feature line and the RNA editing sites in deposited nucleotide sequences. The remediation was carried out on AB254134, AJ006146, AY820131, AY521591, BA000029, DQ645537, DQ984517, X69720, X92735 and Y17812. We could not correct all the errors we encountered, because we could not contact all the depositors. The entries with apparent errors were discarded. For RNA editing on rbcL transcripts, we copied the RNA editing site description into the following entries from a table found in reference ; D14882, D43696, L11055, L11056, and L13485. When we encountered a pair of protein sequences with identical amino acid sequences and identical RNA editing patterns in their mRNAs, we eliminated one of the entries from our dataset.
Protein 3D structures of RNA edited transcripts
We performed a homology search of amino acid sequences predicted from edited mRNAs against amino acid sequences of the protein with known 3D structures in PDB  using BLAST . When the sequence identity was 25% or more, we used the 3D structures in PDB for assigning structural properties of the products encoded by the edited and unedited mRNAs.
Functional residues on proteins
We gathered four types of protein functions based on the literature and protein 3D structures. (1) Ligand-binding residues: Most of the protein 3D structures were determined with their ligands and cofactors. Residues that bind those ligands were determined based on measuring solvent accessibilities of a residue with and without the ligand. When the difference between solvent accessibility of a residue calculated with and without the ligand was non-vanishing, then the residue was assigned as a ligand-binding residue. Solvent accessibility was calculated using a modified method of Sharke and Rupley  with water radius of 1.4Å. (2) Protein-protein interfaces: Most of the protein products of edited mRNAs are components of supramolecules, and their 3D structures were determined in protein complex form. Residues that interact with other subunits were determined by measuring the solvent accessibility of a residue with and without other subunits. (3) Secondary structure: We assigned secondary structures using DSSP . (4) Protein structural core: A structural core was determined using the following procedure; i) Calculate solvent accessibility, find residues with zero accessibility, and calculate all of the carbon atom distances between these residues; if the distance is no more than 4.0Å, then the pair of residues are parts of the protein structural core. ii) Calculate solvent accessibility, find residues with accessibility more than zero but no more than 0.05, calculate the distance between a carbon atom in the residue and a carbon atom in i); if the distances are 4.0Å or less, then the residue is a part of the protein structural core.
Keegan LP, Gallo A, O'Connell MA: The many roles of an RNA editor. Nat Rev Genet. 2001, 2 (11): 869-878.
Gott JM, Emeson RB: Functions and mechanisms of RNA editing. Annu Rev Genet. 2000, 34: 499-531.
Benne R, Burg van den J, Brakenhoff JP, Sloof P, van Boom JH, Tromp MC: Major transcript of the frameshifted coxII gene from trypanosoma mitochondria contains four nucleotides that are not encoded in the DNA. Cell. 1986, 46: 819-826.
Powell LM, Wallis SC, Pease RJ, Edwards YH, Knott TJ, Scott J: A novel form of tissue-specific RNA processing produces apolipoprotein-B48 in intesinte. Cell. 1987, 50: 831-840.
Chen SH, Habib G, Yang CY, Gu ZW, Lee BR, Weng SA, Silberman SR, Cai SJ, Deslypere JP, Rosseneu M, Gotto AM, Li WH, Chan L: Apolipoprotein B-48 is the product of a messenger RNA with an organ-specific in-frame stop codon. Science. 1987, 238: 363-366.
Covello PS, Gray MW: RNA editing in plant mitochondria. Nature. 1989, 341: 662-666.
Hoch B, Maier RM, Appel K, Igloi GL, Kössel H: Editing of a chloroplast mRNA by creation of an initiation codon. Nature. 1991, 353: 178-180.
Hiesel R, Combettes B, Brennicke A: Evidence for RNA editing in mitochondria of all major groups of land plants except the Bryophyta. Proc Natl Acad Sci USA. 1994, 91: 629-633.
Lin S, Zhang H, Spencer DF, Norman JE, Gray MW: Widespread and extensive editing of mitochondrial mRNAs in dinoflagellates. J Mol Biol. 2002, 320: 727-739.
Freyer R, Kiefer-Meyer M-C, Kössel H: Occurrence of plastid RNA editing in all major lineages of land plants. Proc Natl Acad Sci USA. 1997, 94: 6285-6290.
Yoshinaga K, Iinuma H, Masuzawa T, Ueda K: Extensive RNA editing of U to C in addition to C to U substitution in the rbcL transcripts of hornwort chloroplasts and the origin of RNA editing in green plants. Nucleic Acids Res. 1996, 24: 1008-1014.
Kugita M, Yamamoto Y, Fujikawa T, Matsumoto T, Yoshinaga K: RNA editing in hornwort chloroplasts makes more than half the genes functional. Nucleic Acids Res. 2003, 31: 2417-2423.
Giege P, Brennicke A: RNA editing in Arabidopsis mitochondria effects 441 C to U changes in ORFs. Proc Natl Acad Sci USA. 1999, 96: 15324-15329.
Wakasugi T, Hirose T, Horihata M, Tsudzuki T, Kössel H, Sugiura M: Creation of a novel protein-coding region at the RNA level in black pine chloroplasts: the pattern of RNA editing in the gymnosperm chloroplast is different from that in angiosperms. Proc Natl Acad Sci USA. 1996, 93: 8766-8770.
Gray MW, Covello PS: RNA editing in plant mitochondria and chloroplasts. FASEB J. 1993, 7: 64-71.
Bock R: Sense from nonsense: How the genetic information of chloroplasts is altered by RNA editing. Biochimie. 2000, 82: 549-557.
Covello PS, Gray MW: RNA sequence and the nature of the CuA-binding site in cytochrome c oxidase. FEBS lett. 1990, 268: 5-7.
Bonnard G, Grienenberger JM: A gene proposed to encode a transmembrane domain of an ABC transporter is expressed in wheat mitochondria. Mol Gen Genet. 1995, 246: 91-99.
Bock R, Kössel H, Maliga P: Introduction of a heterologous editing site into the tobacco plastid genome: the lack of RNA editing leads to a mutant phenotype. EMBO J. 1994, 13: 4623-4628.
Zito F, Kuras R, Choquet Y, Kössel H, Wollman FA: Mutations of cytochrome b6 in Chlamydomonas reinhardtii disclose the functional significance for a proline to leucine conversion by petB editing in maize and tobacco. Plant Mol Biol. 1997, 33: 79-86.
Sasaki Y, Kozaki A, Ohmori A, Iguchi H, Nagano Y: Chloroplast RNA editing required for functional acetyl-CoA carboxylase in plants. J Biol Chem. 2001, 276: 3937-3940.
Zhang H, Lin S: Mitochondrial cytochrome b mRNA editing in Dinoflagellates: Possible ecological and evolutionary associations?. J Eukaryot Microbiol. 2005, 52 (6): 538-545.
Maier RM, Zeltz P, Kössel H, Bonnard G, Gualberto JM, Grienenberger JM: RNA editing in plant mitochondria and chloroplasts. Plant Mol Biol. 1996, 32: 343-365.
Brennicke A, Marchfelder A, Binder S: RNA editing. FEMS Microbiol Rev. 1999, 23: 297-316.
Berman HM, Henrick K, Nakamura H: Announcing the worldwide Protein Data Bank. Nat Struct Biol. 2003, 10 (12): 980-
Karkehabadi S, Taylor TC, Spreitzer RJ, Andersson I: Altered intersubunit interactions in crystal structures of catalytically compromised ribulose-1,5-bisphosphate carboxylase/oxygenase. Biochemistry. 2005, 44: 113-120.
Kozaki A, Mayumi K, Sasaki Y: Thiol-disulfide exchange between nuclear-encoded and chloroplast-encoded subunits of pea acetyl-CoA carboxylase. J Biol Chem. 2001, 276: 39919-39925.
Tsukihara T, Aoyama H, Yamashita E, Tomizaki T, Yamaguchi H, Shinzawa-Itoh K, Nakashima R, Yaono R, Yoshikawa S: Structures of metal sites of oxidized bovine heart cytochrome c oxidase at 2.8Å. Science. 1995, 269: 1069-1074.
Trumpower BL, Gennis RB: Energy transduction by cytochrome complexes in mitochondrial and bacterial respiration: the enzymology of coupling electron transfer reactions to transmembrane proton translocation. Annu Rev Biochem. 1994, 63: 675-716.
Kelly M, Lappalainen P, Talbo G, Haltia T, Oost van der J, Saraste M: Two cysteines, two histidines, and one methionine are ligands of a binuclear purple copper center. J Biol Chem. 1993, 268: 16781-16787.
Sazanov LA, Burrows PA, Nixon PJ: The plastid ndh genes code for an NADH-specific dehydrogenase: isolation of a complex I analogue from pea thylakoid membranes. Proc Natl Acad Sci USA. 1998, 95: 1319-1324.
Bowler MW, Montgomery MG, Leslie AGW, Walker JE: How azide inhibits ATP hydrolysis by the F-ATPases. Proc Natl Acad Sci USA. 2006, 103: 8646-8649.
Chou PY, Fasman GD: Prediction of the secondary structure of proteins from their amino acid sequence. Adv Enzymol Relat Areas Mol Biol. 1978, 47: 45-148.
Wodak SJ, Janin J: Structural basis of macromolecular recognition. Adv Protein Chem. 2002, 61: 9-73.
Mulligan RM, Chang KLC, Chou CC: Computational analysis of RNA editing sites in plant mitochondrial genomes reveals similar information content and a sporadic distribution of editing sites. Mol Biol Evol. 2007, 24: 1971-1981.
Blaber M, Zhang XJ, Matthews BW: Structural basis of amino acid α helix propensity. Science. 1993, 260: 1637-1640.
Creighton TE: Protein folding. Biochem J. 1990, 270: 1-16.
Loladze VV, Ermolenko DN, Makhatadze GI: Thermodynamic consequences of burial of polar and non-polar amino acid residues in the protein interior. J Mol Biol. 2002, 320: 343-357.
Vos SD, Backmann J, Prevost M, Steyaert J, Loris R: Hydrophobic core manipulations in ribonuclease T1. Biochemistry. 2001, 40: 10140-10149.
Phreaner CG, WIlliams MA, Mulligan RM: Incomplete editing of rps12 transcripts results in the synthesis of polymorphic polypeptides in plant mitochondria. Plant Cell. 1996, 8: 107-117.
Islas-Osuna MA, Silva-Moreno B, Caceres-Carrizosa N, Garcia-Robles JM, Sotelo-Mundo RR, Yepiz-Plascencia GM: Editing of the grapevine mitochondrial cytochrome b mRNA and molecular modeling of the protein. Biochimie. 2006, 88: 431-435.
Zeltz P, Hess WR, Neckermann K, Borner T, Kossel H: Editing of the chloroplast rpoB transcript is independent of chloroplast translation and shows different patterns in barley and maize. EMBO J. 1993, 12: 4291-4296.
Halter CP, Peeters NM, Hanson MR: RNA editing in ribosome-less plastids of iojap maize. Curr Genet. 2004, 45: 331-337.
Tillich M, Lehwark P, Morton BR, Maier UG: The evolution of chloroplast RNA editing. Mol Biol Evol. 2006, 23: 1912-1921.
Covello PS, Gray MW: On the evolution of RNA editing. Trends Genetics. 1993, 9: 265-268.
Malek O, Lättig K, Hiesel R, Brennicke A, Knoop V: RNA editing in bryophytes and a molecular phylogeny of land plants. EMBO J. 1996, 15: 1403-1411.
Singer CE, Ames BN: Sunlight ultraviolet and bacterial DNA base ratios. Science. 1970, 170: 822-826.
Takahashi A, Ohnishi T: The significance of the study about the biological effects of solar ultraviolet radiation using the exposed facility on the international space station. Biol Sci Space. 2004, 18: 255-260.
Berkner LV, Marshall LC: History of major atmospheric components. Proc Natl Acad Sci USA. 1965, 53: 1215-1226. [http://www.pnas.org/content/53/6/1215.full.pdf+html]
Häder DP, Sinha RP: Solar ultraviolet radiation-induced DNA damage in aquatic organisms: potential environmental impact. Mutat Res. 2005, 571 (1-2): 221-233.
Lindahl T: Instability and decay of the primary structure of DNA. Nature. 1993, 362: 709-715.
Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL: GenBank. Nucl Acids Res. 2007, 35: D21-D25.
Kulikova T, Akhtar R, Aldebert P, Althorpe N, Andersson M, Baldwin A, Bates K, Bhattacharyya S, Bower L, Browne P, Castro M, Cochrane G, Duggan K, Eberhardt R, Faruque N, Hoad G, Kanz C, Lee C, Leinonen R, Lin Q, Lombard V, Lopez R, Lorenc D, McWilliam H, Mukherjee G, Nardone F, Pastor MPG, Plaister S, Sobhany S, Stoehr P, Vaughan R, Wu D, Zhu W, Apweiler R: EMBL Nucleotide Sequence Database in 2006. Nucl Acids Res. 2007, 35: D16-D20.
Sugawara H, Abe T, Gojobori T, Tateno Y: DDBJ working on evaluation and classification of bacterial genes in INSDC. Nucl Acids Res. 2007, 35: D13-D15.
He T, Du P, Li Y: DbRES: a web-oriented database for annotated RNA editing sites. Nucl Acids Res. 2007, 35: D141-144.
Picardi E, Regina TMR, Brennicke A, Quagliariello C: REDIdb: the RNA editing database. Nucl Acids Res. 2007, 35: D173-D177.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucl Acids Res. 1997, 25: 3389-3402.
Shrake A, Rupley JA: Environment and exposure to solvent of protein atoms. Lysozyme and insulin. J Mol Biol. 1973, 79: 351-371.
Kabsch W, Sander C: Dictionary of protein secondary structure: pattern recognition of hydrogen-bonded and geometrical features. Biopolymers. 1983, 22: 2577-2637.
KY thanks Ms. Kazuko Kaji and Ms. Atsuko Doi for collecting descriptions of RNA editing sites from the literature, and Mr. Kazuhiro Kobayashi and Mr. Masataka Oikawa for assisting data analyses. KY and MG were supported by Grant-in-Aid for Scientific Research (B) from the Japan Society for the Promotion of Science (JSPS) KAKENHI (18370061).
KY built the database, calculated the correlation and drafted the manuscript. MG started and supervised the study. All authors read and approved the final manuscript.
Electronic supplementary material
Additional file 1: supplementary table 1. A list of residues converted by RNA editing in 52 proteins with known 3D structures. The table also contains secondary structures, solvent accessibility and the species name from which the gene was derived. (XLS 174 KB)
Additional file 2: supplementary table 2. A list of residues converted by RNA editing in 37 proteins without known 3D structures. (XLS 194 KB)
About this article
Cite this article
Yura, K., Go, M. Correlation between amino acid residues converted by RNA editing and functional residues in protein three-dimensional structures in plant organelles. BMC Plant Biol 8, 79 (2008). https://doi.org/10.1186/1471-2229-8-79