Python: module sardcPotTools
tools to aid in setup/analysis of dipolar coupling potential term sardcPot
this module provides functions to simplify the creation, manipulation and
analysis of sardcPot.SARDCPot potential terms.
Functions |
| |
- Rfactor(pot, selection='all', normalize=False)
- R-factor (in percent) for an infinite number of randomly distributed
vectors.
The formula is based on Eq. 3 from Clore+Garrett, JACS 121, 9008 (1999):
R-factor = 100 [1/sqrt(2)] sqrt{ sum[(Di_calc - Di_obs)^2] / sum[Di_obs^2] }
where Di_calc and Di_obs are the ith calculated and observed RDC values,
respectively.
The selection argument can be used to choose a subset of restraints whose
atoms lie in selection.
If normalize is set to True, all values of Di_calc and Di_obs are divided by
sardcPot.SARDCPot_Restraint's Dmax in the above equation.
- addRestraints(rdc, restraints, useSign=True, useDist=False)
- Add XPLOR XDIPO or SARDC style restraints from a string.
One can specify the useSign argument to specify whether the RDC sign will be
significant in evaluating the energy. The useDist argument specifies whether
the 1/r^3 RDC dependence will be taken from the coordinates.
- analyze(potList)
- perform analysis of SARDCPot terms and return nicely formatted summary
- avecScale(pot)
- return pot.rdc.avectorScale()
- chi2(pot, selection='all')
- Compute the Chi^2 value for the specified restraints.
The optional selection argument can be used to choose a subset of
restraints whose atoms lie in selection.
- create_SARDCPot(name, file=0, useSign=True, useDist=False, restraints='', domainSel='not PSEUDO', tensor=None)
- Create an sardcPot.SARDCPot with given name, the filename of an RDC
assignment table and/or a string of assignments.
The file argument can optionally be a sequence of filenames.
domainSel is an atomSel.AtomSel which specifies atoms in an aligning
subunit.
If tensor is specified, it should be a <saTensor>.SATensor object.
One can specify the useSign argument to specify whether the RDC sign will be
significant in evaluating the energy. The useDist argument specifies whether
the 1/r^3 RDC dependence will be taken from the coordinates.
- makeTable(rdc)
- Return the assignment table (a string) corresponding to the
restraints associated with the specified sardcPot.SARDCPot.
- readRestraint(string, useSign, useDist)
- rmsd(pot, selection='all')
- Compute the RMSD value for the specified restraints.
The optional selection argument can be used to choose a subset of
restraints whose atoms lie in selection.
- saupeMatrix(pot, eIndex=None)
- given a SARDCPot term, return the associated Saupe matrix. If eIndex
is specified, return the Saupe matrix associated with the specified
ensemble member. NOT YET TESTED!
- writeNEF(rdc, name)
- Return a formatted NEF record for the restraints in the given SARDCPot
object as a string with the specified saveframe name.
|
Data |
| |
distances = {'CA_C': 1.52, 'CA_CB': 1.52, 'C_N': 1.33, 'N_HN': 1.023} |