Scoring of quaternary structures (QS). The QS scoring is according to the paper
by `Bertoni et al. <https://dx.doi.org/10.1038/s41598-017-09654-8>`_.
.. note ::
Requirements for use:
- A default :class:`compound library <ost.conop.CompoundLib>` must be defined
and accessible via :func:`~ost.conop.GetDefaultLib`. This is set by default
when executing scripts with ``ost``. Otherwise, you must set this with
:func:`~ost.conop.SetDefaultLib`.
- ClustalW must be installed (unless you provide chain mappings)
- Python modules `numpy` and `scipy` must be installed and available
(e.g. use ``pip install scipy numpy``)
def ost.mol.alg.qsscoring.GetContacts |
( |
|
entity, |
|
|
|
calpha_only, |
|
|
|
dist_thr = 12.0 |
|
) |
| |
Get inter-chain contacts of a macromolecular entity.
Contacts are pairs of residues within a given distance belonging to different
chains. They are stored once per pair and include the CA/CB-CA/CB distance.
:param entity: An entity to check connectivity for.
:type entity: :class:`~ost.mol.EntityHandle` or :class:`~ost.mol.EntityView`
:param calpha_only: If True, we only consider CA-CA distances. Else, we use CB
unless the residue is a GLY.
:type calpha_only: :class:`bool`
:param dist_thr: Maximal CA/CB-CA/CB distance to be considered in contact.
:type dist_thr: :class:`float`
:return: A connectivity dictionary. A pair of residues with chain names
*ch_name1* & *ch_name2* (*ch_name1* < *ch_name2*), residue numbers
*res_num1* & *res_num2* and distance *dist* (<= *dist_thr*) are
stored as *result[ch_name1][ch_name2][res_num1][res_num2]* = *dist*.
:rtype: :class:`dict`
Definition at line 870 of file qsscoring.py.