Hi,

I'm currently trying to calculate the structure of a covalently linked
homodimer using xplor v2.33 on Ubuntu 10.04 using only distance , talos
dihedrals, and h-bond restraints.  At some point during the calculation,
there's an error that appears (shown below) about the step size being too
large then too small.  I'm essentially using the same parameters as the
anneal.py found in the eginput>gb1_rdc example in the xplor package.  I
have also used the same anneal.py file for the monomeric protein and
everything worked smoothly.  The only difference between the monomeric and
dimeric protein is the number of restraints, which is doubled.

I also attached my logfile for more information if needed.  Any help would
be really appreciated!

Thank you,
Marie-Laurence Tremblay, PhD candidate
Department of Biochemistry and Molecular Biology
Dalhousie University
Phone:902-494-4812

[stdin](1): execfile('round25.py')
[stdin](1): execfile('round25.py')
 MAKINB: mode   3 found  15831 exclusions,       0 interactions(1-4) and
15831 GB exclusions
 MAKINB: mode   3 found      0 exclusions,       0 interactions(1-4)
and      0 GB exclusions
*--- Dynamics ---- step=      0 ---- time=         0 ---- delta_t=
0.0015999 --*
| E(kin)+E(poten)=     20989.468 E(kin)=      6974.370 temperature=
3457.795 |
|   E(poten)= 14015.0982135       grad=    67.5665364      ANGL=
2326.9185219 |
|       BOND=     9.1756895       CDIH=  2906.9832523      HBDA=
8.1348803 |
|       IMPR=  3231.5519843       RAMA=    -0.5304373       VDW=
3770.5403615 |
|        noe=  1762.3239610
                                  |
*------------------------------------------------------------------------------*
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: large timestep detected. Halving.
InternalDynamics::step: stepsize too small.
Traceback (most recent call last):
  File "<string>", line 2, in <module>
  File "/Applications/xplor-nih-2.33/python/trace.py", line 180, in run
    exec cmd in dict, dict
  File "<string>", line 1, in <module>
  File "round25.py", line 364, in <module>
    averagePotList=potList).run()
  File "/Applications/xplor-nih-2.33/python/simulationTools.py", line 360,
in run
    s.structLoopAction(s)
  File "round25.py", line 324, in calcOneStructure
    cool.run()
  File "/Applications/xplor-nih-2.33/python/simulationTools.py", line 1481,
in run
    s.runIVM(s.bathTemp)
  File "/Applications/xplor-nih-2.33/python/simulationTools.py", line 1534,
in runIVM
    s.ivm.run()
  File "/Applications/xplor-nih-2.33/python/ivm.py", line 513, in run
    (done,s.stepsize_) = PublicIVM.step(s,s.stepsize());
  File "/Applications/xplor-nih-2.33/python/wrappers/publicIVM.py", line
137, in step
    def step(self, *args, **kwargs): return _publicIVM.PublicIVM_step(self,
*args, **kwargs)
_publicIVM.IVMError: IVM error: stepsize too small
PyInterp::command: error executing: >execfile('round25.py')<
 HEAP:   maximum use= 38465598 current use= 24334970
_______________________________________________
Xplor-nih mailing list
Xplor-nih@cake.cit.nih.gov
http://cake.cit.nih.gov/mailman/listinfo/xplor-nih

Reply via email to