Module r_group_asl

This function counts the number of heavy atoms in the structure.

Parameters:

• st (structure.Structure) - Structure

Returns: int

number of heavy atoms

   This function returns important rotatable bonds. Important rotatable bonds are defined when both sides of the
bond have
   either rings or heavy atoms greater than atom_cutoff.

   @type            st: L{structure.Structure}
   @param           st: Structure
   @type   atom_cutoff: C{int}
   @param  atom_cutoff: # of heavy atoms
   @return            : generator for atom pairs in important rotatable bonds
   

Find reactant and product ct from the input, and return touple of reference ligand and mutated ligand structures

Parameters:

• cts - input structures

Returns:

ref_ct, mut_ct

Parameters:

• b0 - first atom index for bond
• b1 - second atom index for bond, this connects to the fragment returned

Returns:

list of atom indice for all atom connected to b1 when bond b0-b1 is removed

Returns:

list of attachment bond pairs

two cts need to have i_fep_subst codes defined ?99-?98 if one of the fragment the linker atom connects contains less than 6 heavy atoms, treat this fragment as functional group

Returns:

list of atom indices Find the core atoms from which a mutation is attached

dummy atoms here includes mutated ring atoms and linker atoms

Returns:

list of dummpy atom indices

combine fused rings from SSSR set

Returns:

list of lists of atoms

given anchor atom and connection graph of a molecule, find the rotable bonds that is cloest (in terms of number of bonds to the anchor atom)

Parameters:

• graph (networkx graph) - molecular graph
• anchor_index - anchor atom index
• rotable_bonds (list of touples) - list of rotable bonds

remove the extra fragments connected to the current hot region.

Parameters:

• ct - structure
• graph (networkx graph) - molecular graph
• anchor_index - anchor atom index
• hot_r (list of integers) - hot region found by finding closes rotable bonds from anchor atom
• rotable_bonds (list of touples) - list of rotable bonds

find rotable bond that is closest to the anchor atom in terms of number of bonds
return the fragment, original rotable bond, bond used to find the fragment,
and list of rotable bonds from which large fragment is cut off (extra_rot).
The last two items are used to find the matching fragment in the other molecule
@return: (hot_r, save_rot_bond, (rot_atom0, rot_atom1), extra_rot)
     hot_r: hot region found
     save_rot_bond: original rotable bond used to find the fragment
     (rot_atom0, rot_atom1): rotable bond in the same atom order passed to get_fragment
     extra_rot: list of extra_rotable bond directly connected to the fragment;
                the fragments attached by these bonds are removed already

given attachment pairs, find the rotable fragments containing the attachment bond return fragments and the rotable bond associated The bigger one of the fragment pair is used to determine which rotabond to use.

Parameters:

• source_ct - source structure
• dest_ct - dest structure
• source_att (touple of integers) - source attachment point
• dest_att (touple of integers) - dest attachment point
• source_rotable_bonds - source rotable bond list
• dest_rotable_bonds - dest rotable bond list

Returns:

(source_hot_frag, dest_hot_frag, (s_rot0, s_rot1), (d_rot0, d_rot1)). The order of atoms in the rotable bonds returned is in accord to the way the fragments are generated.

extend the hot region by moving to the next rotable bond

Parameters:

• total_heavy - total heavy atoms in the hot region so far
• source_ct - source structure
• dest_ct - dest structure
• att_bonds - all attachment bonds for mutations
• source_rot_bonds - rotable bond list of reactant
• dest_rot_bonds - rotable bond list of product
• source_frag_rot_bonds_visite - dictionary of rotable bonds already visited for attchment bonds in reactant
• dest_frag_rot_bonds_visite - dictionary of rotable bonds already visited for attchment bonds in product

This is called whenever the i_fep_mapping is not defined.

Parameters:

• cts - input cts

Returns:

none

@param cts: input structures, s_fep_fragname should be defined for sourec and dest cts
@type cts: Lists of structure
@param max_rest_atoms: maximum allowed number of heavy atoms in hot region, soft limit
@param min_h_atom: minimum heavy atom moved by a rotable bond
@param exclusion_lists: touple of sets of core atoms (hydrogen included)
@param get_fragment: callable to compute fragment given (atom_index0, atom_index1, ct)
@param get_attachment_bonds: callable to compute attathment points given
                             two cts with i_fep_mapping defined (ct0, ct1)
@return: number of heavy atoms assigned

Set the hotregions in the cts according to mutated fragments and maximum allowed number
of atoms.  Core atoms can be exclude via exclusion_lists.

Steps to set REST region
1.  Include all the changed atoms (?99, ?97, ?98)
    If there are too many atoms, stop.  Linker atoms are ignored
    (linker connects two fragments with more than 6 heavy atoms in each of the fragments)
2.  If a rotable bond is also an attachment point, include all atoms in the fragment.
3.  For each changed fragment pair (attachment point), try to include all the atoms up
    to the closest rotatble bond.  How far away the rotable bond is is measure by
    the minimum number of bonds separating the fragment and the rotable bonds.
    Fused rings are excluded.  If the fragment is a ring, the maximum allowed atoms may
    be six more than the allowed number.
4.  Repeat 3, if maximum number of atoms are not reached and it is still possible to add
    more atoms.

n_set

Value:

set_hotregion(cts, max_rest_atoms= MAX_NUM_REST_ATOM, min_h_atom= N_HA TOM)