| Trees | Indices | Help |
|
|---|
|
|
This script will break up a set of input molecules into fragments
based on some simple rules, similar to those described in the original
RECAP paper. Run with -h to see the basic rules.
If -murcko option is specified, then Murco Assembly scaffolds are generated
instead of fragments. This is done by removing terminal from each "branch"
until reacing a ring.
If -recap option is specified, then rules are as follows:
I. Only break the bonds matching the RECAP SMARTS patterns.
1 = Amide
2 = Ester
3 = Amine
4 = Urea
5 = Ether
6 = Olefin
7 = Quarternary nitrogen,
8 = Aromatic nitrogen - aliphatic carbon
9 = Lactam nitrogen - aliphatic carbon
10 = Aromatic carbon - aromatic carbon
11 = Sulphonamide
II. Do not break ring bonds
III. Do not break bonds that would yield one of the following fragments:
hydrogen, methane, ethane, prone, butane.
NOTE: -recap_use can be used to specify which RECAP rules to use. The value
must be a list of comma-separated numbers (no spaces).
Example usage::
from schrodinger.structutils import createfragments
frag_creator = createfragments.FragmentCreator(atoms=options.atoms,
bonds=options.bonds, smarts=options.smarts,
carbon_hetero=options.carbon_hetero,
maxatoms=options.maxatoms,
removedup=options.removedup, murcko=options.murcko,
recap=options.recap, recap_use=recap_list,
complete=options.complete,
add_h=options.add_h, verbose=True)
for struct in my_structures:
# fragments is a list of Structure objects
fragments = frag_creator.fragmentStructure(struct)
Copyright Schrodinger, LLC. All rights reserved.
|
|||
|
_Options A class to hold a list of options. |
|||
|
FragmentCreator This class will break up input structures into fragments based on some simple rules. |
|||
|
_StructureFragmenter Class to create fragments from a given set of _Options and a given schrodinger.structure.Structure. |
|||
|
|||
| bool |
|
||
| set of tuples |
|
||
|
|||
logger = log.get_output_logger(__file__)
|
|||
ATTACHPROP =
|
|||
RECAP_PROP_PREFIX =
|
|||
RECAP_PATTERNS =
|
|||
TERMINAL_PATTERNS =
|
|||
__package__ =
|
|||
|
|||
Checks coordinates of each atom in the specified structures and returns true if they are the same
|
Find all bonds that match the RECAP rules described in J. Chem. Inf. Comput. Sci. 1998, 38, 511-522 Will return a list of [(atom1, atom2, rule_num)].
|
|
|||
RECAP_PATTERNS
|
TERMINAL_PATTERNS
|
| Trees | Indices | Help |
|
|---|
| Generated by Epydoc 3.0.1 on Wed Aug 3 07:59:19 2016 | http://epydoc.sourceforge.net |