Proper RMSD calculation between molecules using the Kuhn-Munkres Hungarian algorithm.
Project description
RMSD
Determining the Root Mean Square Deviation (RMSD) between the Maximum Common Substructures (MCS) of two molecules.
The Kuhn-Munkres Hungarian algorithm allows for a fast match of atoms based on:
- atomic symbol
- Sybyl atom types
- pharmacophoric types (i.e. H-bond donors and acceptors or charges)
Installation
pip install munkres-rmsd
Example
from munkres_rmsd import CalcLigRMSD, AtomType
from munkres_rmsd.RMSD import get_example_molecules
# First, load 3D poses of molecules
mol1, mol2 = get_example_molecules()
# Then compute the RMSD of the best atomic match
rmsd = CalcLigRMSD(mol1, mol2)
print(rmsd) # 10.76150...
Let's use Sybyl atom types to match atoms between the two molecules instead of the default, using atomic elements.
# Then compute the RMSD of the best atomic match
rmsd = CalcLigRMSD(mol1, mol2, AtomType.Sybyl)
print(rmsd) # 11.59752...
Should you prefer pharmacophore types (i.e. H-bond donors & acceptors, charges and others):
# Then compute the RMSD of the best atomic match
rmsd = CalcLigRMSD(mol1, mol2, AtomType.Pharmacophore)
print(rmsd) # 9.49120...
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
Close
Hashes for munkres_rmsd-0.0.1.post1.dev0.tar.gz
Algorithm | Hash digest | |
---|---|---|
SHA256 | a68513d9d5683a44e61a6ff6b1d09e853c28916cafd6d75def8a7c604bcc2f24 |
|
MD5 | 34ef2680f3248731dd7c6c72e50d7fa3 |
|
BLAKE2b-256 | 27b4f824483188e729f477680301d85111f14f9b6c06ec2c9417e15af32d1dc2 |
Close
Hashes for munkres_rmsd-0.0.1.post1.dev0-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 054066a3d928a140caa688016f47b5b8c48b8a020fde46c0fa19caecabaddbcf |
|
MD5 | 8d68f8b0206ea6f9e1a60a57cd55f094 |
|
BLAKE2b-256 | 7af25bb93b1222c68ef6124cd0285c7139db25eaf16d94a29eaa3c1aa42bdd96 |