MolTopolParser, abbreviated from Molecular Topology Parser, is a lightweight, open-source Python package designed for accessing, validating, and processing files or data from various molecular simulation and modeling software. It simplifies interfacing between different simulation software and is pivotal in building computational pipelines for molecular simulations.
-
Python 3.6 or newer
-
pip package manager
To install MolTopolParser, run the following command:
pip install moltopolparser
MolTopolParser requires the following Python packages:
These dependencies are automatically installed during the MolTopolParser installation process.
-
Modular Design: Structured around specific file formats from various simulation software, enhancing maintainability and scalability.
-
Pydantic Dataclasses: Ensures robust data validation and type checking across the package.
-
Parsing as Class Methods: Facilitates easy access to and manipulation of simulation files through well-defined class methods.
-
Exposed Modules and Dataclasses: Enables advanced users to extend functionality or integrate with other tools.
MolTopolParser is designed for developers to manage hierarchical data structures commonly used in molecular simulations. It provides a robust foundation for building advanced simulation tools and pipelines.
-
Accessing and manipulating simulation files beyond mere coordinate data.
-
Interfacing between different molecular simulation software.
-
Building and managing computational pipelines.
-
Simple parsing of structural data.
-
Basic analysis of simulation outputs. (Consider using Biopython , MDAanalysis etc.)
import moltopolparser as mtp
# Validate base data directly
atom_data = {
"id": 10, "atom_type": "C", "resnr": 100,
"residu": "UREA", "atom": "C1", "cgnr": 1, "charge": -0.683,
}
atom = mtp.gmx.MolTopAtom(**atom_data)
print(atom.charge)
# Parsing a single .gro file
gro_file = '../tests/data/gmx/two_water.gro'
gro = mtp.gmx.GroFile.parser(gro_file)
print(gro.box_size)
# Parsing and accessing a top file
top_file = '../tests/data/gmx/membrane-martini-charmmgui/system.top'
sys_top = mtp.gmx.Topology.parser(top_file)
print(sys_top.system)
sys_top.pull_forcefield()
print(sys_top.forcefield.atomtypes[0]) # Output: name='P5' at_num=None mass=72.0 charge=0.0 ptype='A' sigma=0.0 epsilon=0.0
sys_top.pull_molecule_topologies()
print(sys_top.molecule_topologies[0])
We encourage contributions from the community! Whether you are fixing bugs, proposing new features, or improving the documentation, your help is welcome.
-
For Beginners Submit issue requests for bugs or feature suggestions.
-
For Developers Read the concepts.md and follow our Developer's Guidence to understand how to set up your environment and contribute code.
For details on the latest changes and version history, please refer to the full Change Log on GitHub.
MolTopolParser is released under the Apache License. See the LICENSE file for more details.