RELIC toolkit calculates the electrostatic interaction throughout the molecular dynamics trajectory. The interactions between user-defined regions are calculated on a per-residue bases. Relic can read AMBER, CHARMM, DESRES, GAMESS, GROMACS, GROMOS96, GSD, HOOMD and LAMMPS trajectories. Tracking the evolution of electrostatic interactions over the simulation time allows for identification of the residue-residue interactions that are the driving force behind the observed process.
When publishing results obtained using RELIC please cite this publication:
(to be added)
If you do not have setuptools installed, or have a version earlier than 40.0, install it with:
$ pip install --upgrade setuptools
To install the RELIC toolkit, download the zip from GitHub, extract it to a directory of your choice and run:
$ python setup.py install
Alternatively, you can install it with pip:
$ pip install git+https://github.com/mtomin/relictoolkit.git
- Linux
- Python 2.7 and up
- setuptools>=40.0
- numpy
- matplotlib
- six
- requests
- gsd
- MDAnalysis
- decorator
- plotly
- pytz
- traitlets
- validators
- future
- configparser
- duecredit
The installation can be tested by running:
$ python setup.py test
which checks if all the core modules are working properly
The setup.py
script creates the executable file in the /bin directory and the program can be run with:
$ relictoolkit
In addition to the built-in GUI, calculation and plotting scripts can be run as standalone scripts:
$ python relictoolkit_calculate path/to/config/file
$ python relictoolkit_plot path/to/config/file
In this case, it is necessary to manually create the config.ini or config_plot.ini files containing the calculation/plotting parameters.
The examples folder contains a short simulation of a peptide comprised of 15 amino acid residues (HRERFLADKSAPLCG). This tutorial will go through calculating interactions between the first (residues 1-7) and second (residues 8-15) half of the protein.
In the Calculate tab navigate to the topology and trajectory files located in the relictoolkit/examples folder. Set the location of the output file and make sure the filetype "Amber binary trajectory (.ncdf)" is selected. Mask1 and Mask2 define the regions of the studied system. In our case, Mask1 contains residues 1 to 7, whereas Mask2 contains residues 8-15. RELIC uses the MDAnalysis selection syntax, so to select appropriate residues enter:
Mask1: resid 1:7
Mask2: resid 8:15
If left empty, distance cutoff defaults to 5 Å, internal dielectric constant defaults to 4, stride defaults to 1000 and timestep to 1 fs. In this example, we will use a stride of 500 to generate enough data for a decent plot later. A timestep used in the simulations was 2 fs. The interface should now look like this:
Enter the desired number of cores and hit run. The calculation time should not exceed 2 minutes with a single core. The process can be tracked using the logfile:
$ tail relic_logfile.log
Upon selecting the output file obtained in the previous part of the tutorial, the settings pane appears. The default settings can be used (plot against time, include all frames and all residues). Currently RELIC includes two options for data plotting, matplotlib (non-interactive) and plotly (interactive).
After hitting run, the graph is created and automatically opened. The graph can be further customized using Plotly's Chart studio or exported to png file using the action menu bar in the upper right corner.
Matplotlib | Plotly |
---|---|
It can be observed that residues Glu3 and Lys9 contribute most to the interactions between the first and the second half of the peptide. Those interactions are beginning to form at around 1.9 ns, while at that time residues 2 and 4 formed destabilizing interactions, mostly with Lys9.
In case of dependency issues, try creating a clean environment with anaconda or miniconda:
$ conda create myenv python --no-default-packages
$ source activate myenv
And retry the installation.
If you are having issues with the "traditional" setuptools installation, try installing via pip.
If the dependency issues persist, install the requirements with
$ pip install -r requirements.txt
RELIC is distributed under the CC-BY license.