Nucleic Acids Research

This Article Abstract Print PDF (328K) Supplementary Material Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in ISI Web of Science Similar articles in PubMed Alert me to new issues of the journal Add to My Personal Archive Download to citation manager Search for citing articles in: ISI Web of Science (2284) Request Permissions Commercial Re-use Guidelines for Open Access NAR Content Google Scholar Articles by Zuker, M. Search for Related Content PubMed PubMed Citation Articles by Zuker, M. Social Bookmarking          What's this?

Nucleic Acids Research, 2003, Vol. 31, No. 13 3406-3415
© 2003 Oxford University Press

Mfold web server for nucleic acid folding and hybridization prediction

Michael Zuker^*

Department of Mathematical Sciences, Rensselaer Polytechnic Institute, Troy, NY 12180, USA

*Tel: +1 5182766902; Fax: +1 5182764824; Email: zukerm{at}rpi.edu

Received February 14, 2003; Revised and Accepted April 7, 2003

	ABSTRACT

TOP ABSTRACT INTRODUCTION SERVER CONTENT AND ORGANIZATION EQUIPMENT AND ORGANIZATION FUTURE DIRECTIONS CITING THE MFOLD WEB... SUPPLEMENTARY MATERIAL REFERENCES

 
The abbreviated name, ‘mfold web server’, describes a number of closely related software applications available on the World Wide Web (WWW) for the prediction of the secondary structure of single stranded nucleic acids. The objective of this web server is to provide easy access to RNA and DNA folding and hybridization software to the scientific community at large. By making use of universally available web GUIs (Graphical User Interfaces), the server circumvents the problem of portability of this software. Detailed output, in the form of structure plots with or without reliability information, single strand frequency plots and ‘energy dot plots’, are available for the folding of single sequences. A variety of ‘bulk’ servers give less information, but in a shorter time and for up to hundreds of sequences at once. The portal for the mfold web server is http://www.bioinfo.rpi.edu/applications/mfold. This URL will be referred to as ‘MFOLDROOT’.

	INTRODUCTION

TOP ABSTRACT INTRODUCTION SERVER CONTENT AND ORGANIZATION EQUIPMENT AND ORGANIZATION FUTURE DIRECTIONS CITING THE MFOLD WEB... SUPPLEMENTARY MATERIAL REFERENCES

 
The concept of RNA secondary structure began with the work of Doty and Fresco (1,2). The prediction of RNA secondary structure (folding) by energy minimization using nearest neighbor energy parameters began with Tinoco and colleagues (3–6) and also with Delisi and Crothers (7). Efficient algorithms for RNA secondary structure prediction using dynamic programming methods borrowed from sequence alignment were developed independently by a number of people (8–13).

My own early RNA folding programs (12,14) computed a single minimum energy folding of an RNA sequence. They were popular in the 1980s and a modified version was incorporated into the UWGCG (University of Wisconsin Genetics Computer Group) suite of programs (15); the acronym was shortened to GCG when this group left the University of Wisconsin to form a private company. Initially, they used free energy parameters that had been summarized by Salser (16). After 1986, free energies from the Turner group (Doug H. Turner, Department of Chemistry, University of Rochester, Rochester, NY) were used (17).

The ‘mfold’ software for RNA folding was developed in the late 1980s (18). The ‘m’ simply refers to ‘multiple’. The core algorithm predicts a minimum free energy, G, as well as minimum free energies for foldings that must contain any particular base pair. Any base pair, r_i-r_j, between the ith nucleotide and the jth nucleotide that is contained in a folding no more than G from the minimum, is plotted in a triangular plot called the ‘energy dot plot’. The base pair r_i-r_j is plotted in row i and column j of this matrix. The free energy increment, G, is chosen a priori by the user, who selects a ‘percent suboptimality’, P. From this, G is computed to be P/100 |G|. Base pairs within this free energy increment are chosen either automatically, or else by the user, and foldings that contain the chosen base pair are computed. They have minimum free energy conditional on containing the chosen base pair. The description and use of the mfold package has appeared in a number of articles (19–22). The closely related ‘RNAstructure’ program has also been described (23,24).

The Turner group has published numerous articles over the years that detail the development of the RNA folding parameters. A subset of these articles are what I would call ‘major works’ that summarize the current state of the art. Version 1 of the mfold package used free energies that were described by Freier et al. (17). Versions 2.1 to 2.3 used the parameters from Walter et al. (25), although the incorporation of coaxial stacking parameters into the minimization algorithm has not been accomplished. The current version 3 software uses free energy data from Mathews et al. (26).

DNA folding prediction with the mfold software began in 1996, when DNA specific parameters were added to the mfold package through a collaboration with the SantaLucia group (John SantaLuica Jr., Department of Chemistry, Wayne State University, Detroit, MI). These data have been described by SantaLucia (27). The DNA stacking (27), single mismatch (28–32) and dangling end (33) parameters have been measured in the SantaLucia laboratories. The remaining terminal stacking and loop parameters were estimated by SantaLucia and have been incorporated into the mfold package by personal communication from John SantaLucia. They remain unpublished. In 1999, corrections for [Na⁺], [Mg⁺⁺] were incorporated into the mfold package (34).

The mfold web server was first created at Washington University School of Medicine during the fall of 1995. DNA folding parameters were added in the spring of 1996. From 1995 until the fall of 2000, the server ran on SGI workstations or multiprocessor servers (Silicon Graphics, Inc. 1600 Amphitheatre Pkwy., Mountain View, CA 94043), as well as on a dual processor ‘Intel/Solaris’ platform (Pentium 2 processors by Intel Corp., Solaris operating system by Sun Microsystems). The server was moved to Rensselaer Polytechnic Institute (RPI) in October 2000. It ran on a dual processor ‘Intel/Linux’ platform (Linux operating system as developed and marketed by Red Hat, Inc., 1801 Varsity Drive, Raleigh, NC 27606). Since July 2002, the mfold web server has been running on a cluster of 36 dual processor ‘Intel/Linux’ workstations that were donated to the joint RPI-Wadsworth Bioinformatics Center by IBM (IBM Research, P.O. Box 218, Yorktown Heights, NY 10598). This equipment was awarded as an SUR (Shared University Research) grant to RPI and the Wadsworth Center (PI: M. Zuker).

The use of the mfold web server has grown steadily since its inception. The ‘quikfold’ server was added while I was still at Washington University. The remaining servers that will be described have all been added within that past 2 years. The servers have been used extensively by researchers in universities, medical schools, non-profit organizations, US government and military laboratories and by companies all over the world. In addition, the servers are being used in the teaching of computational biology methods.

	SERVER CONTENT AND ORGANIZATION

TOP ABSTRACT INTRODUCTION SERVER CONTENT AND ORGANIZATION EQUIPMENT AND ORGANIZATION FUTURE DIRECTIONS CITING THE MFOLD WEB... SUPPLEMENTARY MATERIAL REFERENCES

 
The mfold web server comprises a number of separate applications that predict nucleic acid folding, hybridization and melting temperatures (T_ms).

The basic ‘mfold’ server: input
The original applications on the mfold web server deal with folding a single RNA or DNA sequence per submission (job). The submission forms for RNA and DNA are separate for historical reasons only. RNA folding came first. These applications may be reached by following the ‘RNA Folding’ or ‘DNA Folding’ hyperlinks from the main portal or entrance page. The URL for the portal will likely remain stable, while those for separate applications will change. The default RNA folding form currently uses the latest version 3.0 free energies (26). These are recommended for most RNA folding. However, there is a link from this page to what is called the ‘RNA mfold version 2.3 server’. This server offers RNA folding using the older, version 2.3 energy parameters (25). Why use older and less accurate parameters? The reason is that we have enthalpies for these older parameters. As with the free energy parameters, the enthalpies were measured at 37°C. However, they are assumed to be constant within the range of temperatures that might occur in vivo or in the laboratory. This enables the server to extrapolate free energies to other temperatures and to fold at these temperatures.

Sequence name.
A sequence name may be typed or pasted (entered) within the ‘Enter a name for your sequence:’ text field. Long names are truncated to 40 characters. Any ASCII characters may be used, including those with octal values greater than 200 Octal (O). The ‘dangerous’ characters, ", <, > and ‘, are converted to ’, <<(253 O), >>(273 O) and ' (264 O), respectively. The character, \, is eliminated. If no name is entered, then the sequence name becomes the ‘Job ID’, which is of the form yyMmmdd-hh-mm-ss, where yy is year, Mmm is month, dd is day, hh is hour, mm is minute and ss is second. For example, a job that comes into the server at 8:23:06 pm local time on 9 February 2003, will be assigned a Job ID of 03Feb09-20-23-06. If two jobs come in during the same second, the second one has the letter ‘a’ appended to it, and so on for more than two jobs during the same second.

Sequence.
A sequence must be entered into the sequence text area box. All characters except for ‘A–Z’ and ‘a–z’ are removed. Lower case characters are converted to upper case. For RNA folding, ‘T’ or ‘t’ are converted to ‘U’, while ‘U’ or ‘u’ are converted to ‘T’ in DNA folding. If, for example, you enter

then the resulting sequence will be

rather than

The letter ‘N’ should be used for an unspecified base. It is not allowed to pair. The letters ‘B’, ‘D’, ‘H’ and ‘V’ denote ‘A’, ‘C’, ‘G’ and ‘U/T’ respectively. These nucleotides may pair only if their 3' neighbor is unpaired. The purpose of this convention is to denote nuclease cleavage of the phosphodiesther bond linking the indicated nucleotide and its 3' neighbor. It is used to constrain folding when nuclease digestion data specify susceptible bonds. In addition, the letters ‘W’, ‘X’, ‘Y’ and ‘Z’ also refer to ‘A’, ‘C’, ‘G’ and ‘U/T’, respectively. These nucleotides, if they pair, should do so only at the end of a helix. Thus, the mfold web server does not support the IUPAC (International Union of Pure and Applied Chemistry) ambiguous DNA character convention (35) shown in Table 1.

