schrodinger.pipeline.stages.pull module

A stage for extracting a subset of ligands/compounds.

The Pull stage identifies a subset of compounds from one set of ligand structures, and extracts all variants of those compounds from a second set of ligand structures.

Copyright Schrodinger, LLC. All rights reserved.

class schrodinger.pipeline.stages.pull.PullStage(*args, **kwargs)

Bases: schrodinger.pipeline.stage.Stage

Stage for extracting a subset of ligands/compounds (and their variants) from a second set of ligand files, given the number or percent to keep from a first set of ligand files. The first set must be ordered; i.e., the ligands to keep appear earliest in the set. Compounds in the second set that don’t appear in the first set are ignored.

The keywords specific to this stage are…

UNIQUEFIELD The ligand property that identifies a compound.
Multiple ligands with the same UNIQUEFIELD value are considered variants of the compound. The default is the structure title (‘s_m_title’).
NUM_TO_KEEP The number of compounds from the first set to
extract from the second set. The compounds kept are those that appear earliest in the first set.
PERCENT_TO_KEEP The percentage of compounds from the first set to
extract from the second set. The compounds kept are those that appear earliest in the first set. Ignored if NUM_TO_KEEP is used.
KEEP_CHARGES If set to True, then after the ligand is pulled from
the original file, the partial charges are set to what they were in the second file (For QPLD). NOTE: If set, only ONE ligand with specified field should be in the first file set. NOTE: solvation_charge property is assummed to be equal to the partial_charge.
CHARGE_PROPERTY If KEEP_CHARGES is True, the charges are taken from
this atom-level property. Default is ‘r_m_charge1’. These charges from the first set are placed in the ‘r_m_charge1’ (partial_charge) and ‘r_m_charge2’ (solvation_charge) properties of the pulled set, not the CHARGE_PROPERTY in the pulled set.

If neither NUM_TO_KEEP or PERCENT_TO_KEEP are present, all compounds in the first set are extracted from the second set. In all cases, every variant of a kept compound from the second set is extracted.

This stage takes one ordered input structure set that identifies the ligands to extract, a second input structure set from which those ligands are extracted, and it generates one output structure file set, each containing about 100,000 structures.

Issues: Keeping the first N compounds may not choose compounds from Glide results in the proper manner. If the variants of a compound are chemically distinct, it’s appropriate to choose the top compounds based on GlideScore. However, if the variants differ only in their conformations (such as when saving multiple poses per ligand), Emodel (not GlideScore) should be used to determine the representative variant of the compound for GlideScore comparison with other compounds. Hopefully, this doesn’t come up often, because it would be very tricky to reflect the proper ordering in the input files if there are both Ligprep-type variants and conformation-variants.

__init__(*args, **kwargs)

See class docstring.

operate()

Extract a subset of the first set compounds from the second set of ligand input structure files. Makes a list of compounds (identified by UNIQUEFIELD) from the first set, pares down the list according to the NUM_TO_KEEP or PERCENT_TO_KEEP setting, and extract all variants of those compounds from the second set. Raises a RuntimeError if any input file cannot be read, if there is a problem accessing the UNIQUEFIELD property of any ligand, if there is a problem writing any output file, or if no ligands are extracted.

pullCompounds()
JobDJOptions()

Returns a dictionary of options to pass to JobDJ: hosts, local, max_retries, default_max_retries, verbosity

__contains__(key)
__len__()
addExpectedInput(position, type, required=True)

A stage can accept one or more pipeio input objects. Use this method to specify the type of input object that is expected at each position.

position - an integer starting at 1. type - structures/grids/etc. required - whether this input always needs to be specified

addExpectedOutput(position, type, always=True)

A stage can return one or more pipeio objects. Use this method to specify the type of object that will be returned and whether or not it will always be produced by the stage.

position - an integer starting at 1. type - structures/grids/etc. always - whether this output is always produced

addOutputFile(filename)

