2 """This script shows how to create a system with multiple copies of the same molecule.
3 We also create some cross-links which take into account the ambiguity.
4 The key to ambiguity is using the same molecule name for ambiguous copies.
5 That way when you perform Selection it automatically finds all relevant molecules.
29 st1 = s.create_state()
30 st2 = s.create_state()
35 m1A = st1.create_molecule(
'ProtA',sequence,chain_id=
'A')
36 m1A.add_representation(m1A,resolutions=[1])
37 m1B = m1A.create_clone(chain_id=
'B')
41 m1C = st1.create_molecule(
'ProtC',sequence,chain_id=
'C')
42 m1C.add_representation(m1C,resolutions=[1])
45 m2A = st2.create_molecule(
'ProtA',sequence,chain_id=
'A')
46 m2A.add_representation(m2A,resolutions=[1])
47 m2C = st2.create_molecule(
'ProtC',sequence,chain_id=
'C')
48 m2C.add_representation(m2C,resolutions=[1])
57 for mol
in (m1A,m1B,m1C,m2A,m2C):
58 dof.create_flexible_beads(mol,
60 dof.create_super_rigid_body(mol)
68 lines =
'''id,mol1,res1,mol2,res2,score
71 tf = tempfile.NamedTemporaryFile(delete=
False, mode=
'w')
77 kw.set_unique_id_key(
"id")
78 kw.set_protein1_key(
"mol1")
79 kw.set_protein2_key(
"mol2")
80 kw.set_residue1_key(
"res1")
81 kw.set_residue2_key(
"res2")
82 kw.set_id_score_key(
"score")
84 xldb.create_set_from_file(tf.name)
96 output_objects.append(xlr)
97 dof.get_nuisances_from_restraint(xlr)
101 for mol
in (m1A,m1B,m1C,m2A,m2C):
104 output_objects.append(cr)
109 output_objects.append(evr1)
112 output_objects.append(evr2)
129 monte_carlo_sample_objects=dof.get_movers(),
130 global_output_directory=
'ambiguity_output/',
131 output_objects=output_objects,
132 monte_carlo_steps=10,
Simplify creation of constraints and movers for an IMP Hierarchy.
Setup cross-link distance restraints from mass spectrometry data.
Restraints for keeping correct stereochemistry.
void show_with_representations(Hierarchy h, std::ostream &out=std::cout)
Traverse through the tree and show atom info, including representations.
Set of Python classes to create a multi-state, multi-resolution IMP hierarchy.
Strings setup_from_argv(const Strings &argv, std::string description, std::string positional_description, int num_positional)
GenericHierarchies get_leaves(Hierarchy mhd)
Get all the leaves of the bit of hierarchy.
Protocols for sampling structures and analyzing them.
Represent the root node of the global IMP.atom.Hierarchy.
This class is needed to convert the keywords from a generic database to the standard ones...
Class for storing model, its restraints, constraints, and particles.
void transform(XYZ a, const algebra::Transformation3D &tr)
Apply a transformation to the particle.
A decorator for a particle with x,y,z coordinates.
General purpose algebraic and geometric methods that are expected to be used by a wide variety of IMP...
Create a restraint between consecutive TempResidue objects or an entire PMI Molecule object...
Create movers and set up constraints for PMI objects.
A class to create an excluded volume restraint for a set of particles at a given resolution.
Restraints for handling cross-linking data.
Python classes to represent, score, sample and analyze models.
A macro to help setup and run replica exchange.
Functionality for loading, creating, manipulating and scoring atomic structures.
this class handles a cross-link dataset and do filtering operations, adding cross-links, merge datasets...
Support for the RMF file format for storing hierarchical molecular data and markup.