View this table: [in this window] [in a new window]   Table 1. IUPAC codes for incompletely specified bases in nucleic acid sequences

 
Constraints.
The text area in the constraints box allow for the optional incorporation of folding constraints. Each constraint consists of a single line in the box that must conform to a rigid format. When constraints are used, a hyperlink labeled ‘Explanation of sequence annotation’ appears on the primary results page.

The various types of constraints are itemized below. Multiple constraints of any form are allowed in any order.

• Force a specific base pair or helix to form. The command
  F i j k
  
• will force the formation of the helix (single base pair if k=1)
  
  The triple (i, j, k) refers to k consecutive base pairs, where r_i·r_j is the exterior closing base pair. If any of these base pairs cannot exist, then an error will be generated and the job will fail. The usual result is an output page that declares ‘Job aborted! No Structure!’. In the text area of the output page, successfully forced base pairs are underlined with (‘and’), while those that are forced to pair but cannot are underlined with ‘!’. Note that isolated base pairs are not allowed by the folding code. That is, even if r_i·r_j is a valid base pair, it will not be allowed to form unless either r_i+1·r_j-1 or r_i-1·r_j+1 is a valid base pair.
  
• Prohibit a specific base pair or helix from forming. The command
  P i j k
  will prohibit every single base pair of the form r_i+h·r_j-h, 0≤h≤k, from occurring. Base pairs that are prohibited that could not form in any case will be underlined with ‘!’ in the text area of the primary output page, but the folding will proceed without error.
  
• Force a string of consecutive bases to pair. The command
  F i 0 k
  (the second to last character is ‘zero’) will force nucleotides r_i, r_i+1, r_i+2,..., r_i+k-1 to pair. This is a single base when k=1. Forcing too many bases to pair or forcing a base labeled ‘N’ to pair will generate a fatal error.
  
• Prohibit a string of consecutive bases from pairing. The command
  P i 0 k
  (the second to last character is ‘zero’) will prevent nucleotides r_i, r_i+1, r_i+2,..., r_i+k-1 from pairing. This is a single base when k=1. Forcing too many bases to be single stranded can generate a fatal error.
  
• Prohibit a string of consecutive bases from pairing with another string. The command
  P i j k l
  
• will prevent the nucleotides r_i, r_i+1, r_i+2,..., r_j from pairing with nucleotides r_k, r_k+1, r_k+2 ,..., r_l (i≤j and k≤l). Note that if i=k and j=l, then the constraint is equivalent to forbidding all base pairs within the segment r_i ,..., r_j.
  

Other folding parameters

• RNA and DNA sequences may be linear or circular. The default is ‘linear’, but ‘circular’ may be chosen using an option button on the folding form. It is worth noting that the algorithm to fold circular nucleic acids is simpler than that for folding linear ones.
  
• The folding temperature is fixed at 37°C for RNA folding using version 3.0 energy rules. For RNA folding with the version 2.3 parameters, or for DNA folding, any integral temperature between 0 and 100°C may be chosen.
  
• Ionic conditions may be altered for DNA folding only. For RNA, the ionic conditions are fixed at [Na⁺]=1 M and [Mg⁺⁺]=0 M. For folding, these are equivalent to physiological conditions. The following constraints apply:
  [Na⁺]≥0.01 M,
  [Mg⁺⁺]≤0.1 M, and
  [Na⁺]≤0.3 M if [Mg⁺⁺]>0 M.
  For the purposes of folding, Na⁺ may be considered equivalent to Li⁺, K⁺ and NH₄⁺, while Mg⁺⁺ is equivalent to Ca⁺⁺.
  
• The percent suboptimality, P, controls the free energy increment, G for displaying base pairs in the energy dot plot and for computing suboptimal foldings. Base pairs that can occur in foldings with free energies ≤G+G will be plotted, and only foldings with free energies ≤G+ G will be computed. Normally, G=P/100|G|, but it is rounded up to 1 kcal/mol or down to 12 kcal/mol if outside this range.
  
• The upper bound on the number of foldings is an absolute limit. The number of foldings computed may be less than this quantity. Normally, it should be less that this number because the number of computed foldings is more appropriately controlled through a proper choice of P and the window parameter.
  
• The window parameter, W, controls the number of foldings that are computed. It may be thought of as a distance parameter. The distance between 2 bp, r_i·r_j and r_i'·r_j' may be defined as max{|i-i'|, |j-j'|}. If k-1 foldings have already been predicted by mfold, the kth folding must have at least W base pairs that are at least a distance W from any of the base pairs in the first k-1 foldings. A new folding is not added to the output list unless this criterion is fulfilled As W increases, the number of predicted foldings decreases. If W is not specified, mfold selects a value by default based on sequence length.
  
• If the maximum distance between paired bases parameter, M, is specified, then any base pair, r_i·r_j, in a folding of a linear molecule must satisfy j-i≤M. In a circular molecule, the condition becomes min{j-i, N +i-j}≤M, where N is the sequence length.
  

Immediate versus batch jobs.
Folding sequences containing up to 800 bases may be done while the user waits. This is the default. At this time, folding results for immediate jobs are erased 24–26 h after they are submitted. For sequences from 801 to 6000 bases in length, the batch option must be selected from the appropriate option button. The user should enter a valid email address in this case, although email addresses are always welcome since they identify users. For batch jobs, pressing the ‘Fold RNA’ or ‘Fold DNA’ submission button takes the user to a notification page that indicates a URL for the results and the email address to which a notification will be sent when the folding is completed. If the link to the results pages is followed immediately, then the target page will be incomplete. It will contain only the sequence and some other input information. When the folding is complete, this page must be ‘refreshed’ or ‘reloaded’ for the results to be seen. The URL for the results remains valid for 48–51 h, after which the results are erased. Folding sequences of up to 10 400 bases is available at special request. Folding 6000 bases currently takes about 1.5 h. The hyperlink labeled ‘View Folding Results’ will enable a user to view any results on the server, provided that they have not been erased and that the user views from the same computer that generated the foldings. That is, the user's IP address must be the same.

Output parameters

• The default, regular image resolution value, is 72x72 pixels per inch² for png and jpg images. The low, medium and high values are 60x60, 110x110 and 200x200, respectively. This parameter has no effect on PostScript output files, as well as the pdf files that can be derived (‘distilled’) from them. These files can be scaled up as much as desired without loss in quality. Jpg and png images will always be displayed at 72x72 pixels per square inch. Changing the resolution will change the size of the entire image. Image sizes are given in Table 2.
  
• Structure format plots are drawn showing individual bases or in an outline mode where bases do not appear. The outline option is suitable for images of very large folded molecules. The default choice of ‘Automatic’ will cause bases to be displayed in foldings of up to 800 nt and an outline to be drawn otherwise. Choosing the ‘Bases’ option for a large sequence might make sense if the user intends to magnify the structure to display a portion in which individual bases can be seen.
  
• The grid lines in the energy dot plot may be turned off.
  
• By default, bases will be numbered according to the length of the sequence. The default values are shown in Table 3. The user is free to alter this value, or to turn off base numbering altogether by choosing a frequency of zero.
  
• The user may alter the automatically chosen orientation of the folded molecular by choosing a rotation angle. Positive values correspond to counter-clockwise.
  
• Structure annotation has been described by Zuker and Jacobson (36). Bases in plotted structures may be annotated by ‘p-num’ values, which represent the number of ways that a base may pair in all foldings within G from the minimum energy. Low values indicate ‘well-defined’ bases. In particular, values of 0 or 1 indicate that a base is always single stranded or always paired to a unique partner, respectively. The number of times that a base is single stranded in the computed foldings is called its ‘ss-count’ number, and structure plots may also be annotated using these numbers. Color schemes are shown at MFOLDROOT/www-NAR03/doc/colors.html. Finally, for RNA folding only, the ‘high-light’ option allows the user to specify regions in the sequence to be highlighted. The selected bases are drawn in green, while the remaining bases are drawn in black. If the structure format is ‘Bases’, then the base characters will be drawn in color. If the structure format is ‘Outline’, then colored dots that represent bases will be superimposed on the ladder-like outline plot. If the default ‘Automatic’ option is used, then colored dots are used if the sequence size is >800. However, if the sequence size is ≤800, then both colored dots and bases are used. The bases appear in black or white inside an appropriately colored dot.
  

View this table: [in this window] [in a new window]   Table 2. Resolutions and images sizes of jpg and png images, the units are pixels per square inch and pixels, respectively

 

View this table: [in this window] [in a new window]   Table 3. Default numbering increments for secondary structures

 

Folding results
Folding results may be viewed by following links from an initial results page. Some input information is reproduced at the top of this page in a text area. An example of the top portion of an initial results page is shown in Figure 1.

View larger version (133K): [in this window] [in a new window]   Figure 1. This figure reproduces part of the initial results page for an RNA folding. Some input is reproduced in the text area. Various hyperlinks lead to the actual folding results.

 
The energy dot plot.
The energy dot plot is available in PostScript, png and jpg form. The ‘Text’ hyperlink leads to a plain text file that gives the basic dot plot information. The ‘istart’, ‘jstart’ and ‘length’ options correspond to the (i, j, k) numbers that define a helix, as defined previously. The ‘energy’ in an integer in units of tenth of a kcal/mol. Level 1 corresponds to helices in optimal foldings. Levels 2–4 correspond to the helices in increasingly suboptimal foldings. The levels from 1 to 4 correspond to the default coloring of the base pairs in the dot plot. The PostScript hyperlink leads to a static PostScript image of the energy dot plot. Letting G be the minimum free energy of a folding containing a base pair, the default color (rgb Hex) of a base pair in the energy dot plot is given in Table 4.

View this table: [in this window] [in a new window]   Table 4. The default colors used in the energy dot plot

 
The png and jpg hyperlinks each lead to a new page containing a png or jpg image of the energy dot plot and buttons that allow the user to interact with or redraw the dot plot. It is called the primary dot plot page. A small ‘Details’ window, 550x40 (pixels) opens in the top left corner of the screen. With the ‘Click on a base pair to identify’ radio button selected (the default), clicking (left click) on the image will cause the (i, j) value of the selected base pair as well as its G value to appear in the Details window. If a blank area of the plot is chosen, the value of G will be declared ‘Undefined’.

