<p>Sorry, one more doubt!<br>
If I decide to use PDB database, is there a way I can run some routines on the template structure(e.g. Shake coordinates by 1A) before it is considered for MRage pipeline after it is downloaded from PDB server?</p>
<p>Thanks,<br>
Kaushik</p>
<div class="gmail_quote">On Nov 20, 2014 11:39 PM, "Kaushik Hatti" <<a href="mailto:hskaushik@gmail.com">hskaushik@gmail.com</a>> wrote:<br type="attribution"><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr">Hi Gabor,<div><br></div><div>Actually the database is not split at chain alone but domain level too. Many multi-domain proteins are split into each structure file containing a single domain eg, 1xxxB_1, 1xxxB_2,... I specifically want to pick homologous structure from this custom database and not the entire PDB. I will try to hack this at my end and will share it here if I succeed.</div><div><br></div><div>Thanks,</div><div>Kaushik</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Nov 20, 2014 at 11:10 PM, Gabor Bunkoczi <span dir="ltr"><<a href="mailto:gb360@cam.ac.uk" target="_blank">gb360@cam.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Kashik,<br>
<br>
hmm, I am not sure I can help you here. MRage uses some general purpose code to fetch PDB entries (from the iotbx.fetch module), and this code is also utilised elsewhere. People would probably frown at me if I modified this without consulting them first, but fortunately, they are also on this mailing list, and will hopefully chime in at some point.<br>
<br>
However, I am wondering whether you could not solve this problem in a different way. In case 1vjhA.pdb happens to be chain A from 1vjh.pdb, you could just change the accession number (blast XML Hit_accession record) to 1vjh_A, and MRage will automatically fetch entry 1vjh and extract chain A. Can you not recompile the database to use these accession numbers? The domain information you can encode with the alignment, i.e. simply omit the residues that do not belong to the domain.<br>
<br>
Best wishes, Gabor<div><div><br>
<br>
<br>
On 2014-11-20 12:26, Kaushik Hatti wrote:<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>
Dear Gabor,<br>
<br>
Thanks so much for the suggestions. I built custom sequence database,<br>
did all-vs-all blast and stored output as XML files. Also understood<br>
how blast output.xml could be used for MRage "homology" scope.<br>
<br>
However, some of my PDBIDs are not 4 but 6 characters long. The<br>
filename includes chainID and domain information as well. Hence, when<br>
I run MRage I get this error:<br>
<br>
XXXXXXXX<br>
Error: Invalid PDB ID '1vjhA.pdb'. IDs must be exactly four<br>
alphanumeric characters, starting with a number from 1-9.<br>
This entry cannot be used: Invalid PDB ID '1vjhA.pdb'. IDs must be<br>
exactly four alphanumeric characters, starting with a number from 1-9.<br>
XXXXXXXX<br>
<br>
I have also placed 1vjhA.pdb under "vj" dir, 2i9yA.pdb under "i9"<br>
etc...<br>
<br>
Am I missing something? How can I overcome this?<br>
<br>
Thanks,<br>
Kaushik<br>
<br>
On Mon, Nov 17, 2014 at 5:34 PM, Gabor Bunkoczi <<a href="mailto:gb360@cam.ac.uk" target="_blank">gb360@cam.ac.uk</a>><br>
wrote:<br>
<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>
Hi Kaushik,<br>
<br>
yes, this is possible, although with some restrictions.<br>
<br>
If you have a local copy of the PDB, you can set the PDB_MIRROR_PDB<br>
environment variable to point to the root directory. This does not<br>
have to contain the whole PDB, as long as the entries you need to<br>
access are present. Entries in the root folder should be divided up<br>
to subfolders that are given by the second and third character of<br>
the accession ID. Only the actual structure files can be fetched in<br>
this way, but this is sufficient for MRage to work.<br>
<br>
As for your query about customising blast, you can build a blast<br>
database against whatever subset you like, and make sure this is the<br>
default. If there are multiple blast executables in the path, the<br>
first one will be picked up. However, it may actually be easier for<br>
you to run the blast search separately, request XML output, and feed<br>
this to MRage's "homology" scope.<br>
<br>
Let me know if you need more details!<br>
<br>
Best wishes, Gabor<br>
<br>
On 2014-11-17 07:47, Kaushik Hatti wrote:<br>
<br>
</div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div><div>
Hi,<br>
<br>
I am new to Phenix. Sorry if my question is naive.<br>
<br>
For a given mtz file, I want to pick templates from a local PDB<br>
database (its actually iTasser template library) and perform<br>
Molecular<br>
Replacement using MRage. I have already installed Blast,<br>
downloaded<br>
iTasser PDB library<br></div></div>
(<a href="http://zhanglab.ccmb.med.umich.edu/library/PDB.tar.bz2" target="_blank">http://zhanglab.ccmb.med.<u></u>umich.edu/library/PDB.tar.bz2</a> [1] [1])<span><br>
and<br>
extracted sequences from both the downloaded PDB template library<br>
set<br>
and target protein for which mtz is available.<br>
<br>
How should I go about building pipeline to perform MR using<br>
templates<br>
from locally downloaded select PDB database?<br>
<br>
Are there any more parameters available for parameter name:<br>
services<br>
under scop: search other than "local" and "ncbi"?<br>
<br>
Thanks in advance,<br>
Kaushik--<br>
<br>
People living deeply have no fear of death - Anais Nin<br>
Caution: I am still the dumbest person I have ever known :-)<br>
<br>
Links:<br>
------<br></span>
[1] <a href="http://zhanglab.ccmb.med.umich.edu/library/PDB.tar.bz2" target="_blank">http://zhanglab.ccmb.med.<u></u>umich.edu/library/PDB.tar.bz2</a> [1]<span><br>
<br>
______________________________<u></u>_________________<br>
phenixbb mailing list<br>
<a href="mailto:phenixbb@phenix-online.org" target="_blank">phenixbb@phenix-online.org</a><br>
</span><a href="http://phenix-online.org/mailman/listinfo/phenixbb" target="_blank">http://phenix-online.org/<u></u>mailman/listinfo/phenixbb</a> [2]<br>
</blockquote><span>
<br>
--<br>
##############################<u></u>####################<br>
<br>
Dr Gabor Bunkoczi<br>
<br>
Cambridge Institute for Medical Research<br>
Wellcome Trust/MRC Building<br>
Addenbrooke's Hospital<br>
Hills Road<br>
Cambridge CB2 0XY<br>
##############################<u></u>####################<br>
</span></blockquote><span>
<br>
--<br>
<br>
People living deeply have no fear of death - Anais Nin<br>
Caution: I am still the dumbest person I have ever known :-)<br>
<br>
<br></span><span>
Links:<br>
------<br>
[1] <a href="http://zhanglab.ccmb.med.umich.edu/library/PDB.tar.bz2" target="_blank">http://zhanglab.ccmb.med.<u></u>umich.edu/library/PDB.tar.bz2</a><br></span>
[2] <a href="http://phenix-online.org/mailman/listinfo/phenixbb" target="_blank">http://phenix-online.org/<u></u>mailman/listinfo/phenixbb</a><br>
</blockquote>
<br>
<br>
</blockquote></div><br><br clear="all"><div><br></div>-- <br><div>People living deeply have no fear of death - Anais Nin<br>Caution: I am still the dumbest person I have ever known :-)<br></div>
</div>
</blockquote></div>