Hi,
Loop routine doesn't use coordinates for the loop residues from the pdb file anyway. All atoms in the residues selected for the loop modeling are aligned on the line that spans two ancoring residues before the loop modeling begins.
Good luck, Bozidar
On Thu, 25 Apr 2002 AA_Canutescu@fccc.edu wrote:
> Dear Bozidar, > > Thank you for help. Your hint(to put '' instead of '0') looks to work so far. > Regarding the fact that coordinates for the loop were zeroed, I did this in > order to get a de novo loop prediction. > > Adrian > > Quoting Bozidar Jerkovic bozidar@oboe.rockefeller.edu: > > > > > Dear Adrian, > > > > A simple answer to your question is: you specified a non-existing > > segment of the protein for the loop optimization: > > SELECTION_SEGMENT = '151:0' '154:0', SELECTION_STATUS = 'initialize' > > > > In your pdb file there is no '0' chain: > > ATOM 1165 CD2 LEU 150 44.348 42.655 43.046 1.00 15.23 C > > ATOM 1166 N GLY 151 0.000 0.000 0.000 1.00 16.38 N > > ATOM 1167 CA GLY 15 0.000 0.000 0.000 1.00 15.57 C > > > > This results in MODELLER selecting the whole protein to perform 'loop' > > routine on it: > > Number of all, selected real atoms : 1480 1480 > > > > Another thing I noticed is that you have compresses all the residues you > > are selecting into one point? Is there any specific reason for this? > > ATOM 1166 N GLY 151 0.000 0.000 0.000 1.00 16.38 N > > ATOM 1167 CA GLY 151 0.000 0.000 0.000 1.00 15.57 C > > ATOM 1168 C GLY 151 0.000 0.000 0.000 1.00 16.31 C > > > > MODELLER complains about this though, but it keeps on: > > serious non-bonded atom clash: 1143 1148 0.000 > > serious non-bonded atom clash: 1143 1149 0.000 > > serious non-bonded atom clash: 1143 1150 0.000 > > > > Energies are way off the charts: > > Initial value of energy before optimization : 477452.2400 > > Current energy : 476905.9000 > > Maximal energy, step : 525495.1200 > > 28 > > Minimal energy, step : 449273.2100 > > 81 > > Average energy, stand.dev. : 506870.8100 > > 15563.8500 > > > > What causes MODELLER to become a memory hog is the fact that you > > unintentionally tried to perform an ab initio calculation of the whole > > protein. I myself almost crashed our main Linux machine in the group, > > running your top file. A job used up 92% of the available memory at one > > point. > > > > After I corrected the .top file: > > SELECTION_SEGMENT = '151:' '154:', SELECTION_STATUS = 'initialize' > > > > Energy returned to normal: > > Initial value of energy before optimization : 27.4512 > > Current energy : 4.5270 > > Maximal energy, step : 67.8339 > > 34 > > Minimal energy, step : 25.2391 > > 363 > > Average energy, stand.dev. : 43.2928 > > 10.7857 > > > > ... and MODELLER produced a beautiful model ;-) > > > > Summa sumarum: > > One has to make sure that when s/he selects a segment of the protein, > > using 'res_ID:chain_ID' convention (which is updated for the 6v2 > > release), that there is an existing chain_ID there. This is actually > > quite a common mistake many users do: their pdb has a chain_ID, and they > > omit it in their top file (just opposite of what you did). > > > > I am still curious why you have compressed all the loop residues into a > > single point (why those residues have all the coordinates equal to > > zero)? It worked, because those resides were brought back into the > > protein in the model MODELLER calculated, which is obviously good, and > > their energy is also fine. If this is some empirical finding, and not > > missing coordinates, I would certainly like to learn more about it. > > > > ...and there is no memory leak, thankfully ;-) > > > > All the best and happy modeling, > > Bozidar > > > > > > > > -----Original Message----- > > From: AA_Canutescu@fccc.edu [mailto:AA_Canutescu@fccc.edu] > > Sent: Wednesday, April 17, 2002 6:49 PM > > To: Bozidar Jerkovic > > Subject: Re: MODELLER > > > > > > Dear Bozidar, > > > > > > > > Yes I am aware of this one. I think it may have to do with static vs. > > > dynamic linking. > > > > It's not obvious to me why could be a linking problem, but anyway it's > > interesting this happens only for certain input files. For the input > > data sets that generate this bug, the bug is repro(reproducible -all the > > time). There is no situation(as far as I tested) in which for the same > > input data to get the software working sometime and not working other > > time. It just doesn't work -all the time - for that input data. > > > > > > > What did you mean by "repro"? > > > > An error that is exactly reproducible under any circumstances. It's used > > in programers' jargon. > > > > > > The other problem, that is actually the most serious one is also repro. > > For certain input data sets, modeller starts allocating huge amounts of > > memory (about 1GB/min on my machine). This indicates a serious memory > > leak. I got my kernal crashed twice today because of this, so I have to > > wait for your patch before using again modeller. > > > > I will attach again an input data set that leads to this memory leak. > > > > Everything is run on RedHat 7.2 on x86 architecture. > > > > Thanks, > > Adrian > > > > > > > > > > Thanks for these reports. It's sometimes hard for us to catch > > > everything. > > > > > > Best, > > > Bozidar > > > > > > > > > On 17/4/02 5:33 PM, "AA_Canutescu@fccc.edu" AA_Canutescu@fccc.edu > > > wrote: > > > > > > > Dear Bozidar, > > > > > > > > I have another problem to report, besides the memory leak one. For > > > some > > > > reasons, I got in certain situations the following error: > > > > > > > > ? FORTRAN Runtime Error: > > > > ? Illegal character in numeric input > > > > ? READ(UNIT=INTERNAL,... > > > > > > > > I wasn't able to track down any cause, and the error message is not > > > very > > > > explanatory to me. I would appreciate if you could give me a hint > > > about this > > > > sitation. This situation is a repro one, like also the memory leak > > > one. I run > > > > everything > > > > > > > > > > > > > > > > Thank you, > > > > > > > > Adrian > > > > > > > > > > > > Quoting Bozidar Jerkovic bozidar@salilab.org: > > > > > > > >> > > > >> Hi, > > > >> > > > >> Loop is definitely supported ;-). I'll see what's going on. > > > >> > > > >> Best, > > > >> Bozidar > > > >> > > > >> > > > >> On 10/4/02 1:18 PM, "AA_Canutescu@fccc.edu" AA_Canutescu@fccc.edu > > > >> wrote: > > > >> > > > >>> Hi, > > > >>> > > > >>> I attached a tar.gz that contains top, log and pdb input file. In > > > >> case 'loop' > > > >>> routine is not anymore supported in mod6v1, what would be a get > > > around > > > >> for my > > > >>> problem? > > > >>> > > > >>> Thanks for your help. > > > >>> > > > >>> Adrian > > > >>> > > > >>> > > > >>> Quoting Bozidar Jerkovic bozidar@salilab.org: > > > >>> > > > >>>> > > > >>>> Hi, > > > >>>> > > > >>>> Can you please send me your .top and .log and all input files and > > > >> I'll > > > >>>> run > > > >>>> it here. I can't actually see which part of .top fails. I presume > > > >> it's > > > >>>> "CALL > > > >>>> ROUTINE = 'loop' # do loop modeling", but I don't wanna > > guess. > > > >>>> > > > >>>> Best, > > > >>>> Bozidar > > > >>>> > > > >>>> > > > >>>> > > > >>>> On 8/4/02 3:24 PM, "AA_Canutescu@fccc.edu" > > > >>>> AA_Canutescu@fccc.edu > > > >>>> wrote: > > > >>>> > > > >>>>> Dear Bozidar, > > > >>>>> > > > >>>>> > > > >>>>> I downloaded the last version of MODELLER and I am running into > > > >> some > > > >>>>> problems. For the time being, I use MODELLER mostly for loop > > > >> modeling. > > > >>>> Using > > > >>>>> the > > > >>>>> previous version(mod6a), I used to create a top file like this, > > > and > > > >>>> everything > > > >>>>> worked ok: > > > >>>>> > > > >>>>> INCLUDE > > > >>>>> SET OUTPUT_CONTROL = 1 0 0 1 0 # Set to 1 1 1 1 0 for > > > >> debugging. > > > >>>>> SET SEQUENCE = '1dysA_47-50' > > > >>>>> SET LOOP_MODEL = '1dysA_47-50.pdb' > > > >>>>> SET LOOP_MD_LEVEL = 'refine2' # refine3 - fastest > > > >>>>> # refine2 - medium > > > >>>>> # refine1 - most thorough SET > > > >>>>> LOOP_STARTING_MODEL = 1 SET LOOP_ENDING_MODEL = 3 # number of > > > >>>>> loop models > > > >>>>> SET RAND_SEED = -38901 # change for different starting models > > > >>>>> CALL ROUTINE = 'loop' # do loop modeling > > > >>>>> > > > >>>>> # Pick model residues that need to be refined. > > > >>>>> > > > >>>>> SUBROUTINE ROUTINE = 'select_loop_atoms' # Uncomment if you > > > >>>>> also > > > >> want > > > >>>> to > > > >>>>> optimize the environment: > > > >>>>> # SET SELECTION_SEARCH = 'SPHERE_SEGMENT', SPHERE_RADIUS = 6 > > > >>>>> PICK_ATOMS SELECTION_SEGMENT = '47:A' '50:A', SELECTION_STATUS = > > > >>>> 'initialize' > > > >>>>> RETURN > > > >>>>> END_SUBROUTINE > > > >>>>> > > > >>>>> > > > >>>>> With the newest version of MODELLER that I downloaded last > > > >>>>> week, > > > >>>> doesn't work > > > >>>>> anymore. I receive an error message in the end of.log file like > > > >> this: > > > >>>>> > > > >>>>> act10___551E> No such routine: > > > >>>>> recover____E> ERROR_STATUS >= STOP_ON_ERROR: 1 1 > > > >>>>> > > > >>>>> Dynamically allocated memory at finish [B,kB,MB]: > > > >>>> 28478873 > > > >>>>> 27811.400 27.160 > > > >>>>> Starting time : > > > >> 2002/04/08 > > > >>>>> 14:20:40.274 > > > >>>>> Closing time : > > > >> 2002/04/08 > > > >>>>> 14:20:51.421 > > > >>>>> Total CPU time [seconds] : > > > >> 11.00 > > > >>>>> > > > >>>>> > > > >>>>> Could you help me please regarding this matter? Could you also > > > tell > > > >> me > > > >>>> if > > > >>>>> there is any dependency needed(some libraries, fortran, etc)? I > > > >>>>> am > > > >>>> using > > > >>>>> RedHat > > > >>>>> 7.2 with kernal 2.4.7-10 > > > >>>>> > > > >>>>> > > > >>>>> > > > >>>>> Thanks, > > > >>>>> > > > >>>>> Adrian Canutescu > > > >>>> > > > >>>> > > > >>> > > > >> > > > >> > > > > > > > > > > > > > > >