schrodinger.application.bindingsite.intfield module

Logic to compute “interaction fields” between (protein) atoms and “probe” particles.

schrodinger.application.bindingsite.intfield.get_hbond_direction(atom)

Returns direction of “ideal” hydrogen bond that would be formed by the given atom. It is not always right, e.g. for O in C-O-H. This needs to be refined or proven irrelevant.

schrodinger.application.bindingsite.intfield.get_interaction_sites(st, atom_indices=None, probes=None, logger=None)

Identifies interaction sites among the protein atoms according to their atom types.

Parameters
  • st (schrodinger.structure.Structure) – Structure.

  • atom_indices (iterable over int) – Iterable over the contributing atom indices.

  • probes (container of str) – Probes of interest.

Returns

List of atoms that interact with the requested probes. Individual atom contributions are given by tuples that hold the -1/0/1 integers (see INTERACTIONS) associated with the corresponding probe.

Return type

list(tuple(schrodinger.structure._StructureAtom, tuple(int)))

class schrodinger.application.bindingsite.intfield.Field(st, atom_indices=None, probes=None, alpha=1.0, r_cut=4.0, a_cut=60.0, logger=None)

Bases: object

Handles computation of the “interaction potentials” generated by the “interaction sites” (protein atoms) acting on “probes”.

__init__(st, atom_indices=None, probes=None, alpha=1.0, r_cut=4.0, a_cut=60.0, logger=None)
Parameters
  • st (schrodinger.structure.Structure) – Structure.

  • atom_indices (iterable over int) – Iterable over the contributing atom indices.

  • probes (container of str) – Probes of interest.

  • alpha (float) – Interaction range (length scale of exponential decay).

  • r_cut (float) – Ignore contributions from atoms further than r_cut from a probe.

  • a_cut (float) – Ignore hydrogen bond interactions for angles exceeding a_cut.

nearest_atom_distance(pos)

Returns distance to the nearest atom that contributes to the potentials.