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...
Finally to allow matching any atom whatever its type:
# Then compute the RMSD of the best atomic match
rmsd = CalcLigRMSD(mol1, mol2, AtomType.Skeleton)
print(rmsd) # 8.96288...
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.post3.dev0.tar.gz
Algorithm | Hash digest | |
---|---|---|
SHA256 | 33bb99c6e219ebf4c4f6a0a0d99cfb3933ddc8d346c70903bf1f55d5a169f0c4 |
|
MD5 | 8678eb8258f260721c35f8cdc958cfd2 |
|
BLAKE2b-256 | 15ed39ee6856aa56d03379975235a709d0ad85ec0b1c294324b1bce4b75626d6 |
Close
Hashes for munkres_rmsd-0.0.1.post3.dev0-py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 7ba4cf4c46c39fab85c7daec60ddfb5a87ad8931eb1b8ad64f729eb5e43e499d |
|
MD5 | 41a618e0fc1f849b4d1d1574d8d55145 |
|
BLAKE2b-256 | 985873ccaa265dc4056a4646d051294ffd85db8af1f8223fd237dd2787331550 |