analyze nucleic acid 3D structures and MD trajectories
Project description
Introduction
Barnaba is a tool for analyzing RNA three-dimensional structures and simulations. Barnaba uses MDtraj to read/write topology and trajectory files, as such it supports several formats including pdb, xtc, trr, dcd, binpos, netcdf, mdcrd, prmtop, and more. Barnaba has been developed by Sandro Bottaro with the crucial help of Giovanni Bussi, Giovanni Pinamonti, Sabine Rei{ss}er and Wouter Boomsma.
This is what you can do with Barnaba:
Calculate eRMSD [1]
Calculate RMSD after optimal alignment
Search for single/double stranded RNA motifs in the PDB database or in simulations [1]
Annotate PDB structures and trajectories with the Leontis-Westhof classification
Produce dynamic secondary structure figures in SVG format
Cluster nucleic acids structures using the eRMSD as a metric distance
Calculate elastic network models for nucleic acids and nucleic acids/protein complexes [2]
Calculate backbone and pucker torsion angles in a PDB structure or trajectory
Back-calculate 3J scalar couplings from PDB structure or trajectory
Score three-dimensional structures using eSCORE [1]
For bugs, questions or comments contact Sandro at sandro dot bottaro (guesswhat) gmail dot com
If you use Barnaba in your work, please cite the following paper:
@article{bottaro2019barnaba,
title={Barnaba: software for analysis of nucleic acid structures and trajectories},
author={Bottaro, Sandro and Bussi, Giovanni and Pinamonti, Giovanni and Rei{\ss}er, Sabine and Boomsma, Wouter and Lindorff-Larsen, Kresten},
journal={RNA},
volume={25},
number={2},
pages={219--231},
year={2019},
publisher={Cold Spring Harbor Lab}
}
The manuscript is also available on biorXiv here: https://www.biorxiv.org/content/10.1101/345678v3
Requirements
- Barnaba requires:
Python >= 3.6
Numpy
Scipy
Mdtraj 1.9
future
Barnaba requires mdtraj (http://mdtraj.org/) for manipulating structures and trajectories. To perform cluster analysis, scikit-learn is required too.
Required packages can be installed using pip, e.g.:
pip install mdtraj
Installation
You can obtain the latest tagged version of barnaba using pip:
pip install barnaba
If you prefer to manage your dependencies with conda you can use:
conda -c conda-forge install barnaba
On MacOS, you can install the same tagged version using the python distributed with MacPorts:
sudo port install py36-barnaba
Just replace 36 with the python version that you prefer to use.
Alternatively, you can find the most recent version of barnaba on Github:
git clone git://github.com/srnas/barnaba.git
then move to the barnaba directory and run the command
pip install -e .
Usage
Barnaba can be either used as a Python library or as a commandline tool. A number of Notebook examples can be found in the examples directory. The notebooks for conducting the analyses and producing the figures in the manuscript can be found in the folder manuscript_figures .
Alternatively, the command-line interface can be found in the bin directory. Here’s a minimal how-to
minimal help: barnaba –help
Calculate the ERMSD between structures
barnaba ERMSD –ref ../test/data/sample1.pdb –pdb ../test/data/sample2.pdb
trajectories can be provided as well, by specifying a topology file
barnaba ERMSD –ref ../test/data/sample1.pdb –top ../test/data/sample1.pdb –trj ../test/data/samples.xtc
other accepted options are shown in a function-specific help
barnaba ERMSD –help
Calculate the RMSD between structures
barnaba RMSD –ref ../test/data/sample1.pdb –pdb ../test/data/sample2.pdb –dump
Find single stranded motif
barnaba SS_MOTIF –query ../test/data/GNRA.pdb –pdb ../test/data/1S72.pdb
Find double stranded motif. l1 and l2 are the lengths of the two strands
barnaba DS_MOTIF –query ../test/data/SARCIN.pdb –pdb ../test/data/1S72.pdb –l1 8 –l2 7
Annotate structures/trajectories according to the Leontis/Westhof classification.
barnaba ANNOTATE –pdb ../test/data/SARCIN.pdb
Produce dynamic secondary-structure figures. It requires as input the files .pairing and .stacking produced with the ANNOTATE command.
barnaba SEC_STRUCTURE –ann outfile.ANNOTATE.stacking.out outfile.ANNOTATE.pairing.out
Calculate backbone/sugar/pseudorotation angles
barnaba TORSION –pdb ../test/data/GNRA.pdb –backbone –sugar –pucker
Calculate J-couplings
barnaba JCOUPLING –pdb ../test/data/sample1.pdb
Calculate elastic network models for RNA and predict SHAPE reactivity. NB: only works with PDB.
barnaba ENM –pdb ../test/data/GNRA.pdb –shape
Calculate relative positions between bases R_ij ang G vectors for pairs within ellipsoidal cutoff
barnaba DUMP –pdb ../test/data/GNRA.pdb –dumpG –dumpR
Extract fragments from structures with a given sequence. NB: only works with PDB.
barnaba SNIPPET –pdb ../test/data/1S72.pdb –seq NNGNRANN
Calculate ESCORE
barnaba ESCORE –ff ../test/data/1S72.pdb –pdb ../test/data/sample1.pdb
References
- [1] Bottaro, Sandro, Francesco Di Palma, and Giovanni Bussi.
“The role of nucleobase interactions in RNA structure and dynamics.” Nucleic acids research 42.21 (2014): 13306-13314.
- [2] Pinamonti, Giovanni, et al.
“Elastic network models for RNA: a comparative assessment with molecular dynamics and SHAPE experiments.” Nucleic acids research 43.15 (2015): 7260-7269.
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
