Package schrodinger :: Package application :: Module ska
[hide private]
[frames] | no frames]

Module ska

Core routines for wrapping SKA backend.

This module provides direct access to the functionality used by the top-level ska program and structalign utility. The usual keyword-value options in the top-level input file are supplied as a python dict. Wrappers are supplied to provide PDB filenames as input as well as maestro CTs which are transformed in place.

Copyright Schrodinger, LLC. All rights reserved.

Classes [hide private]
  _SkaData
Data structure for named backend return values...
Functions [hide private]
 
force_rmdir(dir)
Recursive 'rm -r' using fileutils.force_remove
 
transform_structure(st, u, x)
perform in-place transformation
 
standardize(old_ct, rename, reorder)
Convert non-standard residue names to standard forms where possible (e.g.
list of _SkaData
pairwise_align_ct(query, templist, keywords=None, log=None, debug=False, save_props=False, std_res=False, reorder=False)
Wrapper for a series of pairwise ska jobs with CTs as input.
_SkaData
multiple_align_ct(query, templist, keywords=None, log=None, debug=False, save_props=False, std_res=False, reorder=False)
Wrapper for a single multiple-alignment ska job with CTs as input.
list of _SkaData
pairwise_align_pdb(query, templist, keywords=None, log=None, debug=False, std_res=False, reorder=False)
Wrapper for a series of pairwise ska jobs with PDB files as input.
_SkaData
multiple_align_pdb(query, templist, keywords=None, log=None, debug=False, std_res=False, reorder=False)
Wrapper for a single multiple-alignment ska job with PDB files as input.
_SkaData
run_align(structs, keywords, log=None, debug=False)
Core driver routine, equivalent to $SCHRODINGER/ska with all top-level input options specified in keywords.
Variables [hide private]
  _true_val = re.compile(r'^(1|TRUE|YES)$', re.I)
  _atom_rec = re.compile(r'^((ATOM |HETATM).{24})([ -.\d]{8})([...
  _winlen_vals = [2, 6, 9, 14, 17]
  _trolltop = None
hash(x)
Function Details [hide private]

standardize(old_ct, rename, reorder)

 

Convert non-standard residue names to standard forms where possible (e.g. HID -> HIS)

pairwise_align_ct(query, templist, keywords=None, log=None, debug=False, save_props=False, std_res=False, reorder=False)

 

Wrapper for a series of pairwise ska jobs with CTs as input.

Parameters:
  • query ((seqname, ct) tuple) - fixed reference structure
  • templist (list of (seqname, ct) tuples) - one or more structures to be aligned
  • keywords (dict) - top-level keyword-value input options (other than QUERY_FILE and TEMPLATE)
  • log (logging.logger) - active logger for diagnostic messages
  • debug (bool) - debug flag passed to backend
  • save_props (bool) - True if output data included as CT properties
  • std_res (bool) - True if residue names are translated to standard forms
  • reorder (bool) - True if residues need to be reordered by connectivity
Returns: list of _SkaData
list of alignments corresponding to templist

multiple_align_ct(query, templist, keywords=None, log=None, debug=False, save_props=False, std_res=False, reorder=False)

 

Wrapper for a single multiple-alignment ska job with CTs as input.

Arguments are the same as in pairwise_align_ct above.

Returns: _SkaData
a single multiple alignment

pairwise_align_pdb(query, templist, keywords=None, log=None, debug=False, std_res=False, reorder=False)

 

Wrapper for a series of pairwise ska jobs with PDB files as input.

Parameters:
  • query (string) - fixed reference PDB file
  • templist (list of strings) - one or more PDB files to be aligned
  • keywords (dict) - top-level keyword-value input options (other than QUERY_FILE and TEMPLATE)
  • log (logging.logger) - active logger for diagnostic messages
  • debug (bool) - debug flag passed to backend
  • std_res (bool) - True if residue names are translated to standard forms
  • reorder (bool) - True if residues need to be reordered by connectivity
Returns: list of _SkaData
list of alignments corresponding to templist

multiple_align_pdb(query, templist, keywords=None, log=None, debug=False, std_res=False, reorder=False)

 

Wrapper for a single multiple-alignment ska job with PDB files as input.

Arguments are the same as in pairwise_align_pdb above.

Returns: _SkaData
a single multiple alignment

run_align(structs, keywords, log=None, debug=False)

 

Core driver routine, equivalent to $SCHRODINGER/ska with all top-level input options specified in keywords. Returns alignment, scores, and rotation matrices as well as backend stdout.

Parameters:
  • structs (list of Structures) - structures to align (first in list is query)
  • keywords (dict) - top-level keyword-value input options
  • log (logging.logger) - active logger for diagnostic messages
  • debug (bool) - debug flag passed to backend
Returns: _SkaData
a single multiple alignment

Variables Details [hide private]

_atom_rec

Value:
re.compile(r'^((ATOM  |HETATM).{24})([ -.\d]{8})([ -.\d]{8})([ -.\d]{8\
})(.*)$')