Adds the specified file to the stage’s job control record. File must be specified as local (not absolute) path.

checkFile(file, error='File does not exist:')

Raise exception if specified file does not exist. The message that is printed can be specified.

checkFiles(files, error='File does not exist')

Raise expetion if any file does not exist.

checkInputs()

OVERWRITE: Return False if something is wrong with the input files or the parameter, otherwise return True.

checkParameters()

OVERWRITE: Make sure that all parameters are valid.

checkProducts()

Raises RuntimeError if any of the required products are not installed or the version installed is less that minimum required version. It is possible to override this method. See ligprep.py for example.

clear() → None. Remove all items from D.
copy()
debug(text)

Print a debug line to the log file

dump()

This method dumps all the variables of the Stage to a restart file. Run it every time an important step is performed.

error(text)

Print an error line to the log file

exit(text='')

Print an error line to the log file and exit with code 1

classmethod fromkeys(iterable, value=None)
genFileName(extension=None, filenum=None, start=None, end=None)

Generate a file name to be used by the stage. Returns string:

"<full-stagename>-<filenum><extension>"
"<full-stagename>-<start>_<end><extension>"
"<full-stagename><extension>"
"<full-stagename>", etc.

Depending on given options.

genOutputFileName(position, extension='', filenum=None, start=None, end=None)

Generate a file name to be used by the stage when writing files for the output position <position>. Returns strings:

"<full-varname>-<filenum><extension>"
"<full-varname>-<start>_<end><extension>"
"<full-varname><extension>"
"<full-varname>", etc.

Depending on given options.

get(k[, d]) → D[k] if k in D, else d. d defaults to None.
getAdjustedNJobs(total_mol, min_job_size, max_job_size)

Returns the desired number of subjobs, and adjusts it for the the specified min & max job sizes if the user specified ADJUST option. If the number of desired jobs was specified by the user, the number of available cpus is used or 10, whichever is smaller. Specify the number of input ligands and the smallest and largest desired job sizes (Generally job lengths of 1 minute & 24 hours).

getCleanupRequested()

Stages should clean up after themselves if this returns True

getHostList()

Returns a list of hosts to run the subjobs on. localhost:1 may be in the list as well. Ideally, pass the output to JobDJ. Format: [ (host1,ncpus), (host2,ncpus) ] NOTE: whether or not to run with -LOCAL should be up to the stage. Pass this value to JobDJ.

getHostStr()

Just like getHostList() but instead of returning a list, returns a host string to be passed to the -HOST argument.

getInput(position)

Use in operate() to get the input object for specified position. Returns None if invalid position is specified.

getInputNames()

Return a dictionary of variable name of the inputs at each position. Key:position, value:name

getJobDJ(**kwargs)

Returns a pre-set JobDJ instance for the stage to use. It already has it’s hosts, local, max_retries, max_failures, default_max_retries, and verbosity set.

getLocal()

Whether subjobs should be run with -LOCAL. Pass this value to JobDJ

getMaxRetries()

Return the number of max restarts to use. If -max_retries is specified, returns that value; otherwise if SCHRODINGER_MAX_RETRIES is defined, returns that value; otherwise returns default of 2. Pass this value to JobDJ.

getNCpus()

Returns the total number of processors specified in the host string. For queued hosts with no CPU# specification, 10 is added.

getNJobs()

Returns the requested target number of subjobs, and whether or not to adjust that number if it is unreasonable.

If -NJOBS was not specified, the # of CPUs or 10 is returned (whichever is smaller).

Used by Glide DockingStage and _adjustNJobs()

getName()

Return stagename (jobname of the stage)

getOutput(position)

Returns the output IO object of the stage at specified position. Use this method after running the stage to get its output objects

getOutputName(position)

Return the output name for specified position

getOutputNames()

Return a list of output names for each position (dict)

getRuntimePath(filename)

Return the runtime-path of a file that user specified Prints an error and exits if file does not exist.

