20211221 - 可视化RNA二级结构-主要是arc diagram

Julse

已于 2023-03-20 21:13:51 修改

阅读量3k

点赞数 5

分类专栏： RNA 文章标签： python

于 2021-12-21 22:34:49 首次发布

本文链接：https://blog.csdn.net/Julse/article/details/122073923

版权

RNA 专栏收录该内容

2 篇文章 0 订阅

订阅专栏

文章目录

Structure Editor
成功 VARNA （基于Java的jar）
失败的尝试 - RChie/R4RNA
- 在R中运行Rchie
- 在python中使用Rchie
其他包 arc-diagrams

预期的效果图DSSR介绍文档中的图
在这里插入图片描述

实际实现的图

缺点：base上面的数字看不清楚

Structure Editor

官方网站
https://www.urmc.rochester.edu/rna/
下载链接（需要科学上网才能打开）
https://rochester.app.box.com/s/lnjtj4t8vvi18y25x4b9qydz2k5vsvhe
在这里插入图片描述

成功 VARNA （基于Java的jar）

官方下载VARNA 的jar包
安装 java running，即 jre
然后选这个文件夹下/bin/java.exe 打开jar上述jar包即可显示窗口

GUI上只能显示这一种图，命令行不指定输出目录，即 -o 参数，则也可以打开GUI显示，但是指定算法的命令不起作用
查看相关命令
四种二级结构展示形式，直接输出到文件打开显示
-resolution 110 尝试了多种图像的分辨率，绘制6ZVK，最多设置成了110，成功运行，换成其他RNA又会报错，即使110，放大图像到极限，也看不清楚base上面的数字

在这里插入图片描述

参数algorithm 设置为line算法就可以了

java -cp VARNAv3-93.jar fr.orsay.lri.varna.applications.VARNAcmd -i D:\data\dssr_output\6ZVK\dssr-2ndstrs.dbn -o 6ZVK.png -algorithm line -resolution 50

里面的三种方式都尝试画了一下，以及默认的方式

line
circular

在这里插入图片描述

neview
radiate [默认]

失败的尝试 - RChie/R4RNA

网站里面的图很漂亮
该R包的官方使用手册 https://e-rna.org/r-chie/files/R4RNA_manual_2021.pdf
在这里插入图片描述

首先是 rna-tools 中提到了很多有用的rna工具，其中包括这个绘制二级结构的工具，安装方法
安装R语言， RChie 工具

RChie 下载之后是这样

2.1 根据官方给的命令在R中安装Rchie – 失败

# https://rna-tools.readthedocs.io/en/latest/tools.html#secondary-structure-format-conversion
path_to_file = 'D:\\Program Files\\R\\Rlib\\r-chie\\R4RNA'
install.packages("path_to_file", repos = NULL, type="source")
# Install the optparse and RColorBrewer
install.packages('optparse')
install.packages('RColorBrewer')

在这里插入图片描述
2.2 根据另一个博客先安装BiocManager 再直接安装RChie – 成功

install.packages("BiocManager")
if (!requireNamespace("BiocManager", quietly = TRUE))
    install.packages("BiocManager")

BiocManager::install("R4RNA")
browseVignettes("R4RNA")

在这里插入图片描述