Selecting the ‘Click on the image and redraw the dot plot with options below’ button activates the text area and five option buttons below. Clicking on a point in the dot plot will cause a new image to be created. If the magnification factor is chosen to be >1, then the new image will be magnified about the selected point; otherwise the entire dot plot will be redrawn. The value of the ‘Energy increment in kcal/mole’ option will be an upper bound on G. For example, setting it to zero will cause only optimal base pairs to be plotted. The number of different colors used in the dot plot may be changed. If k is the value of the ‘Filter’ option, then only helices of size at least k will be plotted when the dot plot is redrawn. The ‘Image Width’ controls the resolution (size) of the redrawn image as explained above. If the PostScript ‘Output’ option is chosen, then a static PostScript image is created. If the png or jpg ‘Output’ option is chosen, then clicking on the plot will cause a new window to open containing the redrawn image and the same set of control buttons as before. This is called the secondary dot plot window. New png or jpg images created from this secondary window will be drawn within the secondary window. Similarly, new png or jpg images generated from the primary dot plot window will refresh the secondary window. A typical use would be to magnify about a region in the initial dot plot. Base pairs in the magnified dot plot can then be selected with precision, allowing the user to identify specific helices that may occur in optimal or suboptimal foldings. Closing the secondary window will not affect other windows. Closing the primary window will also close the Details window.

RNAML syntax.
Computed secondary structures are also available in the RNAML (37) format. This format has been proposed as a syntax for exchanging RNA structure information. We are currently offering output that is compatible with the DTD (Document Type Definition) file maintained by the Major group at the Université de Montréal. This file may be found at: http://www-lbit.iro.umontreal.ca/rnaml/current/rnaml.dtd.

Some extra files.
These are raw text files that contain information that might be useful to some users. Annotation by p-num value is available on the server, as described above. The raw values are also made available. The h-num link contains h-num values of each helix in the energy dot plot, sorted from lowest to highest value of h-num. The h-num value of a helix is the ‘well-definedness’ of a helix (38). Helices with lower h-num values tend to be more reliably predicted. The log file for a run is sometimes useful in finding the cause of a failed job.

Bulk downloading of results.
All computed foldings may be downloaded at once in a single zipped file (.zip) or a compressed tar file (.tar.Z). The structures may be taken as PostScript, png or jpg plot files, or as raw text files in different formats.

The text files will be useful to those who wish to create structure plots on their own. Some of the text formats have been described before (21,22). A full description may be found at MFOLDROOT/www-NAR03/doc/structure-format.html. Redrawn structures or structures created using the new ‘sir_graph’ program can only be downloaded individually.

Single strand frequency.
The ss-count file contains explicit statistics on single strandedness of each base in all of the computed foldings. The first line contains the number of foldings that were computed. The ith subsequent line contains the number of the ith base, the base itself and the number of times it is single stranded in all of the computed foldings. The ‘View plot’ button creates a plot of base number versus ss-count. When the ‘Averaging window’ is m>1, then the value that is plotted for the ith base is the average over the m bases centered at i. Fewer bases may be used near the ends of the plot. The magnification factor enables the user to zoom in on a region of interest.

Structure output.
The hyperlink to ‘Structure i’ leads to a portion of a single html file that contains the ith structure in an easy to read text format. This format may be suitable for viewing small foldings. It may easily be pasted into text files. If the sequence is very large and the structure plot has low resolution, it might be more convenient to view a portion of the text file than to magnify a portion of the structure plot.

Selecting PostScript,.ct file, RnaViz ct, Mac ct, GCG connect or XRNA ss will lead to single structures in the designated format. Following the png or jpg link is similar to following the png or jpg links when viewing the energy dot plot. Both lead to a new page containing a png or jpg image of the structure plot and buttons that allow the user to interact with or redraw the folding. It is called the primary structure plot page. A small ‘Loop Free-Energy Decomposition’ window, 550x150 (pixels), opens in the top left corner of the screen. With the ‘Click on a base pair to view related structural details in the Loop Free-Energy Decomposition window’ radio button selected (the default), clicking (left click) on the image will cause the details of the selected base pair to appear in the small window. The window will show the type of loop, including a stack, that is closed by the selected base pair. The free energy of that loop is given, and the specific identity of the base pair is given. If the selected point does not unambiguously select a base pair, then the auxiliary window will print a message instructing the user to try again.

Selecting the ‘Click on the image and redraw the structure with options below’ button activates the text area and four option buttons below. Clicking on a point in the structure plot will cause a new image to be created. If the magnification factor is chosen to be >1, then the new image will be magnified about the selected point; otherwise the entire structure will be redrawn. The ‘annotation’ button allows the user to select p-num or ss-count annotation, even if these were not selected initially. The ‘annotation type’ button has no effect unless an annotation option is selected. In this case, an annotated structure will be drawn with colored base letters, with colored dots or with colored dots containing base letters, corresponding to the choices ‘Character’, ‘Dot’ or ‘Both’, respectively. The ‘Image Width’ controls the resolution (size) of the redrawn image, as explained above. If the PostScript ‘Output’ option is chosen, then a static PostScript image is created. If the png or jpg ‘Output’ option is chosen, then clicking on the plot will cause a new window to open containing the redrawn image. There are no buttons on this page. The ‘Color Table’ option leads to a table that explicitly shows the colors used in the structure plot versus the p-num or ss-count values. The p-num or ss-count file hyperlink leads to the corresponding raw text files. This is called the secondary structure window. New png or jpg images generated from the primary dot plot window will refresh the secondary window. A typical use might be to magnify about a region in the initial structure plot so that a portion of the structure might be printed at a suitable resolution. It could also help the user select base pairs with precision. Closing the secondary window will not affect other windows. Closing the primary window will also close the ‘Loop Free-Energy Decomposition’ window.

New form of structure output.
The structure plots that are drawn by default and that can be downloaded at once are generated for the most part by the ‘naview’ program by Bruccoleri and Heinrich (39). The ‘naview’ program creates a device independent ASCII plot file (with suffix plt2). Software developed for the mfold package at Washington University creates PostScript, jpg and png images from these files. A new drawing program named ‘sir_graph’ was created at Washington University for creating displays of nucleic acid secondary structure. It is an interactive program that runs on a variety of Unix platforms. A non-interactive version, ‘sir_graph_ng’, creates displays for the mfold web server. The ‘Click Here for New Structure Viewing Options’ hyperlink on the primary results page leads to a ‘Structure Viewer’ page that is equivalent to the primary window page described above. This page is identical to the ones that may be reached by following the ‘new’ hyperlinks in the ‘View Individual Structures’ section.

The user has a wider selection of choices with this new software. The external loop, or exterior base pairs, may be drawn in the default way, arranged around a circle, or may be drawn ‘Flat’ with all the helices that meet it being parallel to one another. The ‘Flat_Alt’ option draws successive stems at 180° angles from the previous. The background may be the usual white or black. The base pair symbol may be a ‘dot’ which is usual for the mfold web server, or a ‘line’, which is common in most secondary structure displays. The ‘Algorithm’ may be the ‘Default’ or ‘Simple’ type, both of which avoid overlaps of stems. The ‘Default’ and ‘Simple’ algorithms draw the multi-branch loops differently. The ‘Natural’ option draws every loop as a perfect circle, with the stems coming out at appropriate angles. This may produce pleasing results for small foldings, but will usually produce a hopelessly tangled mess for large foldings. As with the display of structures using the older naview program, the ‘Loop Free-Energy Decomposition’ window opens and allows the user to identify precise loops or stacks, together with their free energies and closing base pairs. Also, the secondary structure window page is available for magnifications about chosen points.

Thermodynamic details.
When temperature cannot be altered, the ‘Thermodynamic Details’ link leads to a portion of an html file that contains the entire decomposition of the particular folding into loops and stacks, together with their free energies and closing base pairs. Consecutive runs of base pairs are summarized as helices. Small portions of this file are displayed in the ‘Loop Free-Energy Decomposition’ window when the user is interacting with a png or jpg plot of a folding.

When temperature may be altered, then a more attractive page appears that gives the free energy, enthalpy, entropy and an estimated T_m. T_m is computed using a simple 2-state model. This assumes that the molecular is either folded as shown or else completely single stranded. Such an assumption is reasonable for short molecules. For larger molecules, it might be useful to refold near the predicted T_m. If the new T_m is ‘significantly’ larger, it indicates that the given folding can rearrange into another folding with a higher melting temperature. The server uses T_m=H/S, where H and S are the structure enthalpy and entropy, respectively. (Because of the units used for entropy, the actual formula is T_m=1000H/S-273.15 in °C.) The usual details table follows.

Structure dot plot.
The ‘Dot plot folding comparisons’ option allows the user to view any subset of the computed foldings in a dot plot. This option is available when two or more foldings are computed. The user may select foldings to be displayed by selecting the corresponding radio buttons. The default is the display of the first two foldings. For convenience, an ‘All’ button has been supplied.

Base pairs that occur in all foldings are colored black. Those that occur in two or more, but not all foldings are colored grey. Otherwise, base pairs are assigned a unique color that depends on the structure. This system breaks down when more than 15 structures are computed.

Selecting the png or jpg option in the first button and pressing the ‘Do the Comparison’ option leads to a primary structure dot plot page that is very similar to the primary energy dot plot page. Clicking on a base pair with the default ‘Click on a base pair to identify’ radio button selected will show the identity of the base pair and list which of the computed structures contain that base pair. When the ‘Click on the image and redraw the dot plot with options below’ radio button is selected, a click on the dot plot will generate a new dot plot. The dot plot will be in a ‘secondary’ structure dot plot window if the png or jpg option is selected. This is very similar to the secondary dot plot window for the energy dot plot. The magnification, output format and image width parameters are the usual ones explained above for the energy dot plot. The ‘Multicolor Overlap’ option that may be selected will have the effect of drawing the otherwise grey dots in a multi-color mode that displays precisely what foldings contain that base pair. In the mfold web server, ‘dots’ are really plotted as (usually small) squares. A collection of trapezoids in different colors are drawn within the square region. These indicate which foldings contain that base pair. The dot plot must be drawn with a high enough magnification for this feature to be seen. It is not necessary when only two foldings have been computed. When three foldings are computed, the gray base pairs that are in structures 1&2, 1&3 and 2&3 are plotted in the lower left triangle in three distinct colors, so they may be identified without using the multicolor option.

