FW: two attached small, pdb files

Forwarded by list-owner

------ Forwarded Message From: Michael Grabe mgrabe@itsa.ucsf.edu Date: Mon, 25 Aug 2003 12:33:05 -0700 To: modeller-care@salilab.org Subject: two attached small, pdb files

Dear Modellers,

Let me describe a minor problem I have encountered with modeller 6V2 on Mac OSX,

What I do:

I take a small alpha helix (25 amino acids) and I copy it and translate it.

I then concatenate the two alpha helix pdb files, separating the first and the second alpha helix by a TER line.

Note, the original helix is missing a capping oxygen at the last residue.

Now I renumber the pdb amino acids so that they go from 1 to 50 instead of 1 to 25 and then 1 to 25 (see helices.ini).

I then start to work with this pdb file without any alignment file.

I run this modeller script on the two alpha helices:

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% INCLUDE

READ_MODEL FILE = 'helices.ini'

READ_TOPOLOGY FILE = '$(LIB)/top_heav.lib' READ_PARAMETERS FILE = '$(LIB)/par.lib'

GENERATE_TOPOLOGY SEQUENCE = 'temp'

WRITE_MODEL FILE = 'refine0.pdb' OPTIMIZE MAX_ITERATIONS = 100 WRITE_MODEL FILE = 'refine1.pdb' ENERGY %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

Now if I look at refine0.pdb generated from modeller, I have a problem.

As I mentioned before amino acids 25 and 50 are missing capping oxygens, but modeller puts them on for me and refine0.pdb has two more atoms than the initial helices.ini.

***Here is the point, the added oxygen to amino acid 25 shifts the labeling of all of the atoms and their positions in the second alpha helix. ***

Look at and compare atoms 201 to 206 in the attached files to see exactly what I mean. Once you see, it will be clear that the optimization that follows has a bit of a tough time since all of the atoms in the second helix have shifted positions.

Any hints about how to make this better?

Thanks, Michael

------------------------------------------------------------------------ -------------------------- Michael Grabe, Ph.D. Post-doctoral Fellow Howard Hughes Medical Institute University of California, San Francisco 533 Parnassus Ave. San Francisco, CA 94143 mgrabe@itsa.ucsf.edu tel: ++ 415.476.6585 http://itsa.ucsf.edu/~mgrabe

------ End of Forwarded Message