getStageDirectory()

Return the directory in which the stage is running

getVerbosity()

Return verbosity of thos stage (for JobDJ)

hasCompleted()

Returns True if this stage’s operate() exited successfully.

hasStarted()

Returns True if this stage has started.

info(text)

Print an info line to the log file

items() → a set-like object providing a view on D's items
iterInputs()

Iterate through input objects: (position, obj)

keys() → a set-like object providing a view on D's keys
lognoret(*args)

Prints specified objects to the stage log file. No EOF return

mainProduct()

If a stage has a main product associated with it, the stage should overwrite this method with a method that returns the product string. For example, the LigPrepStage.mainProduct() will return “ligprep” Used by Pipeline.

outputRequested(position)

Returns True if the user requested optional output at <position>

pop(k[, d]) → v, remove specified key and return the corresponding value.

If key is not found, d is returned if given, otherwise KeyError is raised.

popitem() → (k, v), remove and return some (key, value) pair

as a 2-tuple; but raise KeyError if D is empty.

reportParameters(fh=None)

Print the value of each keyword for this stage to the stream specified as <fh>. Used by Pipeline

requiredProduct(product)

Specify a product that is required for this stage to run; optionally minimum version.

Example: product=”mmshare”

requiredProductRuntime(product)

Similar to requiredProduct() but can be used to specify required products at runtime. For example, ConvertStage doesn’t know what products are required for conversion until runtime. Raises RuntimeError if product is not installed.

run(idle_function=None, restart_file=None, verbosity=None, logfh=None)

Run the stage.

idle_function - function to call when idle

restart_file - file to periodically dump this instance to

verbosity - there are three verbosity levels: “quiet”, “normal”, and “verbose”
“quiet” - only warnings and errors are printed “normal” - stage progress is printed - default “verbose” - additional debugging info is printed

logfh - where to send the loggin output

runJobDJ(jobdj)
setInput(position, name=None, obj=None)

Specify an input to use for this stage. position - input specified is for this position name - Variable name of this IO object obj - the IO object

This method is called by Pipeline.

setJobDJOptions(jobdj)

Use this method to adjust the specified queue.JobDJ instance to the VSW settings.

setJobOptions(subjob_hosts=None, njobs=None, adjust=None, force=None, subjob_local=None, max_retries=None, cleanup=None)

Tell this stage how to run the subjobs

Parameters:
  • subjob_hosts – list of hosts to run subjobs on
  • njobs – number of subjobs to generate. None means determine automatically.
  • adjust – whether to adjust njobs such that job size is within limits
  • force – whether to continue with job if subjobs fail
  • max_retries – number of times to attempt to restart a subjob If not specified, use SCHRODINGER_MAX_RETRIES or 2.
  • cleanup – whether to delete intermediate files
setMainProduct(product)

Specify which product this stage is part of. Will determine which host the subjobs are run on.

setOutput(position, obj)

Use this method at the end of operate() to set the output.

setOutputName(position, varname)

Is called by Pipeline when starting the stage. Tell the stage what name to save each output under.

setdefault(k[, d]) → D.get(k,d), also set D[k]=d if k not in D
update([E, ]**F) → None. Update D from mapping/iterable E and F.

If E present and has a .keys() method, does: for k in E: D[k] = E[k] If E present and lacks .keys() method, does: for (k, v) in E: D[k] = v In either case, this is followed by: for k, v in F.items(): D[k] = v

updateJobdj(jobdj)

Gets called periodically in order to update JobDJ’s hosts. Will ask Pipeline for CPUS when needed, and will tell Pipeline when they are no longer needed.

validateValues(preserve_errors=False)

Validates the stored keywords. This is done by converting <self> to a ConfigObj instance, and calling validate() on it. The validated keywords are then updated back to <self>. This is done as part of Ev:87429

values() → an object providing a view on D's values
warning(text)

Print a warning line to the log file