BACKGROUND
rna editing is a post-transcriptional modification of mrna  <cit> , which may result in the synthesis of proteins that are not directly encoded in the genome.
there are two major types of rna editing in mammals, both of which occur via deamination of a base, either cytidine  or adenosine .
inosine is read by the ribosome  as guanosine, and thus a → i modifications at the mrna level translate into an a → g changes at the genetic code level.
in this work we focus exclusively on a-to-i rna editing, which is catalyzed by enzymes from the adar  family.
adars are double-stranded rna  binding proteins, and thus dsrna is a prerequisite for a-to-i editing  <cit> .
rna editing is a fine-tuning mechanism, capable of changing only a few nucleotides.
both edited and unedited variants of the same transcript may be present in the cell.
a-to-i editing is known to be vital in vertebrates, and important for normal life in invertebrates.
in drosophila, knocking out adar activity causes the flies to exhibit defects in locomotion and mating and to suffer tremors  <cit> .
adar knockout c. elegans worms exhibit chemotaxis defects  <cit> .
in mice, knocking out adar <dig> causes embryonic death and defects in erythropoiesis  <cit> .
adar <dig> -/- mice die shortly after birth and are increasingly seizure prone after postnatal day  <dig>  <cit> .
the lethal phenotype is accounted for by a single editing site resulting in a single amino acid substitution in the glur-b gene.
in addition, alteration of a → i editing has been ascribed to several pathological conditions  <cit> , mainly to neuro-psychiatric conditions such as amyotrophic lateral sclerosis   <cit> , epilepsy  <cit> , major depression disorder  <cit> , and glioblastoma multiforme  <cit> .
reduced a-to-i editing levels have been linked to cancer in various tissues, most strongly to brain tumors.
a correlation between the reduction of adar <dig> and the tumor aggressiveness was observed, and overexpression of adar <dig> and adar <dig> resulted in decreased proliferation rate of the glioblastoma multiforme cell-lines  <cit> .
isolating inosine-containing transcripts from c. elegans and human brain, it has been noticed that most a-to-i editing occurs in non coding regions  <cit> .
genome-wide bioinformatic searches for a-to-i editing sites have enabled the identification of abundant a-to-i editing in the transcriptome of several vertebrates  <cit> .
it was found that editing occurs mainly within repetitive elements.
these repetitive elements are likely to base-pair with a neighboring similar element and form the dsrna structure which is the target of the adar enzymes.
in particular, virtually all a-to-i editing events in human occur specifically within alu repeats.
the alus are a particular set of primate-specific retrotransposons, approximately  <dig> nucleotides in length.
the alus are the most abundant of all transposable elements in primates, making up more than 10% of the human genome, with some  <dig>  million copies.
recent studies  <cit>  have demonstrated that the frequency of a-to-i editing in human is much higher than in mouse, rat, chicken and fly.
this has to do with the abundance and low diversity of the alu elements as compared to similar elements in other genomes  <cit> : since alu is so common in the human genome, there is a high probability that an alu and a counterpart, oppositely oriented alu, exist nearby and are transcribed together.
when the rna transcript folds, these two alus form a helix, thus becoming a target for the dsrna binding adars.
the physiological significance of a-to-i editing within non-coding repetitive elements is still elusive.
several possible mechanisms have been suggested through which editing of a non-coding repetitive element might affect the fate of a transcript: editing may result in insertion or elimination of a splice site, and may theoretically lead to the alteration of transcriptional start and stop codons  <cit> .
hyperedited inosine-containing rnas might be cleaved at specific sites  <cit> .
in addition, inosine containing mrnas were also shown to be retained in the nucleus, suggesting an additional regulatory role for a-to-i editing  <cit> .
however, the validity and scope of this last mechanism has been debated recently  <cit> .
finally, while the molecular significance is yet unclear, editing within alu repeats was shown to be altered in cancerous tissues  <cit> .
a-to-i editing is characterized by a puzzling specificity and selectivity in the adenosines which are edited.
in some substrates, e.g.
the ampa receptor glur-b subunit in mice  <cit>  and the e <dig> sites within an alu repeat in the narf gene  <cit> , rna editing is extremely efficient, editing 100% of transcripts at a specific adenosine.
in others, such as most of the sites in alu repeats, a seemingly random editing pattern is observed, where many adenosines are targeted, with varying editing efficiency.
however, careful analysis reveals that editing in alu repeats is also highly reproducible: the variability among healthy individuals in editing level at a given site within a specific alu repeat is much lower than the site-to-site differences.
sequence preferences for adars have been previously documented.
c and t are overrepresented at the nucleotide 5' to the editing site, while g is underrepresented.
at the nucleotide 3' to the site, g is significantly overrepresented  <cit> .
these motifs are too weak, however, to fully characterize a-to-i editing.
therefore, the question still stands: what controls the editing level at each given site?
adars bind to the rna via double-stranded rna binding motifs.
thus, dsrna is a necessity for a-to-i editing.
indeed, it has been shown for the highly selective r/g editing site within the hairpin of the glutamate receptor subunits mrnas, that the identities of bases in the helical region are evolutionarily conserved, while the bases in the nonhelical part of the hairpin covary so as to maintain their non-helical structure  <cit> .
this distinctive feature demonstrates the importance of the secondary structure to the phenomenon of rna editing.
the internal structure of the dsrna is expected to control the editing efficiency  <cit> .
for example, it has been shown experimentally that internal loops may effectively be equivalent to helix termini in terms of editing efficiency  <cit> .
thus, internal loops along dsrna, if large enough, may act as delimiters separating a large dsrna into many small helices.
since adars deaminate fewer a's in shorter helices, their existence  might be a means to increase the specificity of editing.
it is thus plausible that more features of the secondary structure of an rna molecule play an important role in determining the specificity of adenosine deamination of an adar substrate.
in this paper we will characterize the properties of a-to-i editing sites in terms of their secondary structure properties, their sequence properties, and their thermodynamic properties.
we describe the building of a database of mfold <cit>  foldings used to query these properties, and then display and discuss the results of those queries.
RESULTS
structural analysis
we first look at the editing frequency for each substructure type .
we compare a "test set" of a-to-i editing sites, which we denote by e <dig>  and a control set of sites not known to be edited, denoted by e <dig>  the e <dig> and e <dig> sets are defined with precision in the methods section.
interestingly, while the existence of a helix is well known to be a prerequisite for editing, the overall frequency of e <dig> is actually more than two fold lower in helices  than in interior loops .
as the overwhelming majority of e <dig> sites reside in helices and interior loops, we focus henceforth on these two substructures only.
for clarity, we emphasize here that by "interior loop" we mean only the unpaired nucleotides that form the loop's constituent strands.
size is the number of nucleotides in the structure for bulge, hairpin, interior, junction and strand, and the length of the helix for helices.
the table lists average sizes with 95% confidence intervals.
"strand" here means an rna strand which is not part of a loop .
furthermore, we study the effect of the location of the specific nucleotide within its respective substructure.
we define cepos as the distance of the site  from the closest edge of the substructure it is in .
figs.
<dig> and  <dig> present the frequency of e <dig> sites as a function of cepos.
for helices, one observes a general trend of enhancement of editing as a site lies deeper in the helix.
for interior loops, however, there is dramatic depletion of e <dig> for cepos >  <dig>  in fact, it should be noted that 91% of edited sites in interior loops lie at the very edge of the loop, i.e.
cepos =  <dig>  most of these are in fact a single mismatch within an almost perfect helix .
such mismatches were already implicated as preferred targets of adars, as previous in-vitro data as well as bioinformatic work indicate that ac mismatches are more favorable substrates than a-t pairs  <cit> .
however, it is worthwhile noticing that our analysis shows this trend to persist even for longer interior loops: interior loop strands of length up to five nucleotides are more likely to be edited than the average site in a helix .
for these cepos =  <dig> sites, there is a significant  effect to the direction of the nearest neighboring helix: a-to-i editing frequency is  <dig>  for sites with a helix only in the upstream site,  <dig>  for sites with a helix only in the downstream site, and  <dig>  for sites with helices on both sides.
the above results hold when controlling for the total length of the substructure: we compared e <dig> and e <dig> sites for helices of a given length, and for loops of a given size.
the resulting trend was the same: for e <dig> sites in helices cepos is larger than for e <dig> sites, whereas in interior loops the connection is reversed.
other location variables tested, such as the position relative to the middle of the substructure, or relative to the 5' end, did not result in noticeable results.
sequence analysis
we start with the nucleotide opposite of the editing site.
for helices, it is clear what this means: the "opposite" nucleotide of a site is the nucleotide that pairs with that site .
we expand this idea, however, to sites at the edges of interior loops : for these sites on the most 5'  nucleotide of the loop-strand, the opposite nucleotide is the most 3'  site of the other strand in the loop.
if the site is the only one on its strand, and the opposite strand has more than one nucleotide, the opposite nucleotide is undefined.
we shall refer to the opposite nucleotide as opnuc for short.
there is a very strong enrichment for sites with c on the opposite site: we looked at the frequency of e <dig> for sites with a given opnuc, and obtained a frequency of  <dig> % for c, whereas for a the frequency was  <dig> % and for g,  <dig> %.
this is consistent with  the data presented in  <cit> .
next we look at the statistics of the nucleotides upstream and downstream of the a-to-i editing sites.
in order to avoid biases due to the underlying nucleotide statistics in alu repeats we do not look at the raw distribution of nucleotides but rather at the enrichment factor, i.e.
how much is the editing frequency increased  when the neighboring site is any specific nucleotide.
the enrichment factors are presented in figs  <dig>   <dig> and  <dig> for the two immediate neighbor nucleotides separately, as well as for the joint variable composed of both upstream and downstream neighbors.
overall, the profiles found are similar to those seen in previous large-scale studies of editing  <cit> : t is most preferred upstream and is not preferred downstream, while g is most preferred downstream and least preferred upstream .
however, we do find a significant  difference between the profiles for helices and loops.
for example, the preference for an upstream t is stronger in helices, whereas the preference for a downstream g is stronger in interior loops.
we also calculated the enrichment factors for the joint variable composed of the site's upstream neighbor, downstream neighbor, and opnuc.
the results are displayed in table  <dig>
in addition, we searched for enrichment in the extended neighborhood of the editing sites, looking at  <dig> neighboring nucleotides at both sides of the site .
almost all neighbors show a significantly different nucleotide distribution around edited sites, see tables  <dig>  and  <dig> .
the most significant differences  are observed for neighbors up <dig>  up <dig>  up <dig> and dn <dig> in helices and up <dig>  dn <dig>  up <dig> and up <dig> in interior loops.
we note that while almost all  <dig> neighbors tested show statistically significant difference, it is hard to tell whether these differences are due to adars preference or rather stem from editing hot spots within the alu.
we also present the enrichment factors for seven positions surrounding the editing sites which were reported to show preferences to specific nucleotides when surrounding adar <dig> editing sites  <cit> .
as seen in table  <dig>  the patterns observed here for alu editing are somewhat different: for example, locations dn <dig> and dn <dig> seem to favor g in contrast to the opposite trend reported in  <cit>  for adar <dig> sites.
the differences might be due to the much larger sample we study here.
additionally, it is also possible that editing sites in the coding region, mostly having a functional role, have different characteristics than the ones in alu repeats.
however, these differences could also result from differences between the profiles of adar <dig> and adar <dig>  while the sample of editing sites studied in  <cit>  is biased towards adar <dig> targets, the sample studied here, coming from a wide range of tissues, represents a different mix of the two enzymes, with larger weight of adar <dig>  moreover, the different splice-variants of the adars possibly have varying editing efficiencies and site preferences.
the mix of these variants occurring in our in-vivo sample, could also lead to slight variations in the preferences observed as compared to results of in-vitro studies.
the neighbor index gives the location of the nucleotide relative to the edited  a, where negative values correspond to upstream nucleotides.
the total number of e <dig> sites here was  <dig> , and about  <dig>  e <dig> sites .
the neighbor index gives the location of the nucleotide relative to the edited  a, where negative values correspond to upstream nucleotides.
the total number of e <dig> sites here was  <dig> , and about  <dig>  e <dig> sites .
nucleotide distributions at certain locations around editing sites, reported to exhibit nucleotide biases  <cit> .
for each of the sites, we present the probability to have a given nucleotide n when the  <dig> location adenosine is edited, divided by the probability of that nucleotide regardless of whether the  <dig> adenosine is edited: p/p.
values >  <dig> indicate enrichment of n for edited sites, and <  <dig> indicate depletion.
upx = x nucleotides upstream of site  <dig>  dnx = x nucleotides downstream.
the probabilities are given for editing sites in helices and interior loops separately, but are very similar for both.
for comparison, we present the patterns reported in  <cit> .
thermodynamic calculations
finally, we study the effect of thermodynamic stability on editing efficiency.
for each genomic neighborhood, we look at the thermodynamic average over all the low free-energy structures.
the laws of statistical mechanics give us the probability that the rna is in a specific secondary structure n,
where t is the temperature in degrees kelvin, kb is boltzmann's constant and z is defined by the sum
where the label n runs over all available foldings, and gn are the respective free energies.
in practice, we only use those folds generated by mfold which are expected to be all folds relevant at human-body temperature.
we may now, for example, calculate the probability of some particular site to be in a helix,
where  is the indical function, defined by
similarly, one may calculate the probabilities for all other substructures.
let s denote the set of possible substructures,
we define a site's structural entropy to be
where  is the frequency of site i being in substructure of type s. this entropy is a measure of the thermodynamic volatility of the site's substructure: if a site is always in the same substructure , it will have zero structural entropy.
if, however, the site's substructure fluctuates, for example between a helical structure and a loop structure, it will have higher structural entropy.
the structural entropy of a site with equal probability to be in two difference substructures is ln =  <dig> .
the highest possible structural entropy is of a site which spends equal time in each of the possible substructures.
figs.
<dig> and  <dig> show the frequency of e <dig> as a function of the structural entropy, for sites whose lowest free-energy structure is helix or interior-loop separately.
interestingly, a-to-i editing is enriched for sites with low structural entropy, i.e.
having a well-defined low energy micro-structure.
a wobbling state, fluctuating between two or more possible structures is less well edited.
this holds regardless of the ground-state structure, but the effect is stronger for interior loops: sites with a well-defined interior loops structure are twice more likely to be edited compared with sites whose ground state structure is also an interior loop but having even 1% probability to be in other structures.
analysis of a large dataset of secondary structures of putatively edited alu repeats reveals that structural motifs are indeed important in determining the a-to-i editing efficiency of a given site.
most notably, we highlight the strong preference for editing of adenosines within short symmetrical internal loops.
moreover, the microstructure also has modest but noticeable effect on the cis-preferences of the adars.
long perfect dsrna duplexes are often considered to be the best target for editing by adars.
here we find that sites adjacent to the edge of helices  are even more efficiently edited.
averaged over our whole database, adenosines deep within  long  perfect helices are indeed edited more efficiently than the average adenosine in a helix - we find  <dig> such sites, with editing frequency  <dig> %, compared to only  <dig> % for the average helix site.
however, this is still lower than the average frequency for interior loops,  <dig> %.
moreover, focusing on single a-c mismatches within a helix  raises the frequency to 19%.
finally, choosing also the optimal neighbors, i.e.
t upstream and g downstream, raises the editing frequency as observed in our database to 37% !
we stress again that these frequencies should not be regarded as the true editing efficiency, but rather as a relative measure.
yet, one is able to conclude that the best way to engineer an efficient editing site is not to put it in a long perfect duplex, but rather to be in a single mismatch within a duplex.
interestingly, the 100% edited e <dig> site in the narf gene  <cit> , fits nicely with these engineering rules - it is a cepos =  <dig> site in a symmetric loop, with c opposite to it and t and g in the upstream and downstream sites, respectively.
however, the strand length there is  <dig> and not the optimal  <dig>  an editing site that fully adheres to the above "rules" is the amber/w one of the hepatitis delta virus antigenome   <cit> .
this site is critical for the virus to assemble viral particles and to be infectious  <cit> .
given the high adaptivity of viruses, it is not surprising to find that this site fits with all of the above: it is located in a single a-c mismatch within a helix , and has t and g as its immediate neighbors.
however, the q/r site in glurb does not fit to our observations.
it lies within a rather long  helix, with cepos =  <dig>  with c  upstream and g downstream  <cit> .
yet, this site is also 100% edited.
apparently, there is still much more to learn about the characteristics of editing by adars, beyond the information presented in the present study.
one possible explanation is that this site in known to be edited only by adar <dig>  <cit> .
the two editing enzymes adar <dig> and adar <dig> are known to have overlapping, but distinct, preferences  <cit> .
however, our approach does not allow us to distinguish between them.
it was recently shown that editing of mouse b <dig> and b <dig> sines is mediated by both enzymes  <cit> .
some sites within these repeats are adar <dig> specific, some are adar <dig> specific and some are edited by both.
it is not yet clear which enzyme is the main one in terms of alu editing in human.
since our database is based on mrna sequences from a wide range of tissues, it is possible that it characterized mainly the profile of the widely-expressed adar <dig> rather than that of adar <dig> which is expressed mainly in the brain.
it is thus likely that some of the preferences identified in this work characterize adar <dig> and are therefore not present in the glurb adar2-specific site.
the discrepancies between nucleotide distributions around the editing sites reported above and those reported by  <cit>  for adar <dig> sites might also attest for differences between the adar <dig> profile and the one characterizing our dataset, probably a mix of the two enzymes, with larger weight of adar <dig>
in an attempt to estimate the differences between the two enzymes, we compared  <dig> editing sites supported by  <dig> brain mrnas, where both adar <dig> and adar <dig> are present, and  <dig> sites residing in  <dig> non-brain mrnas, presumably edited mostly by adar <dig> .
the patterns observed were similar but not identical.
for example,  <dig> of the  <dig> brain sites residing in a helix  had a g in the dn <dig> position, compared to  <dig> out of  <dig>  in non-brain sites in a helix .
however, differences were not statistically significant upon bonferroni-correcting for multiple testing.
thus, a larger and better dataset  is required in order to study the small tissue differences between the preferences of the two adar enzymes.
CONCLUSIONS
using a dataset of  <dig>  editing sites within alu repeats, we analyzed the editing preferences.
we found that the micro-structure a site resides in affects its editing frequency.
in addition, the sequence motifs characterizing editing sites vary with the micro-structure.
we have also shown that structural entropy and thermodynamic stability play a role in determining editing efficiency.
we find that the probability of a nucleotide fluctuating between a number of possible structures to be edited is lower than the weighted average of the probabilities for each possible structure alone.
this provides a hint as to the rate of the a-to-i editing process compared to the relaxation time scales controlling the transition between the possible folds.
taken together, the results presented here could be of help in revealing the complex nature of the adars editing profile.

 rna editing is a post-transcriptional modification of mrna, which may result in the synthesis of proteins that are not directly encoded in the genome. there are two major types of rna editing in mammals, both of which occur via deamination of a base, either cytidine or adenosine. inosine is read by the ribosome as guanosine, and thus a → i modifications at the mrna level translate into an a → g changes at the genetic code level. in this work we focus exclusively on a-to-i rna editing, which is catalyzed by enzymes from the adar family. adars are double-stranded rna binding proteins, and thus dsrna is a prerequisite for a-to-i editing. rna editing is a fine-tuning mechanism, capable of changing only a few nucleotides. both edited and unedited variants of the same transcript may be present in the cell. a-to-i editing is known to be vital in vertebrates, and important for normal life in invertebrates. in drosophila, knocking out adar activity causes the flies to exhibit defects in locomotion and mating and to suffer tremors. adar knockout c. elegans worms exhibit chemotaxis defects. in mice, knocking out adar causes embryonic death and defects in erythropoiesis. adar -/- mice die shortly after birth and are increasingly seizure prone after postnatal day. the lethal phenotype is accounted for by a single editing site resulting in a single amino acid substitution in the glur-b gene. in addition, alteration of a → i editing has been ascribed to several pathological conditions, mainly to neuro-psychiatric conditions such as amyotrophic lateral sclerosis, epilepsy, major depression disorder, and glioblastoma multiforme. reduced a-to-i editing levels have been linked to cancer in various tissues, most strongly to brain tumors. a correlation between the reduction of adar and the tumor aggressiveness was observed, and overexpression of adar and adar resulted in decreased proliferation rate of the glioblastoma multiforme cell-lines. isolating inosine-containing transcripts from c. elegans and human brain, it has been noticed that most a-to-i editing occurs in non coding regions. genome-wide bioinformatic searches for a-to-i editing sites have enabled the identification of abundant a-to-i editing in the transcriptome of several vertebrates. it was found that editing occurs mainly within repetitive elements. these repetitive elements are likely to base-pair with a neighboring similar element and form the dsrna structure which is the target of the adar enzymes. in particular, virtually all a-to-i editing events in human occur specifically within alu repeats. the alus are a particular set of primate-specific retrotransposons, approximately nucleotides in length. the alus are the most abundant of all transposable elements in primates, making up more than 10% of the human genome, with some million copies. recent studies have demonstrated that the frequency of a-to-i editing in human is much higher than in mouse, rat, chicken and fly. this has to do with the abundance and low diversity of the alu elements as compared to similar elements in other genomes: since alu is so common in the human genome, there is a high probability that an alu and a counterpart, oppositely oriented alu, exist nearby and are transcribed together. when the rna transcript folds, these two alus form a helix, thus becoming a target for the dsrna binding adars. the physiological significance of a-to-i editing within non-coding repetitive elements is still elusive. several possible mechanisms have been suggested through which editing of a non-coding repetitive element might affect the fate of a transcript: editing may result in insertion or elimination of a splice site, and may theoretically lead to the alteration of transcriptional start and stop codons. hyperedited inosine-containing rnas might be cleaved at specific sites. in addition, inosine containing mrnas were also shown to be retained in the nucleus, suggesting an additional regulatory role for a-to-i editing. however, the validity and scope of this last mechanism has been debated recently. finally, while the molecular significance is yet unclear, editing within alu repeats was shown to be altered in cancerous tissues. a-to-i editing is characterized by a puzzling specificity and selectivity in the adenosines which are edited. in some substrates, e.g. the ampa receptor glur-b subunit in mice and the e sites within an alu repeat in the narf gene, rna editing is extremely efficient, editing 100% of transcripts at a specific adenosine. in others, such as most of the sites in alu repeats, a seemingly random editing pattern is observed, where many adenosines are targeted, with varying editing efficiency. however, careful analysis reveals that editing in alu repeats is also highly reproducible: the variability among healthy individuals in editing level at a given site within a specific alu repeat is much lower than the site-to-site differences. sequence preferences for adars have been previously documented. c and t are overrepresented at the nucleotide 5' to the editing site, while g is underrepresented. at the nucleotide 3' to the site, g is significantly overrepresented. these motifs are too weak, however, to fully characterize a-to-i editing. therefore, the question still stands: what controls the editing level at each given site? adars bind to the rna via double-stranded rna binding motifs. thus, dsrna is a necessity for a-to-i editing. indeed, it has been shown for the highly selective r/g editing site within the hairpin of the glutamate receptor subunits mrnas, that the identities of bases in the helical region are evolutionarily conserved, while the bases in the nonhelical part of the hairpin covary so as to maintain their non-helical structure. this distinctive feature demonstrates the importance of the secondary structure to the phenomenon of rna editing. the internal structure of the dsrna is expected to control the editing efficiency. for example, it has been shown experimentally that internal loops may effectively be equivalent to helix termini in terms of editing efficiency. thus, internal loops along dsrna, if large enough, may act as delimiters separating a large dsrna into many small helices. since adars deaminate fewer a's in shorter helices, their existence might be a means to increase the specificity of editing. it is thus plausible that more features of the secondary structure of an rna molecule play an important role in determining the specificity of adenosine deamination of an adar substrate. in this paper we will characterize the properties of a-to-i editing sites in terms of their secondary structure properties, their sequence properties, and their thermodynamic properties. we describe the building of a database of mfold foldings used to query these properties, and then display and discuss the results of those queries.  structural analysis we first look at the editing frequency for each substructure type. we compare a "test set" of a-to-i editing sites, which we denote by e and a control set of sites not known to be edited, denoted by e the e and e sets are defined with precision in the methods section. interestingly, while the existence of a helix is well known to be a prerequisite for editing, the overall frequency of e is actually more than two fold lower in helices than in interior loops. as the overwhelming majority of e sites reside in helices and interior loops, we focus henceforth on these two substructures only. for clarity, we emphasize here that by "interior loop" we mean only the unpaired nucleotides that form the loop's constituent strands. size is the number of nucleotides in the structure for bulge, hairpin, interior, junction and strand, and the length of the helix for helices. the table lists average sizes with 95% confidence intervals. "strand" here means an rna strand which is not part of a loop. furthermore, we study the effect of the location of the specific nucleotide within its respective substructure. we define cepos as the distance of the site from the closest edge of the substructure it is in. figs. and present the frequency of e sites as a function of cepos. for helices, one observes a general trend of enhancement of editing as a site lies deeper in the helix. for interior loops, however, there is dramatic depletion of e for cepos > in fact, it should be noted that 91% of edited sites in interior loops lie at the very edge of the loop, i.e. cepos = most of these are in fact a single mismatch within an almost perfect helix. such mismatches were already implicated as preferred targets of adars, as previous in-vitro data as well as bioinformatic work indicate that ac mismatches are more favorable substrates than a-t pairs. however, it is worthwhile noticing that our analysis shows this trend to persist even for longer interior loops: interior loop strands of length up to five nucleotides are more likely to be edited than the average site in a helix. for these cepos = sites, there is a significant effect to the direction of the nearest neighboring helix: a-to-i editing frequency is for sites with a helix only in the upstream site, for sites with a helix only in the downstream site, and for sites with helices on both sides. the above results hold when controlling for the total length of the substructure: we compared e and e sites for helices of a given length, and for loops of a given size. the resulting trend was the same: for e sites in helices cepos is larger than for e sites, whereas in interior loops the connection is reversed. other location variables tested, such as the position relative to the middle of the substructure, or relative to the 5' end, did not result in noticeable results. sequence analysis we start with the nucleotide opposite of the editing site. for helices, it is clear what this means: the "opposite" nucleotide of a site is the nucleotide that pairs with that site. we expand this idea, however, to sites at the edges of interior loops: for these sites on the most 5' nucleotide of the loop-strand, the opposite nucleotide is the most 3' site of the other strand in the loop. if the site is the only one on its strand, and the opposite strand has more than one nucleotide, the opposite nucleotide is undefined. we shall refer to the opposite nucleotide as opnuc for short. there is a very strong enrichment for sites with c on the opposite site: we looked at the frequency of e for sites with a given opnuc, and obtained a frequency of % for c, whereas for a the frequency was % and for g, %. this is consistent with the data presented in. next we look at the statistics of the nucleotides upstream and downstream of the a-to-i editing sites. in order to avoid biases due to the underlying nucleotide statistics in alu repeats we do not look at the raw distribution of nucleotides but rather at the enrichment factor, i.e. how much is the editing frequency increased when the neighboring site is any specific nucleotide. the enrichment factors are presented in figs and for the two immediate neighbor nucleotides separately, as well as for the joint variable composed of both upstream and downstream neighbors. overall, the profiles found are similar to those seen in previous large-scale studies of editing: t is most preferred upstream and is not preferred downstream, while g is most preferred downstream and least preferred upstream. however, we do find a significant difference between the profiles for helices and loops. for example, the preference for an upstream t is stronger in helices, whereas the preference for a downstream g is stronger in interior loops. we also calculated the enrichment factors for the joint variable composed of the site's upstream neighbor, downstream neighbor, and opnuc. the results are displayed in table in addition, we searched for enrichment in the extended neighborhood of the editing sites, looking at neighboring nucleotides at both sides of the site. almost all neighbors show a significantly different nucleotide distribution around edited sites, see tables and. the most significant differences are observed for neighbors up up up and dn in helices and up dn up and up in interior loops. we note that while almost all neighbors tested show statistically significant difference, it is hard to tell whether these differences are due to adars preference or rather stem from editing hot spots within the alu. we also present the enrichment factors for seven positions surrounding the editing sites which were reported to show preferences to specific nucleotides when surrounding adar editing sites. as seen in table the patterns observed here for alu editing are somewhat different: for example, locations dn and dn seem to favor g in contrast to the opposite trend reported in for adar sites. the differences might be due to the much larger sample we study here. additionally, it is also possible that editing sites in the coding region, mostly having a functional role, have different characteristics than the ones in alu repeats. however, these differences could also result from differences between the profiles of adar and adar while the sample of editing sites studied in is biased towards adar targets, the sample studied here, coming from a wide range of tissues, represents a different mix of the two enzymes, with larger weight of adar moreover, the different splice-variants of the adars possibly have varying editing efficiencies and site preferences. the mix of these variants occurring in our in-vivo sample, could also lead to slight variations in the preferences observed as compared to results of in-vitro studies. the neighbor index gives the location of the nucleotide relative to the edited a, where negative values correspond to upstream nucleotides. the total number of e sites here was, and about e sites. the neighbor index gives the location of the nucleotide relative to the edited a, where negative values correspond to upstream nucleotides. the total number of e sites here was, and about e sites. nucleotide distributions at certain locations around editing sites, reported to exhibit nucleotide biases. for each of the sites, we present the probability to have a given nucleotide n when the location adenosine is edited, divided by the probability of that nucleotide regardless of whether the adenosine is edited: p/p. values > indicate enrichment of n for edited sites, and < indicate depletion. upx = x nucleotides upstream of site dnx = x nucleotides downstream. the probabilities are given for editing sites in helices and interior loops separately, but are very similar for both. for comparison, we present the patterns reported in. thermodynamic calculations finally, we study the effect of thermodynamic stability on editing efficiency. for each genomic neighborhood, we look at the thermodynamic average over all the low free-energy structures. the laws of statistical mechanics give us the probability that the rna is in a specific secondary structure n, where t is the temperature in degrees kelvin, kb is boltzmann's constant and z is defined by the sum where the label n runs over all available foldings, and gn are the respective free energies. in practice, we only use those folds generated by mfold which are expected to be all folds relevant at human-body temperature. we may now, for example, calculate the probability of some particular site to be in a helix, where is the indical function, defined by similarly, one may calculate the probabilities for all other substructures. let s denote the set of possible substructures, we define a site's structural entropy to be where is the frequency of site i being in substructure of type s. this entropy is a measure of the thermodynamic volatility of the site's substructure: if a site is always in the same substructure, it will have zero structural entropy. if, however, the site's substructure fluctuates, for example between a helical structure and a loop structure, it will have higher structural entropy. the structural entropy of a site with equal probability to be in two difference substructures is ln =. the highest possible structural entropy is of a site which spends equal time in each of the possible substructures. figs. and show the frequency of e as a function of the structural entropy, for sites whose lowest free-energy structure is helix or interior-loop separately. interestingly, a-to-i editing is enriched for sites with low structural entropy, i.e. having a well-defined low energy micro-structure. a wobbling state, fluctuating between two or more possible structures is less well edited. this holds regardless of the ground-state structure, but the effect is stronger for interior loops: sites with a well-defined interior loops structure are twice more likely to be edited compared with sites whose ground state structure is also an interior loop but having even 1% probability to be in other structures. analysis of a large dataset of secondary structures of putatively edited alu repeats reveals that structural motifs are indeed important in determining the a-to-i editing efficiency of a given site. most notably, we highlight the strong preference for editing of adenosines within short symmetrical internal loops. moreover, the microstructure also has modest but noticeable effect on the cis-preferences of the adars. long perfect dsrna duplexes are often considered to be the best target for editing by adars. here we find that sites adjacent to the edge of helices are even more efficiently edited. averaged over our whole database, adenosines deep within long perfect helices are indeed edited more efficiently than the average adenosine in a helix - we find such sites, with editing frequency %, compared to only % for the average helix site. however, this is still lower than the average frequency for interior loops, %. moreover, focusing on single a-c mismatches within a helix raises the frequency to 19%. finally, choosing also the optimal neighbors, i.e. t upstream and g downstream, raises the editing frequency as observed in our database to 37%! we stress again that these frequencies should not be regarded as the true editing efficiency, but rather as a relative measure. yet, one is able to conclude that the best way to engineer an efficient editing site is not to put it in a long perfect duplex, but rather to be in a single mismatch within a duplex. interestingly, the 100% edited e site in the narf gene, fits nicely with these engineering rules - it is a cepos = site in a symmetric loop, with c opposite to it and t and g in the upstream and downstream sites, respectively. however, the strand length there is and not the optimal an editing site that fully adheres to the above "rules" is the amber/w one of the hepatitis delta virus antigenome. this site is critical for the virus to assemble viral particles and to be infectious. given the high adaptivity of viruses, it is not surprising to find that this site fits with all of the above: it is located in a single a-c mismatch within a helix, and has t and g as its immediate neighbors. however, the q/r site in glurb does not fit to our observations. it lies within a rather long helix, with cepos = with c upstream and g downstream. yet, this site is also 100% edited. apparently, there is still much more to learn about the characteristics of editing by adars, beyond the information presented in the present study. one possible explanation is that this site in known to be edited only by adar. the two editing enzymes adar and adar are known to have overlapping, but distinct, preferences. however, our approach does not allow us to distinguish between them. it was recently shown that editing of mouse b and b sines is mediated by both enzymes. some sites within these repeats are adar specific, some are adar specific and some are edited by both. it is not yet clear which enzyme is the main one in terms of alu editing in human. since our database is based on mrna sequences from a wide range of tissues, it is possible that it characterized mainly the profile of the widely-expressed adar rather than that of adar which is expressed mainly in the brain. it is thus likely that some of the preferences identified in this work characterize adar and are therefore not present in the glurb adar2-specific site. the discrepancies between nucleotide distributions around the editing sites reported above and those reported by for adar sites might also attest for differences between the adar profile and the one characterizing our dataset, probably a mix of the two enzymes, with larger weight of adar in an attempt to estimate the differences between the two enzymes, we compared editing sites supported by brain mrnas, where both adar and adar are present, and sites residing in non-brain mrnas, presumably edited mostly by adar. the patterns observed were similar but not identical. for example, of the brain sites residing in a helix had a g in the dn position, compared to out of in non-brain sites in a helix. however, differences were not statistically significant upon bonferroni-correcting for multiple testing. thus, a larger and better dataset is required in order to study the small tissue differences between the preferences of the two adar enzymes.  using a dataset of editing sites within alu repeats, we analyzed the editing preferences. we found that the micro-structure a site resides in affects its editing frequency. in addition, the sequence motifs characterizing editing sites vary with the micro-structure. we have also shown that structural entropy and thermodynamic stability play a role in determining editing efficiency. we find that the probability of a nucleotide fluctuating between a number of possible structures to be edited is lower than the weighted average of the probabilities for each possible structure alone. this provides a hint as to the rate of the a-to-i editing process compared to the relaxation time scales controlling the transition between the possible folds. taken together, the results presented here could be of help in revealing the complex nature of the adars editing profile. 
