Problem when reconstruct the protein with know structure!
Hi all I divided the protein with know structure (bovine rhodopsin receptor) into three part: N-terminus region, transmembrane region and C-terminus region without change their initial coordinates, and I try to use the three parts to reconstruct the original structure of the protein. First, I aligned the three parts as the following way: (sorry, last letter I submited a wrong alignment,the following is I used) >P1;N-terminus structureX:N-terminus aaaaaaaaaa------------------------------* >P1;tm-region structureX:tm-region ----------bbbbbbbbbbbbbbbbbbbb----------* >P1;C-terminus structureX:C-terminus ------------------------------cccccccccc* >P1;original sequence:original aaaaaaaaaabbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbbb*
then I used the following script to generate the original structure of my protein: from modeller import * from modeller.automodel import * env = environ() a = automodel(env, alnfile='bhro_mult.ali', knowns=('bhro_n','bhro_tm','bhro_c'), sequence='bhro') a.starting_model = 1 a.ending_model =1 a.make()
However, the generated model was different from the original one, looking closely, I found the difference mainly occurred in the N-terminus and C-terminus regions, their coordinates have changed greatly when compared with the original structure. I was missing something important? Or, did it in a wrong way? I am new to modeller, any suggestions are appreciated!
Thanks in advance
yours Xin Liu
liu xin wrote: > I divided the protein with know structure (bovine rhodopsin receptor) > into three part: N-terminus region, transmembrane region and C-terminus > region without change their initial coordinates, and I try to use the > three parts to reconstruct the original structure of the protein. ... > However, the generated model was different from the original one, > looking closely, I found the difference mainly occurred in the > N-terminus and C-terminus regions, their coordinates have changed > greatly when compared with the original structure.
Modeller will try to build models that resemble the template whenever it is aligned with your target sequence, but these restraints are on internal coordinates, not Cartesian coordinates, so there could well be rotation or translation of your three segments. In particular, by breaking the sequence into three, you have prevented Modeller from using the information present in your original structure about the distances between these three segments.
Ben Webb, Modeller Caretaker
participants (2)
-
liu xin
-
Modeller Caretaker