Generate coarse-grained molecular dynamics models from atomistic trajectories.
Project description
PyCGTOOL
Generate coarse-grained molecular dynamics models from atomistic trajectories.
The aim of this project is to provide a tool to aid in parametrising coarse-grained (CG) molecular mechanics models. PyCGTOOL generates coarse-grained models from atomistic simulation trajectories using a user-provided mapping. Equilibrium values and force constants of bonded terms are calculated by Boltzmann Inversion of bond distributions collected from the input trajectory.
Alternatively map-only mode (behaving similarly to MARTINIZE) may be used to generate initial coordinates to use with existing CG topologies such as the MARTINI lipid models. For instance, a pre-equilibrated atomistic membrane may be used to create starting coordinates for a MARTINI membrane simulation.
PyCGTOOL makes it easy to test multiple variations in mapping and bond topology by making simple changes to the config files.
If you find PyCGTOOL useful, please cite our JCIM paper (https://doi.org/10.1021/acs.jcim.7b00096) and the code itself (https://doi.org/10.5281/zenodo.598143).
@article{Graham2017,
author = {James A. Graham and Jonathan W. Essex and Syma Khalid},
doi = {10.1021/acs.jcim.7b00096},
issn = {1549-9596},
issue = {4},
journal = {Journal of Chemical Information and Modeling},
month = {4},
pages = {650-656},
title = {PyCGTOOL: Automated Generation of Coarse-Grained Molecular Dynamics Models from Atomistic Trajectories},
volume = {57},
url = {https://pubs.acs.org/doi/10.1021/acs.jcim.7b00096},
year = {2017},
}
Install
PyCGTOOL requires Python 3.6 or higher and may be installed using pip:
pip install pycgtool
MDTraj on macOS
On some versions macOS, with some versions of the Clang compiler, MDTraj may fail to load GROMACS XTC simulation trajectories. If you encounter this issue, make sure you have the latest version of MDTraj.
For more information see MDTraj/#1572.
Usage
Input to PyCGTOOL is an atomistic simulation trajectory in the form of a topology (e.g. PDB, GRO, etc.) and a trajectory file (e.g. XTC, DCD, etc.), along with two custom files: MAP and BND. These files provide the atomistic-to-CG mapping and bonded topology respectively.
Example files are present in the test/data directory. The format of these files is described in the full documentation.
For more information, see the tutorial. It is important to perform validation of any new parameter set - a brief example is present at the end of the tutorial.
For a full list of options, see the documentation or use:
pycgtool -h
Generate a Model
To generate a CG model from an atomistic simulation:
pycgtool <topology file> <trajectory file> -m <MAP file> -b <BND file>
Map Only
To use PyCGTOOL to convert a set of atomistic simulation coordinates to CG coordinates:
pycgtool <topology file> -m <MAP file>
Or to convert a complete simulation trajectory:
pycgtool <topology file> <trajectory file> -m <MAP file>
Maintainers
James Graham (@jag1g13)
Contributing
If you experience problems using PyCGTOOL or wish to see a new feature added please open an issue or submit a PR.
To help develop PyCGTOOL, you can create a fork of this repository, clone your fork and install PyCGTOOL using:
poetry install
This will install PyCGTOOL in editable mode (similar to pip install -e .
) along with all the necessary runtime and development dependencies.
The Makefile at the root of the repository contains targets for running unit and integration tests (make test
) and linting (make lint
).
License
GPL-3.0 © James Graham, University of Southampton
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Hashes for pycgtool-2.0.0b1-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 78647b74e30a2f0d25d8ecb57b7d4d9fdf85408744eae4e1fc4fd9118904eae8 |
|
MD5 | b9c57db98e21da1f8031cedf6b67f87c |
|
BLAKE2b-256 | 3439d185c03ec38bcdcfe71aacde247b661c5d1cd4b2fea1556f8e115b3d7c01 |