finally, we show that thermodynamic stability of the site is important for its editing. surprisingly, we find that interior loops, and especially the nucleotides at their edges, are more likely to be edited than helices. several studies have identified weak sequence motifs characterizing the editing sites, but these alone do not account for the large diversity observed. in addition, the sequence motifs characterizing editing sites vary with the micro-structure. analysis of a large dataset of editing events reveals more information on sequence and structural motifs characterizing the a-to-i editing process alu repeats in the human transcriptome undergo massive adenosine to inosine rna editing. we focus on structural aspects, studying the double-stranded rna structure of the alu repeats, and show the editing frequency of a given site to depend strongly on the micro-structure it resides in.
 alu repeats in the human transcriptome undergo massive adenosine to inosine rna editing. this process is selective, as editing efficiency varies greatly among different adenosines. several studies have identified weak sequence motifs characterizing the editing sites, but these alone do not account for the large diversity observed.  here we build a dataset of editing sites and use it to characterize editing preferences. we focus on structural aspects, studying the double-stranded rna structure of the alu repeats, and show the editing frequency of a given site to depend strongly on the micro-structure it resides in. surprisingly, we find that interior loops, and especially the nucleotides at their edges, are more likely to be edited than helices. in addition, the sequence motifs characterizing editing sites vary with the micro-structure. finally, we show that thermodynamic stability of the site is important for its editing.  analysis of a large dataset of editing events reveals more information on sequence and structural motifs characterizing the a-to-i editing process 