IMP Reference Guide
develop.49f026e7cd,2024/12/04
The Integrative Modeling Platform
|
Setup cross-link distance restraints at atomic level The "atomic" aspect is that it models the particle uncertainty with a Gaussian. More...
Inherits RestraintBase.
Setup cross-link distance restraints at atomic level The "atomic" aspect is that it models the particle uncertainty with a Gaussian.
The noise in the data and the structural uncertainty of cross-linked amino-acids is inferred using Bayes' theory of probability
Definition at line 515 of file crosslinking.py.
Public Member Functions | |
def | __init__ |
Constructor. More... | |
def | create_restraints_for_rmf |
create dummy harmonic restraints for each XL but don't add to model Makes it easy to see each contribution to each XL in RMF More... | |
def | get_best_stats |
return the probability, best distance, two coords, and possibly the psi for each xl More... | |
def | get_particles_to_sample |
Get the particles to be sampled by the IMP.pmi.sampler object. More... | |
def | load_nuisances_from_stat_file |
Read a stat file and load all the sigmas. More... | |
def | plot_violations |
Create CMM files, one for each state, of all cross-links. More... | |
def IMP.pmi.restraints.crosslinking.AtomicCrossLinkMSRestraint.__init__ | ( | self, | |
root_hier, | |||
xldb, | |||
atom_type = 'NZ' , |
|||
length = 10.0 , |
|||
slope = 0.01 , |
|||
nstates = None , |
|||
label = None , |
|||
nuisances_are_optimized = True , |
|||
sigma_init = 5.0 , |
|||
psi_init = 0.01 , |
|||
one_psi = True , |
|||
filelabel = None , |
|||
weight = 1.0 |
|||
) |
Constructor.
Automatically creates one "sigma" per cross-linked residue and one "psis" per pair. Other nuisance options are available.
root_hier | The root hierarchy on which you'll do selection |
xldb | CrossLinkDataBase object |
atom_type | Can either be "NZ" or "CA" |
length | The XL linker length |
slope | Linear term to add to the restraint function to help when far away |
nstates | The number of states to model. Defaults to the number of states in root. |
label | The output label for the restraint |
nuisances_are_optimized | Whether to optimize nuisances |
sigma_init | The initial value for all the sigmas |
psi_init | The initial value for all the psis |
one_psi | Use a single psi for all restraints (if False, creates one per XL) |
filelabel | automatically generated file containing missing/included/excluded cross-links will be labeled using this text |
weight | Weight of restraint |
Definition at line 540 of file crosslinking.py.
def IMP.pmi.restraints.crosslinking.AtomicCrossLinkMSRestraint.create_restraints_for_rmf | ( | self | ) |
create dummy harmonic restraints for each XL but don't add to model Makes it easy to see each contribution to each XL in RMF
Definition at line 754 of file crosslinking.py.
def IMP.pmi.restraints.crosslinking.AtomicCrossLinkMSRestraint.get_best_stats | ( | self, | |
limit_to_state = None , |
|||
limit_to_chains = None , |
|||
exclude_chains = '' , |
|||
use_CA = False |
|||
) |
return the probability, best distance, two coords, and possibly the psi for each xl
limit_to_state | Only examine contributions from one state |
limit_to_chains | Returns the particles for certain "easy to visualize" chains |
exclude_chains | Even if you limit, don't let one end be in this list. Only works if you also limit chains |
use_CA | Limit to CA particles |
Definition at line 973 of file crosslinking.py.
def IMP.pmi.restraints.crosslinking.AtomicCrossLinkMSRestraint.get_particles_to_sample | ( | self | ) |
Get the particles to be sampled by the IMP.pmi.sampler object.
Definition at line 781 of file crosslinking.py.
def IMP.pmi.restraints.crosslinking.AtomicCrossLinkMSRestraint.load_nuisances_from_stat_file | ( | self, | |
in_fn, | |||
nframe | |||
) |
Read a stat file and load all the sigmas.
This is potentially quite stupid. It's also a hack since the sigmas should be stored in the RMF file. Also, requires one sigma and one psi for ALL XLs.
Definition at line 809 of file crosslinking.py.
def IMP.pmi.restraints.crosslinking.AtomicCrossLinkMSRestraint.plot_violations | ( | self, | |
out_prefix, | |||
max_prob_for_violation = 0.1 , |
|||
min_dist_for_violation = 1000000000.0 , |
|||
coarsen = False , |
|||
limit_to_chains = None , |
|||
exclude_chains = '' |
|||
) |
Create CMM files, one for each state, of all cross-links.
will draw in GREEN if non-violated in all states (or if only one state) will draw in PURPLE if non-violated only in a subset of states (draws nothing elsewhere) will draw in RED in ALL states if all violated (if only one state, you'll only see green and red)
out_prefix | Output xlink files prefix |
max_prob_for_violation | It's a violation if the probability is below this |
min_dist_for_violation | It's a violation if the min dist is above this |
coarsen | Use CA positions |
limit_to_chains | Try to visualize just these chains |
exclude_chains | Try to NOT visualize these chains |
Definition at line 833 of file crosslinking.py.