A package to caculate residue-residue contacts from MD trajectories
Project description
trajcontacts
A Package for Calculating Residue-Residue Contacts from MD Trajectories
Currently, there are two modes for contact determination:
-
Normal Contact Analysis (norm): The presence of a contact between two residues in macromolecular systems, such as proteins and nucleic acids, is established when the minimum inter-residue heavy atom pair distance falls within a specified cutoff (4-5 Å) for a majority of the simulation time (e.g., 75% or more of snapshots). This program can be utilized to extract residue contacts from long MD simulation trajectories as well as PDB structures.
-
Continuous Contact Analysis (cont): As described in Botello-Smith et al., Journal of Chemical Theory and Computation 15.4 (2019): 2116-2126 and Westerlund et al., The Journal of Chemical Physics 153.14 (2020).
This is a command-line program for Linux systems and relies on the MDTraj simulation analysis package (https://www.mdtraj.org/). It uses mdtraj.load()
for trajectory/topology inputs, and inter-atomic distances are calculated using mdtraj.distance()
. For efficiency and to avoid memory-related errors, save the trajectories without unnecessary atoms (solvent/ions or hydrogen atoms) and divide the trajectory into multiple chunks, then write their names (or pathnames) in a file (for example, name that file 'traj.dat') - see below.
Installation
You can install trajcontacts by running:
pip install trajcontacts
After installation, you can access detailed options by typing:
trajcontacts -h
For a simple usage example, running the following command:
trajcontacts -p topologyfile -f file_containing_names/paths_to_coordinate_files -n 10
will execute the program on 10 processors and extract contacts with a 4.5 Ångstrom cutoff. Both 'topologyfile' (specified in the command line) and 'coordinatefile' (in the file containing names of coordinate files) should be specified as the same for extracting contacts from a PDB file.
Example
trajcontacts -p 3sn6.pdb -f traj.dat -n 10
where 'traj.dat' is a file with a single line containing '3sn6.pdb'.
How to Cite:
Please cite the following research when using trajcontacts:
Madhu, Midhun K., Kunal Shewani, and Rajesh K. Murarka. "Biased Signaling in Mutated Variants of β2-Adrenergic Receptor: Insights from Molecular Dynamics Simulations." bioRxiv (2023): 2023-09.
Copyright:
Copyright belongs to the Computational Biophysics and Soft Matter Group, IISER Bhopal. Visit our website: https://home.iiserb.ac.in/~rkm/.
Project details
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distributions
Hashes for trajcontacts-1.0.0-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 2d9b0a8813307926ebce44d718c0116abafe01a444d17cc702adcdc39977da26 |
|
MD5 | 1d42b802017fe6089313d4f30811da80 |
|
BLAKE2b-256 | 038fe3f2fe078734b95925dfa7bb4f81d677187f45b9f4dc995834cc4688c8e1 |
Hashes for trajcontacts-1.0.0-py2.py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 2a9527e2ce8d7321efe3b1f0c22c027bd16e61e174e86dee18a6e41fbfad2e02 |
|
MD5 | d73576d4b8be869f46896e4c0e21f40a |
|
BLAKE2b-256 | c51d478f2ebce74594d26e82565d08a338faec346ddd9073b301e95f14c57a42 |