December 2016 - modeller_usage

EXTREME BOND LENGTHS IN OPTIMIZED MODELS
by Mahesh Velusamy 19 Dec '16

19 Dec '16

Dear All, I have tried to do build single point mutant models (n=100) along with separate optimization step on the neighbor residues around the mutated residue.I can able to get 98 models without any issue out of 100 models but in some models, the bond lengths of longer residues(ARG,LYS,TYR) became extreme than original value. For example: IN ARG residue, the distance between (CZ & NH1) or (CZ & NH2) was extreme (3.4Å) than the original value (1.5Å ) which makes the bonds disappear. I really don't understand that what makes these issue and I use the default MD&SA parameters which used in mutate_model.py with Lenard-Jones potential. *I don't get this issue if I change the restraint from Lenard-Jones into soft-sphere potential. So, kindly give me some suggestion about the parameters that I can tune to resolve this issue. Thanking you in advance With Regards Mahesh V

1 0

Difference between dynamic_sphere and lennard-jones potential restraints ?
by Mahesh Velusamy 12 Dec '16

12 Dec '16

Dear All, This might be very basic question and also I have read this from the manual https://salilab.org/modeller/9v8/manual/node102.html but still I am curious to know the answer for the following questions: 1.what is the difference between these two restraints in terms (when to use) usage ? 2.Which one do we need to prefer with respect to mutate_model.py? 3.Why dynamic_sphere kept as False in this script https://salilab.org/modeller/wiki/Mutate%20model ? 4.Is it good idea if I use dynamic_sphere in mutate_model.py ? Thanking you in advance With Regards Mahesh V

2 1

ERROR IN MODELS ASSESSMENT
by Mahesh Velusamy 12 Dec '16

12 Dec '16

Dear All, This is to inform you that I have introduced new residue type called BGC which has a Glycosidic linkage between ASP and Beta-D-GLUCOSE as mentioned in the following protocol https://salilab.org/modeller/6v2/FAQ.html#2. I can access this newly introduced residue successfully throughout the modeling but after the models generation I get the following warning message for the single models assessment. >> Model assessment by DOPE potential iatmcls_286W> MODEL atom not classified: BGC:N BGC iatmcls_286W> MODEL atom not classified: BGC:CA BGC iatmcls_286W> MODEL atom not classified: BGC:CB BGC iatmcls_286W> MODEL atom not classified: BGC:CG BGC iatmcls_286W> MODEL atom not classified: BGC:OD1 BGC iatmcls_286W> MODEL atom not classified: BGC:OD2 BGC iatmcls_286W> MODEL atom not classified: BGC:C BGC iatmcls_286W> MODEL atom not classified: BGC:O BGC iatmcls_286W> MODEL atom not classified: BGC:C1 BGC iatmcls_286W> MODEL atom not classified: BGC:O5 BGC iatmcls_286W> MODEL atom not classified: BGC:C5 BGC iatmcls_286W> MODEL atom not classified: BGC:C2 BGC iatmcls_286W> MODEL atom not classified: BGC:O2 BGC iatmcls_286W> MODEL atom not classified: BGC:C3 BGC iatmcls_286W> MODEL atom not classified: BGC:O3 BGC iatmcls_286W> MODEL atom not classified: BGC:C4 BGC iatmcls_286W> MODEL atom not classified: BGC:O4 BGC iatmcls_286W> MODEL atom not classified: BGC:C6 BGC iatmcls_286W> MODEL atom not classified: BGC:O6 BGC iatmcls_286W> MODEL atom not classified: ALA:OXT ALA preppdf_453W> No fixed restraints selected; there may be some dynamic ones. preppdf_454W> Restraints file was probably not read; use restraints.append(). Here, I just want to know what does these statements implies "No fixed restraints selected; there may be some dynamic ones. Restraints file was probably not read; use restraints.append()" Even, I added all these "BGC" atoms in "/atmcls-mf.lib" file but still I am getting this warning. So, Kindly let me know that is it necessary to worry about this warning ? If not can you guide me how to append this "BGG new residue" restraints? Thanking you

2 1

Is it possible to control atom movements of specific residue during optimization?
by Mahesh Velusamy 07 Dec '16

07 Dec '16

Dear All, I have been doing separate optimization steps on selected region where i experienced that certain residue(ARG,LYS,TYR) atoms get distorted (like bonds between two atoms were distorted). So is their anyway to control the random deviation of certain residues atoms. Thanking you in advance With Regards Mahesh V

2 1

Modeling loop using template from PDB file
by Markus Heller 07 Dec '16

07 Dec '16

Hello, I have a structure (A) which is missing a stretch of a dozen or so residues. I also have the structure of the protein from a different organism which has said stretch of residues resolved (B). Now I'd like to take the coordinates of A and model the missing residues using the relevant coordinates (i.e. for these residues only) from B. How would I go about doing this? Thanks Markus

2 3

using crosslinking information as restrain
by ncuspywangxuedong 07 Dec '16

07 Dec '16

Hey all, I was using modeller in order to obtain a model satisfying crosslink restrains. Now I know the crosslinked sites, but I do not know how to add crosslink restrains. Could you tell me how to do it? It is best If you can give specific examples to me! thanks a lot!!!! Su Mu

2 1