[IMP-users] coarse models with MCCG sampler : biased ?
To: Help and discussion for users of IMP <>
Subject: [IMP-users] coarse models with MCCG sampler : biased ?
From: Benjamin SCHWARZ <>
Date: Wed, 14 Sep 2011 17:33:54 +0200
Reply-to: Help and discussion for users of IMP <>
Hi guys,
I am currently interested in the coarse modeling of a small (10 subunits) system for which I have nothing but a few interactomics data as well as some (parts of) structures to help me to assess the size and number of beads to represent some of the subunits (details on my OS as well as on the modeled system can be found below.).
Thanks to the nup84 example, I was rapidly able to generate a large number of configurations and to visualize it (thanks for the MCCG sampler and the IMP.display.PymolWriter). I am however puzzled by the look of the obtained configurations : having represented my subunits as a succession of beads (I explicitly defined my restraints between successive balls), I expected my subunits to accommodate various types of shapes whereas they seem to have a tendency to adopt preferentially a straight aligned rod-like shape ( An example of a pymol file containing 1000 configurations can be downloaded for visualization at alnitak.u-strasbg.fr/~schwarz/Perso/1000-configs.pym).
I am now wondering if these configurations are really random generated, or if -as I fear- I have a bias in my sampling scheme, and how I can circumvent the trouble.
Here is one explanations I could find to this apparent behavior : I don't know how the excluded volume is handled by the MCCG sampler, but my guess is that the CG step helps a lot in "locally" optimizing this particular restraint. Now, if all balls are considered in the computation of derivatives (and not only those that really present a common intersection), it is intuitively quite probable that the system will ultimately push successive balls the furthest away from each other, hence favoring rod-like shapes upon epileptic-snake-like ones. Does this assumption makes any sense ?
I am running IMP kernel SVN 10177 with Boost.FileSystem on Mac OS
My system is composed of
- 10 subunits (divided in 30 beads to represent specific domains)
- 29 restraints to enforce pair connectivity between domains of a same subunits, or domains of different subunits for which a connectivity is known
- 2 ambiguous connections between several domains or subunits
- 1 excluded volume restraint between every 30 particles in the system.
--Ben (S)