The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.

TmCalculator

v1.0.9

Genome-wide nucleic acid melting temperature (Tm) profiling and multi-omics integration. Results are returned as GRanges objects, so Tm can be used directly as a quantitative genomic feature alongside ATAC-seq, RNA-seq, ChIP-seq and other assays.

1. install

install.packages("TmCalculator")

install dev version from github

pak::pkg_install("JunhuiLi1017/TmCalculator@dev")

2. usage and examples

Please see the vignetts for the details.

library(TmCalculator)

seqs <- to_genomic_ranges("AAAATTTTTTTCCCCCCCCCCCCCCGGGGGGGGGGGGTGTGCGCTGC")
tm_calculate(seqs, method = "tm_nn", nn_table = "DNA_NN_SantaLucia_2004", Na = 50)

3. thermodynamic parameter sets

Twenty-seven nearest-neighbor parameter sets are available, in two families.

Reference-salt sets were fitted at a single reference sodium concentration. Other conditions are reached through the salt_method correction formulas.

Duplex Sets
DNA/DNA DNA_NN_Breslauer_1986, DNA_NN_Sugimoto_1996, DNA_NN_Allawi_1998, DNA_NN_SantaLucia_2004 (default)
RNA/RNA RNA_NN_Freier_1986, RNA_NN_Xia_1998, RNA_NN_Chen_2012
RNA/DNA RNA_DNA_NN_Sugimoto_1995

Condition-specific sets were fitted directly at the sodium concentration shown, by melting-temperature optimization. They are intended to replace salt correction rather than be corrected. When the requested Na matches the concentration a set was fitted at, salt correction is skipped automatically; when it does not, the correction is applied with a warning.

Duplex Sets Fitted at
DNA/DNA DNA_NN_Weber_2015 1020 mM
DNA/DNA DNA_NN_Weber_OW04_69 / _119 / _220 / _621 / _1020 69–1020 mM
RNA/RNA RNA_NN_Weber_VIF_71 / _121 / _221 / _621 / _1021 71–1021 mM
RNA/RNA RNA_NN_Weber_FIF_71 / _121 / _221 / _621 / _1021 71–1021 mM
RNA/DNA RNA_DNA_NN_Weber_2019_FT, RNA_DNA_NN_Weber_2019_VH 1000 mM
RNA/DNA RNA_DNA_NN_Weber_2019_LS 100 mM

For RNA, the VIF (variable initiation factors) sets gave better cross-validation than FIF. For RNA/DNA hybrids at high salt, ..._FT was the best-performing set in the source study.

# Fitted at 100 mM, so no salt correction is applied on top of it
res <- tm_calculate(seqs, method = "tm_nn",
                    nn_table = "RNA_DNA_NN_Weber_2019_LS", Na = 100)
res$options[["Salt correction applied"]]                    # FALSE
res$options[["Parameter set fitted at [Na+] (mM)"]]         # 100

Pick the set whose fitted salt is closest to your experimental condition rather than correcting a distant one. See ?tm_nn for the full list and citations.

4. launch an R shiny application

using R function TmCalculatorShiny::TmCalculator_shiny()

5. citation

If you use the melting-temperature-optimized parameter sets, please also cite the source studies:

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.