Official government website icon

bit	Xplor-NIH	VMD-XPLOR

Pydoc: package pasd

pasd

index

Stubs for future PASD NOE assignment support.

Package Contents

formatPipp
formatSTAR
formatSparky

Functions

covalentNeighbors(atom, range): return list of atoms within range bonds of argument atom

findUnassignedAtoms(shiftList, flexibleRegion='', tclOutput=False, simulation=None, expectedNonassigned='\n (name o*) or \n (name c) or \n (name s*) or\n (resn thr and name hg1) or \n (resn ser and name hg) or \n (resn lys and (name hz* nz)) or \n (resn tyr and (name hh cg cz)) or \n (resn arg and (name he hh* cz nh* ne)) or\n (resn gln and name cd) or \n (resn glu and name cd) or \n (resn asn and name cg) or \n (resn asp and name cg) or\n (resn phe and name cg) or \n (resn trp and (name cg cd2 ce2)) or \n (resn his and name cg) or\n (resn pro and name n) or\n (name ht*)'): Given a list of shift entries, find atoms with no assignment

expectedNonassigned - default atoms which we don't expect to see
                      (exchangable protons, oxygens, sulfurs, or
                      heavyatoms with no attached protons)
flexibleRegion      - selection specifying atoms which we don't expect
                      to see due to flexibility.
tclOutput           - if True, return an informational string,
                      else return a named tuple with string member
                      message a list of atoms missingAtoms.

hasAssignmentWithinBondedNeighborhood(peak, range): Return True if the specified peak contains any assignments containing
a from/to pair separated by range bonds or fewer.

nefShifts(nef, name=None, verbose=True, tclOutput=False): Read chemical shift data from the global NEF cif.Cif
object, which must have been previously populated by calling readNEF().

Returns a list of two-membered lists:
[chemical shift value, selection string]

peakAssignLikelihoodsFromStructs(structureData, peaks, violCutoff=0.5, sim=None): #
# Given a list of <filename, atomPosArray> pairs of converged structures
# and a list of peaks,
# use the coords to determine the likelihood of each peakAssignment
#
# Could be replaced with something that would calculate previous likelihood
# of any pair of ShiftAssignments, whether they're used in a PeakAssignment
# or not.
#

pipeShifts(pipe, verbose=True, tclOutput=False): Read chemical shift data from NMRPipe table format.

pipe is a string specifying either the path of the table or, directly, its
contents.

removeBondedPeaks(pot, bondRange=2): Remove all peaks from the < pasdPot>.PASDPot argument pot which
contain any assignments with one or more from/to pair separated by
bondRange or few bonds.

structureMetrics(coordsByFilename, peaks, pasdPots, violCutoff=0.5, completenessWeight=0.0, inverseBound=4.0, inverseMethylCorrection=0.0, verbose=True, sim=None): Given a list of <filename, atomPosArray> pairs from grabPDBfiles,
check each one's violations of a given set of peaks
and return a list of the named tuple StructureMetrics for each structure.
This list is sorted by the structure's score.

updateUser(mess): Output the message to stdout if pasd.updateUser_enabled is set to
True.

xeasyPeaks(filename, noePot, fromProtonColumnName, toProtonColumnName, heightColumnName, fromHeavyatomColumnName=None, toHeavyatomColumnName=None, idColumnName=None, name=None, verbose=True, tclOutput=True): Read XEASY format peaks into pasdPot.PASDPot

arguments:

  filename                - name of XEASY peaks file
  noePot                  - a PASDPot instance
  fromProtonColumnName    - string label for from- Proton
  toProtonColumnName      - string label for to- Proton
  fromHeavyatomColumnName - string label for from- heavyatom
  toHeavyatomColumnName   - string label for to- heavyatom
  heightColumnName             name                    - prefix to use for peak names. In PASD,
                            these must be different for each peaklist.
  verbose                 - whether to produce verbose output
  tclOutput               - If True, return a valid TCL list string

Not robustly tested.

xeasyShifts(filename, verbose=True, tclOutput=True, simulation=None, segmentName=None): Read XEASY format chemical shifts. Not robustly tested.

Data
		aveExp = 6 distanceBins = [(1.8, 2.7, 0.2), (1.8, 3.3, 0.5), (1.8, 5.0, 0.8), (1.8, 6.0, 1.0)] longRangeResidCutoff = 6 nMono = 1 updateUser_enabled = False