The other servers
The ‘quikfold’ server: folding many sequences at once.
Many users are not interested in the sophisticated viewing options that are available on the mfold web server. Some may be content with the easy to read text format or with the ‘ct format’ files that may be downloaded and used locally to create secondary structure plots. Many wish to fold many sequences at once under the same conditions. For example, a user may wish to predict foldings and melting temperatures for hundreds of short molecular beacons.

For this reason, the ‘quikfold server’ was created. It uses the same ‘nafold2’ program that folds sequences in the regular server, but in a different, multiple molecule mode. Only immediate jobs are allowed, with each sequence containing no more than 600 nucleotides. The theoretical upper limit for the number of sequences that may be submitted at once is 25 000, but I would not advise submitting nearly that many! The upper limit depends on the Internet, browser settings (such as timeout) and other factors. I do know, for example, that 1000 sequences of length 100 each should pose no problem.

The input page contains buttons and text areas that have already been defined. The single new item is that one extra character is required for the sequence input box. Each of the input sequences must be separated by (at least) one semi-colon, ‘;’. All other characters, except for ‘A–Z’ and ‘a–z’, are removed in a preprocessing step.

The quikfold server gives text, ct, RNAML and thermodynamic details as output. In each case, the output is in a single, possibly large text file. The sequences are named by taking the given or default sequence name, and adding ‘_i’ to each name, where i ranges from 1 to the number of sequences. The folding results are erased within 40 min to several hours. An individual sequence that cannot fold will not abort the entire job. It will simply be skipped over in the output.

The ‘zipfold’ server: prediction of minimum free energy only.
Some users want only the minimum folding energy of a sequence. For this reason, the zipfold server was created. The input page is similar to that of the quikfold server. Sequence lengths may be up to 800. In one instance, a user submitted many thousands of jobs; each containing 50 randomized sequences of length 500. The failure rate was less than 1%.

Because speed is important and so little information is being requested, the underlying code was simplified to run much faster. At the time of this writing, 378 tRNAs were processed in just over 11 s. The actual performance depends on the server load and difficult to control Internet related factors.

The underlying code has been stripped of its ability to handle constraints. Moreover, the ‘fill algorithm’ has been truncated so that only an optimal energy or a single optimal structure may be computed. Newer, experimental versions, streamline the sequence input as well.

All the output appears on a single, primary results page. The G values are arranged in a list. Those sequences that cannot fold are given a large positive folding energy. The current value is 10 000.

The ‘T_m’ server for single stranded nucleic acids: prediction of two state melting temperature.
Some users want a little more that just a minimum folding energy. They also desire an estimated T_m. In this case, the T_m server may be used. Only version 2.3 RNA folding parameters are available for this server, since enthalpies are required for estimating T_m. The input is the same as with the zipfold server.

The T_m server uses the same simplified folding program as the zipfold server, except that a minimum energy folding is computed in this case. The enthalpy, H, of this folding is then computed using the appropriate nearest neighbor parameters. From this, it is easy to compute an entropy, S and then T_m, using a 2-state model as described above.

In fact, this server was created to ‘service’ the ‘OligoArray 1’ program of Rouillard et al. (40). The OligoArray program ‘hits’ the server with single DNA sequences using a direct request created by a Java application. Subsequent versions of OligoArray will abandon this inefficient procedure.

The hybridization server: hybridization of two strands.
The underlying mfold software has had numerous small changes added over the years to accommodate ‘special requests’ or applications. Many users have used the regular mfold software, either in ‘stand alone’ form or on the web, to simulate the hybridization of two strands of RNA or DNA. They have done so by taking two sequences, A=a₁a₂a₃···a_m and B=b₁b₂b₃···b_m and creating a single sequence, S=a₁a₂a₃···a_mn₁n₂···n_kb₁b₂b₃···b_m, by concatenating A with B using some non-pairing characters, n₁···n_k as a linker.

There are two problems with this approach. The first, less serious problem, is that the k linker residues end up in a loop of some sort. If A and B hybridize perfectly, then this loop will almost certainly be a hairpin loop. An example of such a folding is given in Figure 2. The problem is that an incorrect hairpin free energy is applied to loop ‘2’. Instead, it should be treated as an exterior loop, like ‘1’. In addition, an initiation free energy, I, needs to be added. This is 4.1 kcal/mol for RNA at 37°C and 1.96 kcal/mol for DNA at 37°C.

View larger version (9K): [in this window] [in a new window]   Figure 2. Two complementary DNA sequences have been concatenated by linking them with six Ns. The two loops are labeled ‘1’ and ‘2’. Loop ‘1’ is an exterior loop, which is correct, but loop ‘2’ is a hairpin loop instead of an exterior loop.

 
In order to accommodate such situations, the following feature was added to mfold. When three consecutive Ls occur in a sequence (L for Linker), they are recognized as a linker. If they occur in a hairpin loop, then this loop is treated instead as an exterior loop and I is added as well. The output still contains the linker residues, but the value of G is now correct.

The linker residues may be forced to be in a hairpin loop by forcing a simple hybridization of A with B. This may be accomplished by using the two constraints

These constraints forbid all intramolecular base pairs, so that the only allowed base pairs link A with B, forcing, among other things, the linker residues to be in a hairpin loop.

The more serious problem, however, is that any estimated T_m is nonsense, since total nucleic acid concentration, C, must be taken into account when two strands hybridize. The correct 2-state estimate for T_m is

where f=2 if AB and f=4 if A=B.

All of this is accomplished automatically by the hybridization server. The two sequences, A and B, are entered into the sequence box text area, separated by a semi-colon, ‘;’. The total nucleic acid concentration must be given. The result is a simple output page containing the job ID, G, H, S and T_m in a single line. No units are given, although they are understood to be kcal/mol, kcal/mol, cal/(K · mol) and °C, respectively.

Access
The mfold web server is open to all users. No restrictions are applied to commercial users. However, users should be aware that the server is not secure and that data flowing both in and out may be detected by others. Moreover, query information is stored indefinitely in log files on the server. These log files are treated as confidential information, although gross statistics on usage are collected and disseminated. Furthermore, some of the submissions are selected as examples for teaching, but only if database searches reveal that the sequence is already in a public database.

The regular mfold server may not be hit by ‘GET’ requests. Only the ‘POST’ operation is permitted. The other applications that have been described may also be hit with ‘GET’ requests. Hitting these other applications directly using, for example, Java applications is permitted and will continue to be permitted until such time as the server is overloaded. This has not happened yet. At Washington University, perhaps 100–200 folding jobs were submitted each work day, with far fewer on weekends. For the month of December, 2002 (for example), I counted 21 688 submissions to the regular mfold web server, 3128 to quikfold, 35 904 to zipfold, 193 428 to T_m and 542 to the hybridization server.

	EQUIPMENT AND ORGANIZATION

