[RNABuilder-help] Re: MMB RNA threading

[Samuel Flores]

The problem of residue 309 is with atom C5', or more precisely the angle formed with O5' and C4'.  it's a bit obtuse -- 153 degrees.  If it's made a little more acute (basically by moving C5'  a little in the negative x,y,z directions) the problem disappears.  

I am working on a fix for this issue but it's actually a bit subtle.  in the meantime moving the atoms in question a bit should help you work around the issue.

Sam


> I think I answered my own questions. I temporarily removed chains (B, C) from the last.1.pdb file. I will thread them individually and merge them later into a single pdb file. I assume coordinates will be retained.
> 
> Right now I loaded the entire rna sequene in the input file and MMB doesn't seem to be limited by 3k.
> 
> If this performs well I will cite you all over the place :)
> 
> uzay
> ________________________________________
> From: UGUR UZAY Sezen
> Sent: Thursday, February 14, 2013 10:27 AM
> To: Samuel Flores
> Subject: RE: MMB RNA threading
> 
> I managed to thread 3000nt so far. I had to delete another residue #1623. I have a few questions.
> 
> [1] MMB tutorial says 3000nt is the upper limit for threading. My sequence is 3388nt long. How do you suggest to get over this problem?
> [2] 3iz9 has 3 chains. As you see in the picture (attached) the two chains extend out without being threaded. Is there a way to thread multiple chains using the same command file?
> 
> So far so good!
> 
> uzay
> ________________________________________
> From: Samuel Flores [samuelfloresc at gmail.com]
> Sent: Wednesday, February 13, 2013 8:55 PM
> To: UGUR UZAY Sezen
> Subject: Re: MMB RNA threading
> 
> you can delete parts of your last.1.pdb to see which pieces caused the problem.
> 
> I would like to try to figure out why those residues are problematic -- in the meantime this way can help you at least do part of your model.
> 
> 
> 
> sam
> 
> 
>> I deleted 309 + 310 and voila! it worked until 1500's. Now I encounter the same error. How can you tell the problem bases in pdb file of 3iz9? I attached the command file for you to see.
>> 
>> I appreciate your input.
>> 
>> uzay
>> ________________________________________
>> From: Samuel Flores [samuelfloresc at gmail.com]
>> Sent: Tuesday, February 12, 2013 11:55 PM
>> To: UGUR UZAY Sezen
>> Subject: Re: MMB RNA threading
>> 
>> I'm at a conference and am having trouble troubleshooting this properly.  However I can tell you that if you delete residues 309 and 310 from last.1.pdb you remove the problem.  perhaps there is some geometric issue with these residues?
>> 
>> 
>> sam
>> 
>>> Great!
>>> 
>>> Here's the command file and the alignment. I commented out most of the sequences while playing around.
>>> 
>>> I appreciate your guidance.
>>> 
>>> Uzay
>>> 
>>> Uğur Uzay Sezen
>>> Post-Doctoral researcher
>>> Plant Genome Mapping Laboratory
>>> University of Georgia
>>> Center for Applied Genetic Technologies
>>> 111 Riverbend Rd., Athens, GA 30602
>>> uzay at uga.edu
>>> tel: 706-583-0170
>>> fax: 706-583-8120
>>> http://www.plantgenome.uga.edu/
>>> http://naturedocumentaries.org/
>>> 
>>> ________________________________________
>>> From: Samuel Flores [samuelfloresc at gmail.com]
>>> Sent: Tuesday, February 12, 2013 4:55 PM
>>> To: UGUR UZAY Sezen
>>> Cc: rnatoolbox-help at simtk.org
>>> Subject: Re: MMB RNA threading
>>> 
>>> Hi Uzay,
>>> 
>>> Happy to help.  Can you send me your command file?  I will take a look.
>>> 
>>> Sam
>>> 
>>>> You have been contacted by a member of Simtk.org.
>>>> 
>>>> From: Uzay Sezen
>>>> Email: uzay at uga.edu
>>>> 
>>>> The e-mail was sent via the "contact" link on Simtk.org by the above member. It is not screened by Simtk.org or its staff.
>>>> 
>>>> 
>>>> *** MESSAGE BEGINS ***
>>>> 
>>>> I successfully run MMB using the examples in tutorial. I want to thread Sorghum ribosomal RNA on wheat model (3IZ9) but can't get it to work even in short (less than 120nt) sequences in the command input file. Can anybody provide me feedback?
>>>> 
>>>> All the best,
>>>> 
>>>> Uzay
>>>> 
>>>> uzay at uga.edu
>>> 
>>> <commands.3iz9-threading.dat><Sbi26S-Tae28Salign.txt>
>> 
>> 
>> <commands.3iz9-threading.dat>
> 
> 
>