Hi Foggy,

Here is another thing you can try:  

Ignore the chain B completely for now, just you would a ligand.  Run autobuild with rebuild_in_place=False, and then afterwards with rebuild_in_place=True to get the best model you can. Even without 10% of the model I'd expect you might get a bit better model than you have now that way.

Then focus on chain B. If you cannot see something fairly clearly in your map then autobuild won't be able to either.  So if there is clear density for chain B, autobuild should build it using the methods Kostya and I suggested before.  You can also use phenix.cut_out_density to just cut out the density for chain B and then you can build into that density with phenix.find_helices_strands or phenix.build_one_model, then put it back in the unit cell.

All the best,
Tom T




From: [email protected] [[email protected]] on behalf of Konstantin Kogan [[email protected]]
Sent: Sunday, September 14, 2014 12:18 AM
To: Zhang Foggy
Cc: [email protected]
Subject: Re: [phenixbb] [ccp4bb] Autobuild in specific area

Hi Foggy,

I had similar situation recently. I would suggest to build manually or in coot (place a helix) at least one helix that you see clearly. Then check that protein a and protein b (helix) are different chains in model pdb. Then provide the sequences of both chains. In this situation autobuild will certainly rebuild this helix with real sequence and will try to build more (don't forget to select 'build outside model' option). I would also increase the number of rebuilding cycles to something 20-30. In my case it took around 7 cycles until it started to put any new amino acid.

Kostya

On Sep 14, 2014, at 7:51 AM, Zhang Foggy wrote:

Hi, Tom, 

Thanks for your great suggestions. 
Here is my awkward situation: protein-A has 900 amino acids, while protein-B only 70. When I tired your first suggestion, I found that the program only built up the protein-A part. No extra model was built up at all.
Additionally, there is no published model for protein-B. Sequence alignment and homolog modeling indicate that there are about 40 aa forming two alpha-helix with the rest unknown. I can only see clear density on one alpha-helix at the moment, and the R/Rfree values are 0.37/0.43 (2.5 A, P1). Regarding to your second suggestion, I do not find this option though GUI yet, will try tomorrow. Do you have any other suggestions on my situation? Thank you in advance.

Best,

Foggy

2014-09-13 21:03 GMT-05:00 Terwilliger, Thomas Charles <[email protected]>:
Hi Foggy,

This is a situation that autobuild does not address very well yet.  Yes, you want rebuild_in_place=False.  And yes, if you do that it will remove your model and rebuild it. It will also try to build into the remainder of the map.

In practice, you already have pretty much the best that autobuild is going to produce.  Delete the protein A chains that autobuild built, and put in your original protein A. Now you have your original model plus whatever autobuild was able to build.

An alternative is:  select "consider_main_chain_list=True" .  This will take your input model A and at each step it will check to see if parts of that model are better than any other model that it has built. This can help for keeping the input model that you put in.

All the best,
Tom T


From: CCP4 bulletin board [[email protected]] on behalf of Zhang Foggy [[email protected]]
Sent: Saturday, September 13, 2014 10:40 AM
To: [email protected]
Subject: [ccp4bb] Autobuild in specific area

Hi, All,

I am trying to use Autobuild to build a protein-protein complex model. I obtained partial phase solution (protein-A) through MR, and can see some density of protein-B near protein-A. I would like to use Autobuild to build protein-B model into this area. I load data file, initial density map file, sequence file (both of protein A and B), and protein A model, and run Autobuild. However, The "rebuild-in-place" function confused me a lot, no matter I set it to True or False. When I set to True, it only rebuilt protein A model without adding any additional residues; while I set it to False, it deleted the input protein-A model, and rebuilt a new one... Could somebody tell me how to build the protein B model with keeping the input protein A model in Autobuild? Thanks.
Foggy

_______________________________________________
phenixbb mailing list
[email protected]
http://phenix-online.org/mailman/listinfo/phenixbb