TOP ABSTRACT INTRODUCTION SERVER CONTENT AND ORGANIZATION EQUIPMENT AND ORGANIZATION FUTURE DIRECTIONS CITING THE MFOLD WEB... SUPPLEMENTARY MATERIAL REFERENCES

 
The current web server is running on equipment donated to RPI by IBM Research in the fall of 2001. The detailed hardware specifications are given in Table 5. The total value of the donated equipment comes to about $214 000. We contributed another $5000–$10 000 for a switch and a custom designed power source. All equipment was originally assembled and housed at the Academy of Electronic Media (http://www.academy.rpi.edu/). The server is expected to be moved to the Voorhees Computing Center during the spring or summer of 2003.

View this table: [in this window] [in a new window]   Table 5. The mfold web server is a major application on the RPI-Wadsworth Bioinformatics Center web site. It is housed on a high-performance web server consisting of a dual 1 GHz Pentium processor with 4 Gb memory and a 73 Gb disk, connected to a cluster of 35 other dual 1 GHz Pentiums, each with 1 Gb of memory and a 36 Gb disk

 

	FUTURE DIRECTIONS

TOP ABSTRACT INTRODUCTION SERVER CONTENT AND ORGANIZATION EQUIPMENT AND ORGANIZATION FUTURE DIRECTIONS CITING THE MFOLD WEB... SUPPLEMENTARY MATERIAL REFERENCES

 
It is our intention to keep the RNAML output up-to-date and in accord with the DTD file: http://www-lbit.iro.umontreal.ca/rnaml/current/rnaml.dtd. Within the next month or two, the RNAML output will contain the ‘x, y’ co-ordinates that correspond to the structure plots produced by the ‘sir_graph’ program, when they are available. The next year or two should see the incorporation of important, novel applications that are currently being developed.

Hybridization with partition functions
New algorithms for the folding and hybridization of two separate strands of RNA or DNA have been created (41). Initial versions of the software are now being carefully reprogrammed and updated. We are investigating theoretical and experimental ways to improve the accuracy of predications. We are also beginning to consider dealing with the thermodynamics of hybridization of probes to nucleic acids immobilized on chips (42).

Database searching
We are developing a new ‘BLAST-like’ algorithm for searching nucleic acid databases that is based on the computer science technique of hashing, for speed, but searches for regions of best complementarity between two sections of DNA or RNA. Scoring is based on nearest neighbor free energy parameters for RNA or DNA. An early version is now being used to search for putative binding sites for micro-RNAs (miRNAs) in databases of mRNAs, and to search for (undesirable) alternative binding sites for gene specific DNA probes.

Development of user interface
The current user interface for the mfold web server is now several years old and would benefit from either partial or total rewriting. The submission of data through forms is still reasonable. A ‘browse’ option for uploading sequence files could easily be added. It would be useful to make it easier to enter constraint information. Even more valuable would be a preprocessor that would check the constraint information for consistency and reasonableness. Up to now, many constrained foldings have failed because of minor errors in constraints.

The output pages could benefit substantially by being redesigned. As a first step, the relatively new ‘sir_graph’ program should replace the older naview program for creating secondary structure plots. It would be desirable to create an interactive structure drawing feature. Such a feature already exists in the stand alone version of the ‘sir_graph’ program. We are fortunate here at RPI to have the cutting edge Academy of Electronic Media; a group that might collaborate with us in creating a really first rate user interface for the web server.

	CITING THE MFOLD WEB SERVER

TOP ABSTRACT INTRODUCTION SERVER CONTENT AND ORGANIZATION EQUIPMENT AND ORGANIZATION FUTURE DIRECTIONS CITING THE MFOLD WEB... SUPPLEMENTARY MATERIAL REFERENCES

 
Authors who make use of the mfold web server should cite this article as a general reference and should also include the URL to the entrance page, http://www.bioinfo.rpi.edu/applications/mfold. The web server pages will list additional articles for citation that relate to the free energy parameters that are used and the underlying software.

	SUPPLEMENTARY MATERIAL

TOP ABSTRACT INTRODUCTION SERVER CONTENT AND ORGANIZATION EQUIPMENT AND ORGANIZATION FUTURE DIRECTIONS CITING THE MFOLD WEB... SUPPLEMENTARY MATERIAL REFERENCES

 
Supplementary Material is available at NAR Online and can also be found at MFOLDROOT/www-NAR03/supp.

	ACKNOWLEDGEMENTS

 
This work was supported, in part, by grant #GM54250 from the National Institutes of Health. I thank Gerald Johns for hardware and systems support in the early days of the server at Washington University. All of the graphics applications (except for naview) were written by Darrin Stewart, who also wrote the ‘cgi’ scripts for interactive viewing of plots. I thank Art Sanderson (Vice President of Research at RPI), for connecting me with the Academy of Electronic Media and for supporting this project; Bill Shumway, for initiating and facilitating interactions with IBM Research; and Alex Yu, who has done so much work in assembling the hardware, organizing the web server layout, porting applications and in keeping the server running day in and day out. Finally, I thank IBM Research for the SUR grant that gave us a thirty-fold increase in computer power for housing this valuable resource.

	REFERENCES

TOP ABSTRACT INTRODUCTION SERVER CONTENT AND ORGANIZATION EQUIPMENT AND ORGANIZATION FUTURE DIRECTIONS CITING THE MFOLD WEB... SUPPLEMENTARY MATERIAL REFERENCES

 

1. Doty,P., Boedtker,H., Fresco,J.R., Haselkorn,R. and Litt,M. (1959) Secondary structure in ribonucleic acids. Proc. Natl Acad. Sci. USA, 45, 482–499.[Free Full Text]

2. Fresco,J.R., Alberts,B.M. and Doty,P. (1960) Some molecular details of the secondary structure of ribonucleic acid. Nature, 188, 98–101.[CrossRef][Medline]

3. Borer,P.N., Dengler,B., Tinoco,I.,Jr and Uhlenbeck,O.C. (1974) Stability of ribonucleic acid double-stranded helices. J. Mol. Biol., 86, 843–853.[CrossRef][Web of Science][Medline]

4. Tinoco,I.,Jr and Uhlenbeck,O.C. (1971) Estimation of secondary structure in ribonucleic acids. Nature, 230, 362–367.[CrossRef][Medline]

5. Tinoco,I.,Jr, Borer,P.N., Dengler,B., Levine,M.D., Uhlenbeck,O.C., Crothers,D.M. and Gralla,J. (1973) Improved estimation of secondary structure in ribonucleic acids. Nature New Biol., 246, 40–41.[Web of Science][Medline]

6. Uhlenbeck,O.C., Borer,P.N., Dengler,B. and Tinoco,I.,Jr (1973) Stability of RNA hairpin loops: A₆-C_m-U₆. J. Mol. Biol., 73, 483–496.[CrossRef][Web of Science][Medline]

7. Delisi,C. and Crothers,D.M. (1971) Prediction of RNA secondary structure. Proc. Natl Acad. Sci. USA, 68, 2682–2685.[Abstract/Free Full Text]

8. Waterman,M.S. and Smith,T.M. (1978) RNA secondary structure: a complete mathematical analysis. Math. Biosci., 42, 257–266.[CrossRef][Web of Science]

9. Waterman,M.S. (1978) Secondary structure of single-stranded nucleic acids. In Rota,G.-C. (ed.), Studies in Foundations and Combinatorics number 1 in Advances in Mathematics, Supplementary Studies. Academic Press, NY, pp. 167–212.

10. Nussinov,R., Pieczenik,G., Griggs,J.R. and Kleitman,D.J. (1978) Algorithm for loop matchings. SIAM J. Appl. Math., 35, 68–82.[CrossRef]

11. Nussinov,R. and Jacobson,A.B. (1980) Fast algorithm for predicting the secondary structure of single-stranded RNA. Proc. Natl Acad. Sci. USA, 77, 6309–6313.[Abstract/Free Full Text]

12. Zuker,M. and Stiegler,P. (1981) Optimal computer folding of large RNA sequences using thermodynamics and auxiliary information. Nucleic Acids Res., 9, 133–148.[Abstract/Free Full Text]

13. Sankoff,D., Kruskal,J.B., Mainville,S. and Cedergren,R.J. (1983) Fast algorithms to determine RNA secondary structures containing multiple loops. chapter 3. In Sankoff,D. and Kruskal,J.B. (eds), Time Warps, String Edits, and Macromolecules: The Theory and Practice of Sequence Comparison, Addison-Wesley Reading, MA, pp. 93–120.

14. Zuker,M. (1989) Computer prediction of RNA structure. Methods Enzymol., 180, 262–288.[Web of Science][Medline]

15. Devereux,J., Haeberli,P. and Smithies,O. (1984) A comprehensive set of sequence analysis programs for the VAX. Nucleic Acids Res., 12, 387–395.[Abstract/Free Full Text]

16. Salser,W. (1977) Globin mRNA sequences: analysis of base pairing and evolutionary implications. Cold Spring Harbor Symp. Quant. Biol., 42, 985–1002.[Web of Science]

17. Freier,S.M., Kierzek,R., Jaeger,J.A., Sugimoto,N., Caruthers,M.M., Neilson,T. and Turner,D.H. (1986) Improved free-energy parameters for predictions of RNA duplex stability. Proc. Natl Acad. Sci. USA, 83, 9373–9377.[Abstract/Free Full Text]

18. Zuker,M. (1989) On finding all suboptimal foldings of an RNA molecule. Science, 244, 48–52.[Abstract/Free Full Text]

19. Jaeger,J.A., Turner,D.H. and Zuker,M. (1989) Improved predictions of secondary structures for RNA. Proc. Natl Acad. Sci. USA, 86, 7706–7710.[Abstract/Free Full Text]

20. Jaeger,J.A., Turner,D.H. and Zuker,M. (1990) Predicting optimal and suboptimal secondary structure for RNA. Methods Enzymol., 183, 281–306.[Web of Science][Medline]

21. Zuker,M. (1994) Prediction of RNA secondary structure by energy minimization. chapter 23. In Griffin,A.M. and Griffin,H.G. (eds), Computer Analysis of Sequence Data, 25, Part II, Humana Press, Inc., Totowa, NJ, pp. 267–294.[CrossRef]

22. Zuker,M., Mathews,D.H. and Turner,D.H. (1999) Algorithms and thermodynamics for RNA secondary structure prediction: A practical guide. In Barciszewski,J. and Clark,B.F.C. (ed.), RNA Biochemistry and Biotechnology, number 70 in NATO Science Partnership Sub-Series: 3: High Technology. chapter 2, Kluwer Academic Publishers Dordrecht, The Netherlands, pp. 11–43.

23. Mathews,D.H., Andre,T.C., Kim,J., Turner,D.H. and Zuker,M. (1998) An updated recursive algorithm for RNA secondary structure prediction with improved free energy parameters. chapter 15. In Leontis,N.B. and SantaLucia,J.,Jr (eds), American Chemical Society Symposium Series 682, American Chemical Society Washington, DC, pp. 246–257.

24. Mathews,D.H., Turner,D.H. and Zuker,M. (2000) RNA secondary structure prediction. chapter 11.2. In Beaucage,S., Bergstrom,D.E., Glick,G.D. and Jones,R.A. (eds), Current Protocols in Nucleic Acid Chemistry, John Wiley & Sons New York, NY, pp. 1–10.

25. Walter,A.E., Turner,D.H., Kim,J., Lyttle,M.H., Müller,P., Mathews,D.H. and Zuker,M. (1994) Coaxial stacking of helixes enhances binding of oligoribonucleotides and improves predictions of RNA folding. Proc. Natl Acad. Sci. USA, 91, 9218–9222.[Abstract/Free Full Text]

26. Mathews,D.H., Sabina,J., Zuker,M. and Turner,D.H. (1999) Expanded sequence dependence of thermodynamic parameters improves prediction of RNA secondary structure. J. Mol. Biol., 288, 911–940.[CrossRef][Web of Science][Medline]

27. SantaLucia,J.,Jr (1998) A unified view of polymer, dumbell, and oligonucleotide DNA nearest-neighbor thermodynamics. Proc. Natl Acad. Sci. USA, 95, 1460–1465.[Abstract/Free Full Text]

28. Allawi,H.T. and SantaLucia,J.,Jr (1997) Thermodynamics and NMR of internal G·T mismatches in DNA. Biochemistry, 36, 10581–10594.[CrossRef][Medline]

29. Allawi,H.T. and SantaLucia,J.,Jr (1998) Nearest neighbor thermodynamic parameters for internal G·A mismatches in DNA. Biochemistry, 37, 2170–2179.[CrossRef][Medline]

30. Allawi,H.T. and SantaLucia,J.,Jr (1998) Thermodynamics of internal C·T mismatches in DNA. Nucleic Acids Res., 26, 2694–2701.[Abstract/Free Full Text]

31. Allawi,H.T. and SantaLucia,J.,Jr (1998) Nearest-neighbor thermodynamics of internal A·C mismatches in DNA: Sequence dependence and pH effects. Biochemistry, 37, 9435–9444.[CrossRef][Medline]

32. Peyret,N., Senevirtne,P.A., Allawi,H.T. and SantaLucia,J.,Jr (1999) Nearest-neighbor thermodynamics and NMR of DNA sequences with internal A·A, C·C, G·G, and T·T mismatches. Biochemistry, 38, 3468–3477.[CrossRef][Medline]

33. Bommarito,S., Peyret,N. and SantaLucia,J.,Jr (2000) Thermodynamic parameters for DNA sequences with dangling ends. Nucleic Acids Res., 28, 1929–1934.[Abstract/Free Full Text]

34. Peyret,N. (2000) Prediction of Nucleic Acid Hybridization: Parameters and Algorithms. PhD Theris, Wayne State University Department of Chemistry, Detroit, MI.

35. Cornish-Bowden,A. (1985) Nomenclature for incompletely specified bases in nucleic acid sequences: recommendations 1984. Nucleic Acids Res., 13, 3021–3030.[Free Full Text]

36. Zuker,M. and Jacobson,A. (1998) Using reliability information to annotate RNA secondary structures. RNA, 4, 669–679.[Abstract]

37. Waugh,A., Gendron,P., Altman,R., Brown,J.W., Case,D., Gautheret,D., Harvey,S.C., Leontis,N., Westbrook,J., Westhof,E. et al. (2002) RNAML: A standard syntax for exchanging RNA information. RNA, 8, 707–717.[Abstract]

38. Zuker,M. and Jacobson,A.G. (1995) ‘Well-determined’ regions in RNA secondary structure prediction. application to small and large subunit rRNA. Nucleic Acids Res., 23, 2791–2798.[Abstract/Free Full Text]

39. Bruccoleri,R. and Heinrich,G. (1988) An improved algorithm for nucleic acid secondary structure display. Comput. Appl. Biosci., 4, 167–173.[Abstract/Free Full Text]

40. Rouillard,J.-M., Herbert,C.J. and Zuker,M. (2002) OligoArray: Genome-scale oligonucleotide design for microarrays. Bioinformatics, 18, 486–487.[Abstract/Free Full Text]

41. Dimitrov,R.A. and Zuker,M. (2003) Prediction of hybridization and melting for double stranded nucleic acids. Biophys. J., in press.

42. Fotin,A.V., Drobyshev,A.L., Proudnikov,D.Y., Perov,A.N. and Mirzabekov,A.D. (1998) Parallel thermodynamic analysis of duplexes on oligodeoxyribonucleotide microchips. Nucleic Acids Res., 26, 1515–1521.[Abstract/Free Full Text]

CiteULike    Connotea    Del.icio.us    What's this?

This article has been cited by other articles:

				F. Martinez, J. Marques, M. L. Salvador, and J.-A. Daros Mutational analysis of eggplant latent viroid RNA processing in Chlamydomonas reinhardtii chloroplast J. Gen. Virol., December 1, 2009; 90(12): 3057 - 3065. [Abstract] [Full Text] [PDF]

				T. Nitta, A. Hofacre, S. Hull, and H. Fan Identification and Mutational Analysis of a Rej Response Element in Jaagsiekte Sheep Retrovirus RNA J. Virol., December 1, 2009; 83(23): 12499 - 12511. [Abstract] [Full Text] [PDF]

				L. R. Cruz-Vera, R. Yang, and C. Yanofsky Tryptophan Inhibits Proteus vulgaris TnaC Leader Peptide Elongation, Activating tna Operon Expression J. Bacteriol., November 15, 2009; 191(22): 7001 - 7006. [Abstract] [Full Text] [PDF]

				P. Friebe and R. Bartenschlager Role of RNA Structures in Genome Terminal Sequences of the Hepatitis C Virus for Replication and Assembly J. Virol., November 15, 2009; 83(22): 11989 - 11995. [Abstract] [Full Text] [PDF]

				L. Li, J. Victoria, A. Kapoor, O. Blinkova, C. Wang, F. Babrzadeh, C. J. Mason, P. Pandey, H. Triki, O. Bahri, et al. A Novel Picornavirus Associated with Gastroenteritis J. Virol., November 15, 2009; 83(22): 12002 - 12006. [Abstract] [Full Text] [PDF]

				Y. Wei, C. Qin, T. Jiang, X. Li, H. Zhao, Z. Liu, Y. Deng, R. Liu, S. Chen, M. Yu, et al. Translational Regulation by the 3' Untranslated Region of the Dengue Type 2 Virus Genome Am J Trop Med Hyg, November 1, 2009; 81(5): 817 - 824. [Abstract] [Full Text] [PDF]

				E. De Paoli, A. Dorantes-Acosta, J. Zhai, M. Accerbi, D.-H. Jeong, S. Park, B. C. Meyers, R. A. Jorgensen, and P. J. Green Distinct extremely abundant siRNAs associated with cosuppression in petunia RNA, November 1, 2009; 15(11): 1965 - 1970. [Abstract] [Full Text] [PDF]

				E. de Wit, S. E.V. Linsen, E. Cuppen, and E. Berezikov Repertoire and evolution of miRNA genes in four divergent nematode species Genome Res., November 1, 2009; 19(11): 2064 - 2074. [Abstract] [Full Text] [PDF]

				D. H. Turner and D. H. Mathews NNDB: the nearest neighbor parameter database for predicting stability of nucleic acid secondary structure Nucleic Acids Res., October 30, 2009; (2009) gkp892v1. [Abstract] [Full Text] [PDF]

				A. Vijaykumar, A. Saini, and N. Jawali Phylogenetic Analysis of Subgenus Vigna Species Using Nuclear Ribosomal RNA ITS: Evidence of Hybridization among Vigna unguiculata Subspecies J. Hered., October 27, 2009; (2009) esp084v1. [Abstract] [Full Text] [PDF]

				J. W. Snow and S. H. Orkin Translational Isoforms of FOG1 Regulate GATA1-interacting Complexes J. Biol. Chem., October 23, 2009; 284(43): 29310 - 29319. [Abstract] [Full Text] [PDF]

				B. P. Higgins, A. C. Popkowski, P. R. Caruana, and A. C. Karls Site-Specific Insertion of IS492 in Pseudoalteromonas atlantica J. Bacteriol., October 15, 2009; 191(20): 6408 - 6414. [Abstract] [Full Text] [PDF]

				R. Moreno, S. Marzi, P. Romby, and F. Rojo The Crc global regulator binds to an unpaired A-rich motif at the Pseudomonas putida alkS mRNA coding sequence and inhibits translation initiation Nucleic Acids Res., October 13, 2009; (2009) gkp825v1. [Abstract] [Full Text] [PDF]

				X. Xu, K. Tsumagari, J. Sowden, R. Tawil, A. P. Boyle, L. Song, T. S. Furey, G. E. Crawford, and M. Ehrlich DNaseI hypersensitivity at gene-poor, FSH dystrophy-linked 4q35.2 Nucleic Acids Res., October 9, 2009; (2009) gkp833v1. [Abstract] [Full Text] [PDF]

				G. Leon and P. H. Roy Potential Role of Group IIC-attC Introns in Integron Cassette Formation J. Bacteriol., October 1, 2009; 191(19): 6040 - 6051. [Abstract] [Full Text] [PDF]

				S. Liang and G. J. Connell An electrochemiluminescent aptamer switch for a high-throughput assay of an RNA editing reaction RNA, October 1, 2009; 15(10): 1929 - 1938. [Abstract] [Full Text] [PDF]

				Y. Kino, C. Washizu, Y. Oma, H. Onishi, Y. Nezu, N. Sasagawa, N. Nukina, and S. Ishiura MBNL and CELF proteins regulate alternative splicing of the skeletal muscle chloride channel CLCN1 Nucleic Acids Res., October 1, 2009; 37(19): 6477 - 6490. [Abstract] [Full Text] [PDF]

				M. L. Yeung, Y. Bennasser, K. Watashi, S.-Y. Le, L. Houzet, and K.-T. Jeang Pyrosequencing of small non-coding RNAs in HIV-1 infected cells: evidence for the processing of a viral-cellular double-stranded RNA hybrid Nucleic Acids Res., October 1, 2009; 37(19): 6575 - 6586. [Abstract] [Full Text] [PDF]

				L. Leclere, P. Schuchert, C. Cruaud, A. Couloux, and M. Manuel Molecular Phylogenetics of Thecata (Hydrozoa, Cnidaria) Reveals Long-Term Maintenance of Life History Traits despite High Frequency of Recent Character Changes Syst Biol, October 1, 2009; 58(5): 509 - 526. [Abstract] [Full Text] [PDF]

				K. A. Corbin-Lickfett, I-H. B. Chen, M. J. Cocco, and R. M. Sandri-Goldin The HSV-1 ICP27 RGG box specifically binds flexible, GC-rich sequences but not G-quartet structures Nucleic Acids Res., September 26, 2009; (2009) gkp793v1. [Abstract] [Full Text] [PDF]

				G. Deikus and D. H. Bechhofer Bacillus subtilis trp Leader RNA: RNase J1 ENDONUCLEASE CLEAVAGE SPECIFICITY AND PNPase PROCESSING J. Biol. Chem., September 25, 2009; 284(39): 26394 - 26401. [Abstract] [Full Text] [PDF]

				J. A. Mertz, A. B. Chadee, H. Byun, R. Russell, and J. P. Dudley Mapping of the Functional Boundaries and Secondary Structure of the Mouse Mammary Tumor Virus Rem-responsive Element J. Biol. Chem., September 18, 2009; 284(38): 25642 - 25652. [Abstract] [Full Text] [PDF]

				J. M. Peters, R. A. Mooney, P. F. Kuan, J. L. Rowland, S. Keles, and R. Landick Rho directs widespread termination of intragenic and stable RNA transcription PNAS, September 8, 2009; 106(36): 15406 - 15411. [Abstract] [Full Text] [PDF]

				Y. Zhang and G. Hong Post-transcriptional regulation of NifA expression by Hfq and RNase E complex in Rhizobium leguminosarum bv. viciae Acta Biochim Biophys Sin, September 1, 2009; 41(9): 719 - 730. [Abstract] [Full Text] [PDF]

				R. S. P. Horler and C. K. Vanderpool Homologs of the small RNA SgrS are broadly distributed in enteric bacteria but have diverged in size and sequence Nucleic Acids Res., September 1, 2009; 37(16): 5465 - 5476. [Abstract] [Full Text] [PDF]

				Y. Terui, K. Higashi, Y. Tabei, H. Tomitori, K. Yamamoto, A. Ishihama, K. Igarashi, and K. Kashiwagi Enhancement of the Synthesis of RpoE and StpA by Polyamines at the Level of Translation in Escherichia coli under Heat Shock Conditions J. Bacteriol., September 1, 2009; 191(17): 5348 - 5357. [Abstract] [Full Text] [PDF]

				S. Chirayil, R. Chirayil, and K. J. Luebke Discovering ligands for a microRNA precursor with peptoid microarrays Nucleic Acids Res., September 1, 2009; 37(16): 5486 - 5497. [Abstract] [Full Text] [PDF]

				A. M. Wahid, V. K. Coventry, and G. L. Conn The PKR-binding domain of adenovirus VA RNAI exists as a mixture of two functionally non-equivalent structures Nucleic Acids Res., September 1, 2009; 37(17): 5830 - 5837. [Abstract] [Full Text] [PDF]

				A. Martin, H. Adam, M. Diaz-Mendoza, M. Zurczak, N. D. Gonzalez-Schain, and P. Suarez-Lopez Graft-transmissible induction of potato tuberization by the microRNA miR172 Development, September 1, 2009; 136(17): 2873 - 2881. [Abstract] [Full Text] [PDF]

				L. S. Spruill and P. J. McDermott Role of the 5'-untranslated region in regulating translational efficiency of specific mRNAs in adult cardiocytes FASEB J, September 1, 2009; 23(9): 2879 - 2887. [Abstract] [Full Text] [PDF]

				Z. Ma, W.-O. Lui, A. Fire, and S. S. Dadras Profiling and Discovery of Novel miRNAs from Formalin-Fixed, Paraffin-Embedded Melanoma and Nodal Specimens J. Mol. Diagn., September 1, 2009; 11(5): 420 - 429. [Abstract] [Full Text] [PDF]

				J. S. Lin and K. P. McNatty Aptamer-Based Regionally Protected PCR for Protein Detection Clin. Chem., September 1, 2009; 55(9): 1686 - 1693. [Abstract] [Full Text] [PDF]

				S.-I. Yun, Y.-J. Choi, B.-H. Song, and Y.-M. Lee 3' cis-Acting Elements That Contribute to the Competence and Efficiency of Japanese Encephalitis Virus Genome Replication: Functional Importance of Sequence Duplications, Deletions, and Substitutions J. Virol., August 15, 2009; 83(16): 7909 - 7930. [Abstract] [Full Text] [PDF]

				A. Phillips and L. Cardelli A programming language for composable DNA circuits J R Soc Interface, August 6, 2009; 6(Suppl_4): S419 - S436. [Abstract] [Full Text] [PDF]

				L. Bernad, N. Duran-Vila, and S. F. Elena Effect of citrus hosts on the generation, maintenance and evolutionary fate of genetic variability of citrus exocortis viroid J. Gen. Virol., August 1, 2009; 90(8): 2040 - 2049. [Abstract] [Full Text] [PDF]

				C. H. Logue, C. F. Bosio, T. Welte, K. M. Keene, J. P. Ledermann, A. Phillips, B. J. Sheahan, D. J. Pierro, N. Marlenee, A. C. Brault, et al. Virulence variation among isolates of western equine encephalitis virus in an outbred mouse model J. Gen. Virol., August 1, 2009; 90(8): 1848 - 1858. [Abstract] [Full Text] [PDF]

				G. Leon and P. H. Roy Group IIC intron mobility into attC sites involves a bulged DNA stem-loop motif RNA, August 1, 2009; 15(8): 1543 - 1553. [Abstract] [Full Text] [PDF]

				T. Tezuka, H. Hara, Y. Ohnishi, and S. Horinouchi Identification and Gene Disruption of Small Noncoding RNAs in Streptomyces griseus J. Bacteriol., August 1, 2009; 191(15): 4896 - 4904. [Abstract] [Full Text] [PDF]

				Y. Luo and J. D. Helmann Extracytoplasmic Function {sigma} Factors with Overlapping Promoter Specificity Regulate Sublancin Production in Bacillus subtilis J. Bacteriol., August 1, 2009; 191(15): 4951 - 4958. [Abstract] [Full Text] [PDF]

				Y. Morita, C. Gilmour, D. Metcalf, and K. Poole Translational Control of the Antibiotic Inducibility of the PA5471 Gene Required for mexXY Multidrug Efflux Gene Expression in Pseudomonas aeruginosa J. Bacteriol., August 1, 2009; 191(15): 4966 - 4975. [Abstract] [Full Text] [PDF]

				M. E. Christiansen and B. M. Znosko Thermodynamic characterization of tandem mismatches found in naturally occurring RNA Nucleic Acids Res., August 1, 2009; 37(14): 4696 - 4706. [Abstract] [Full Text] [PDF]

				K. Kikuchi, M. Fukuda, T. Ito, M. Inoue, T. Yokoi, S. Chiku, T. Mitsuyama, K. Asai, T. Hirose, and Y. Aizawa Transcripts of unknown function in multiple-signaling pathways involved in human stem cell differentiation Nucleic Acids Res., August 1, 2009; 37(15): 4987 - 5000. [Abstract] [Full Text] [PDF]

				G. Kollberg, M. Tulinius, A. Melberg, N. Darin, O. Andersen, D. Holmgren, A. Oldfors, and E. Holme Clinical manifestation and a new ISCU mutation in iron-sulphur cluster deficiency myopathy Brain, August 1, 2009; 132(8): 2170 - 2179. [Abstract] [Full Text] [PDF]

				C. Johnson, A. Kasprzewska, K. Tennessen, J. Fernandes, G.-L. Nan, V. Walbot, V. Sundaresan, V. Vance, and L. H. Bowman Clusters and superclusters of phased small RNAs in the developing inflorescence of rice Genome Res., August 1, 2009; 19(8): 1429 - 1440. [Abstract] [Full Text] [PDF]

				S. Nakamura A Novel Virtual Spectrometry: Visualized Regulatory Motifs on ADM, rPol{beta} and CD83 mRNAs in Human-friendly Manners J. Biochem., August 1, 2009; 146(2): 251 - 261. [Abstract] [Full Text] [PDF]

				A. Schmidt, T. Schwerd, W. Hamm, J. C. Hellmuth, S. Cui, M. Wenzel, F. S. Hoffmann, M.-C. Michallet, R. Besch, K.-P. Hopfner, et al. 5'-triphosphate RNA requires base-paired structures to activate antiviral signaling via RIG-I PNAS, July 21, 2009; 106(29): 12067 - 12072. [Abstract] [Full Text] [PDF]

				A. Leiherer, C. Ludwig, and R. Wagner Uncoupling Human Immunodeficiency Virus Type 1 gag and pol Reading Frames: Role of the Transframe Protein p6* in Viral Replication J. Virol., July 15, 2009; 83(14): 7210 - 7220. [Abstract] [Full Text] [PDF]

				C. Romero-Lopez, R. Diaz-Gonzalez, A. Barroso-delJesus, and A. Berzal-Herranz Inhibition of hepatitis C virus replication and internal ribosome entry site-dependent translation by an RNA molecule J. Gen. Virol., July 1, 2009; 90(7): 1659 - 1669. [Abstract] [Full Text] [PDF]

				C. H. Williams, M. Panayiotou, G. D. Girling, C. I. Peard, S. Oikarinen, H. Hyoty, and G. Stanway Evolution and conservation in human parechovirus genomes J. Gen. Virol., July 1, 2009; 90(7): 1702 - 1712. [Abstract] [Full Text] [PDF]

				T. J. Gardiner, C. P. Christov, A. R. Langley, and T. Krude A conserved motif of vertebrate Y RNAs essential for chromosomal DNA replication RNA, July 1, 2009; 15(7): 1375 - 1385. [Abstract] [Full Text] [PDF]

				A. Ely, T. Naidoo, and P. Arbuthnot Efficient silencing of gene expression with modular trimeric Pol II expression cassettes comprising microRNA shuttles Nucleic Acids Res., July 1, 2009; 37(13): e91 - e91. [Abstract] [Full Text] [PDF]

				S. J. Schroeder Advances in RNA Structure Prediction from Sequence: New Tools for Generating Hypotheses about Viral RNA Structure-Function Relationships J. Virol., July 1, 2009; 83(13): 6326 - 6334. [Full Text] [PDF]

				M. Vargiolu, D. Fusco, I. Kurelac, D. Dirnberger, R. Baumeister, I. Morra, A. Melcarne, R. Rimondini, G. Romeo, and E. Bonora The Tyrosine Kinase Receptor RET Interacts in Vivo with Aryl Hydrocarbon Receptor-Interacting Protein to Alter Survivin Availability J. Clin. Endocrinol. Metab., July 1, 2009; 94(7): 2571 - 2578. [Abstract] [Full Text] [PDF]

				Y. Yao, Y. Zhao, L. P. Smith, M. Watson, and V. Nair Novel MicroRNAs (miRNAs) Encoded by Herpesvirus of Turkeys: Evidence of miRNA Evolution by Duplication J. Virol., July 1, 2009; 83(13): 6969 - 6973. [Abstract] [Full Text] [PDF]

				L. Cheng and K. C. Keiler Correct Timing of dnaA Transcription and Initiation of DNA Replication Requires trans Translation J. Bacteriol., July 1, 2009; 191(13): 4268 - 4275. [Abstract] [Full Text] [PDF]

				L. Zargarian, I. Kanevsky, A. Bazzi, J. Boynard, F. Chaminade, P. Fosse, and O. Mauffret Structural and dynamic characterization of the upper part of the HIV-1 cTAR DNA hairpin Nucleic Acids Res., July 1, 2009; 37(12): 4043 - 4054. [Abstract] [Full Text] [PDF]

				H. Wang, R. Li, and Y. Hu The Alternative Noncoding Exons 1 of Aromatase (Cyp19) Gene Modulate Gene Expression in a Posttranscriptional Manner Endocrinology, July 1, 2009; 150(7): 3301 - 3307. [Abstract] [Full Text] [PDF]

				J. Waldispuhl, S. Devadas, B. Berger, and P. Clote RNAmutants: a web server to explore the mutational landscape of RNA secondary structures Nucleic Acids Res., July 1, 2009; 37(suppl_2): W281 - W286. [Abstract] [Full Text] [PDF]

				M. R. Sharma, T. M. Booth, L. Simpson, D. A. Maslov, and R. K. Agrawal Structure of a mitochondrial ribosome with minimal RNA PNAS, June 16, 2009; 106(24): 9637 - 9642. [Abstract] [Full Text] [PDF]

				J.-Y. Lin, S.-R. Shih, M. Pan, C. Li, C.-F. Lue, V. Stollar, and M.-L. Li hnRNP A1 Interacts with the 5' Untranslated Regions of Enterovirus 71 and Sindbis Virus RNA and Is Required for Viral Replication J. Virol., June 15, 2009; 83(12): 6106 - 6114. [Abstract] [Full Text] [PDF]

				K. M. Gustin, B.-J. Guan, A. Dziduszko, and D. A. Brian Bovine Coronavirus Nonstructural Protein 1 (p28) Is an RNA Binding Protein That Binds Terminal Genomic cis-Replication Elements J. Virol., June 15, 2009; 83(12): 6087 - 6097. [Abstract] [Full Text] [PDF]

				K. C. Pang, M. E. Dinger, T. R. Mercer, L. Malquori, S. M. Grimmond, W. Chen, and J. S. Mattick Genome-Wide Identification of Long Noncoding RNAs in CD8+ T Cells J. Immunol., June 15, 2009; 182(12): 7738 - 7748. [Abstract] [Full Text] [PDF]

				D. J. Scanlan, M. Ostrowski, S. Mazard, A. Dufresne, L. Garczarek, W. R. Hess, A. F. Post, M. Hagemann, I. Paulsen, and F. Partensky Ecological Genomics of Marine Picocyanobacteria Microbiol. Mol. Biol. Rev., June 1, 2009; 73(2): 249 - 299. [Abstract] [Full Text] [PDF]

				S. A. Yatsenko, E. K. Brundage, E. K. Roney, S. W. Cheung, A. C. Chinault, and J. R. Lupski Molecular mechanisms for subtelomeric rearrangements associated with the 9q34.3 microdeletion syndrome Hum. Mol. Genet., June 1, 2009; 18(11): 1924 - 1936. [Abstract] [Full Text] [PDF]

				F. B. Stabell, N. J. Tourasse, and A.-B. Kolsto A conserved 3' extension in unusual group II introns is important for efficient second-step splicing Nucleic Acids Res., June 1, 2009; 37(10): 3202 - 3214. [Abstract] [Full Text] [PDF]

				M. Nongkhlaw, P. Dutta, J. W. Hockensmith, S. S. Komath, and R. Muthuswami Elucidating the mechanism of DNA-dependent ATP hydrolysis mediated by DNA-dependent ATPase A, a member of the SWI2/SNF2 protein family Nucleic Acids Res., June 1, 2009; 37(10): 3332 - 3341. [Abstract] [Full Text] [PDF]

				R. Hoshina and N. Imamura Phylogenetically Close Group I Introns with Different Positions among Paramecium bursaria Photobionts Imply a Primitive Stage of Intron Diversification Mol. Biol. Evol., June 1, 2009; 26(6): 1309 - 1319. [Abstract] [Full Text] [PDF]

				R. Cheng, C. Miao, Q. Gong, Y. Gu, X. Lu, F. Han, and W. Yu Specific gene silencing by artificial trans-encoded small noncoding RNAs in bacteria Nucleic Acids Res., May 27, 2009; (2009) gkp447v1. [Abstract] [Full Text] [PDF]

				Z. Wang, K. Treder, and W. A. Miller Structure of a Viral Cap-independent Translation Element That Functions via High Affinity Binding to the eIF4E Subunit of eIF4F J. Biol. Chem., May 22, 2009; 284(21): 14189 - 14202. [Abstract] [Full Text] [PDF]

				M. Wels, R. S. Bongers, J. Boekhorst, D. Molenaar, M. Sturme, W. M. de Vos, R. J. Siezen, and M. Kleerebezem Large Intergenic Cruciform-Like Supermotifs in the Lactobacillus plantarum Genome J. Bacteriol., May 15, 2009; 191(10): 3420 - 3423. [Abstract] [Full Text] [PDF]

				D. Stoddart, A. J. Heron, E. Mikhailova, G. Maglia, and H. Bayley From the Cover: Single-nucleotide discrimination in immobilized DNA oligonucleotides with a biological nanopore PNAS, May 12, 2009; 106(19): 7702 - 7707. [Abstract] [Full Text] [PDF]

				M. Wieland, M. Gfell, and J. S. Hartig Expanded hammerhead ribozymes containing addressable three-way junctions RNA, May 1, 2009; 15(5): 968 - 976. [Abstract] [Full Text] [PDF]

				L. Childs, Z. Nikoloski, P. May, and D. Walther Identification and classification of ncRNA molecules using graph properties Nucleic Acids Res., May 1, 2009; 37(9): e66 - e66. [Abstract] [Full Text] [PDF]

				H.-M. Chen and S.-H. Wu Mining small RNA sequencing data: a new approach to identify small nucleolar RNAs in Arabidopsis Nucleic Acids Res., May 1, 2009; 37(9): e69 - e69. [Abstract] [Full Text] [PDF]

				W. B. Langdon, G. J. G. Upton, and A. P. Harrison Probes containing runs of guanines provide insights into the biophysics and bioinformatics of Affymetrix GeneChips Brief Bioinform, May 1, 2009; 10(3): 259 - 277. [Abstract] [Full Text] [PDF]

				D. Xu and H. Shi Composite RNA aptamers as functional mimics of proteins Nucleic Acids Res., May 1, 2009; 37(9): e71 - e71. [Abstract] [Full Text] [PDF]

				M. Caporale, F. Arnaud, M. Mura, M. Golder, C. Murgia, and M. Palmarini The Signal Peptide of a Simple Retrovirus Envelope Functions as a Posttranscriptional Regulator of Viral Gene Expression J. Virol., May 1, 2009; 83(9): 4591 - 4604. [Abstract] [Full Text] [PDF]

				G. C. Saunders, S. Cawthraw, S. J. Mountjoy, A. C. Tout, A. R. Sayers, J. Hope, and O. Windl Ovine PRNP untranslated region and promoter haplotype diversity J. Gen. Virol., May 1, 2009; 90(5): 1289 - 1293. [Abstract] [Full Text] [PDF]

				A. Cimerman, D. Pacifico, P. Salar, C. Marzachi, and X. Foissac Striking Diversity of vmp1, a Variable Gene Encoding a Putative Membrane Protein of the Stolbur Phytoplasma Appl. Envir. Microbiol., May 1, 2009; 75(9): 2951 - 2957. [Abstract] [Full Text] [PDF]

				A. San Millan, F. Depardieu, S. Godreuil, and P. Courvalin VanB-Type Enterococcus faecium Clinical Isolate Successively Inducibly Resistant to, Dependent on, and Constitutively Resistant to Vancomycin Antimicrob. Agents Chemother., May 1, 2009; 53(5): 1974 - 1982. [Abstract] [Full Text] [PDF]

				K. T. Gagnon, S.-Y. Ju, M. B. Goshe, E. S. Maxwell, and S. Franzen A role for hydrophobicity in a Diels-Alder reaction catalyzed by pyridyl-modified RNA Nucleic Acids Res., May 1, 2009; 37(9): 3074 - 3082. [Abstract] [Full Text] [PDF]

				B. Szarzynska, L. Sobkowiak, B. D. Pant, S. Balazadeh, W.-R. Scheible, B. Mueller-Roeber, A. Jarmolowski, and Z. Szweykowska-Kulinska Gene structures and processing of Arabidopsis thaliana HYL1-dependent pri-miRNAs Nucleic Acids Res., May 1, 2009; 37(9): 3083 - 3093. [Abstract] [Full Text] [PDF]

				S. L. Fishman, S. H. Factor, C. Balestrieri, X. Fan, A. M. DiBisceglie, S. M. Desai, G. Benson, and A. D. Branch Mutations in the Hepatitis C Virus core Gene Are Associated with Advanced Liver Disease and Hepatocellular Carcinoma Clin. Cancer Res., May 1, 2009; 15(9): 3205 - 3213. [Abstract] [Full Text] [PDF]

This Article Abstract Print PDF (328K) Supplementary Material Alert me when this article is cited Alert me if a correction is posted Services Email this article to a friend Similar articles in this journal Similar articles in ISI Web of Science Similar articles in PubMed Alert me to new issues of the journal Add to My Personal Archive Download to citation manager Search for citing articles in: ISI Web of Science (2284) Request Permissions Commercial Re-use Guidelines for Open Access NAR Content Google Scholar Articles by Zuker, M. Search for Related Content PubMed PubMed Citation Articles by Zuker, M. Social Bookmarking          What's this?

Online ISSN 1362-4962 - Print ISSN 0305-1048

Copyright © 2009 Oxford University Press

Oxford Journals Oxford University Press

• Site Map
• Privacy Policy
• Frequently Asked Questions

Other Oxford University Press sites: Oxford University Press Oxford Journals China Oxford Journals Japan American National Biography Booksellers' Information Service Children's Fiction and Poetry Children's Reference Corporate & Special Sales Dictionaries Dictionary of National Biography Digital Reference English Language Teaching Higher Education Textbooks Humanities International Education Unit Law Medicine Music Online Products Oxford English Dictionary Reference Rights and Permissions Science School Books Social Sciences Very Short Introductions World's Classics