GenTopo: A python toolkit for generating gromacs topology template.
Source
Code link: GitHub - masrul/GenTopo
Description
For non-standard Force-field (which cannot be invoked through pdb2gmx), this tool can be used to generate gromacs topology template. It uses slightly different approach than x2top. It requires a PDB file of the molecule, which has two extra entries after the 80th column/character respectively atomTypes and charges. An example PDB file can be found in Example/test.pdb. It is recommended to have connectivity information in the PDB file, but not mandatory. If connectivity is missing, then the code will generate connectivity information using the heuristic of Van der Waals radius of atoms. In addition to internal coordinates, it also creates all required bondTypes, angleTypes, and dihedralTypes. GenTopo is independent of particular force filed, it simply creates template. User needs to provide required parameters once template is created. An example topology template created by GenTopo can be found in Example/topol.top.
# File name: Example/test.py
# python test.py
from GenTopo.Coord import PDBobj
from GenTopo.Graph import MolGraph
from GenTopo.GMXTopo import Topo
# load a PDB file
mol = PDBobj("test.pdb")
# Create a Molecular Graph for generating bond, angles, dihedral etc.
graph = MolGraph(mol, guessImpropers=True)
# Now create gromacs topology by passing mol and graph
gmx = Topo(mol, graph)
# write topology file
gmx.write("topol.top")
Features
-
PDB to gromacs template
-
Generates internal coordinates
- Bonds
- Angles
- Dihedrals
- Pairs (1-4)
- Improper dihedrals (for aromatic molecule)
-
Generates types
- Bonds types
- Angle types
- Dihedral types
-
Support periodic molecule
Installation:
git clone https://github.com/masrul/GenTopo
cd GenTopo
pip install .
Contact Person: Masrul Huda
Email:: mmh568@msstate.edu